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ABSTRACT 


Viscoelastic  Characterization  of  a Nonlinear, 

Glass/Epoxy  Composite  Including  the 
Effects  of  Damage.  (December  1974) 

Scott  Williams  Beckwith,  B.S.,  Texas  A&M  University. 

M.S.,  California  Institute  of  Technology 
Chairman  of  Advisory  Committee:  Dr.  R.  A.  Schapery 

4' 

Isothermal  creep  and  recovery  tests  were  conducted  on  an 
epoxy  resin  and  a glass  fiber-reinforced  composite  made  from  the 
same  bulk  resin.  The  glass/epoxy  which  was  studied  included  uni- 
directional and  laminated  (angle-ply)  composites  as  well  as  samples 
removed  from  a Minuteman  III  solid  rocket  motor  case.  The  creep 

and  recovery  tests  were  carried  out  at  a series  of  stress  levels 

■■s 

well  into  the  nonlinear  region  at  temperatures  oc  20,  75  and  140°’T 
for  several  fiber  angles.  Both  the  epoxy  and  glass/epoxy  were 
found  to  be  thermorheologically  complex  materials  with  a creep 
compliance  which  may  be  represented  by  a power  law  in  time. 

The  linear  viscoelastic  principal  creep  compliances  were 
determined  for  the  glass/epoxy  using  fourth-order  tensor  transfor- 
mations. Using  the  Halpin-Tsai  relationships  and  the  "rule  of  mix- 
tures", the  principal  creep  compliances  were  compared  with  those 
predicted  by  micromechanics.  The  experimental  results  were  found 
to  agree  very  well  with  the  Halpin-Tsai  model  except  at  the  highest 
temperatures  and  were  within  the  upper  and  lower  theoretical  bounds 


i v 

'’on  compliance.  , Lvcn  at  low  stress  levels  the  presence  ot  mu.ro- 
, crack  growth  was  found  to  produce  appreciable  softening  at  the 
v.  highest  temperatures. 

The  nonlinear  properties  were  found  to  depend  primarily  on  the 
stress  normal  to  the  fiber,  suggesting  a crock  opening' inode  as  the 
essential  mechanism  of  growth.  Multiple  cycles  of  creep  arid  recovery 
showed  a disproportionate  amount,  of  damage  during  the  first  cycle. 
Crack  growth  was  fuUmi  Lu  reduce  mure  rapidly  and  with  less  softening 
effect  in  the  laminated  (to)  composites,  probably  as  a result  of  the 
interfacial  harrier  between  the  layers.  In  general,  the  oil  angle 
composites  exhibit  considerable  softening  due  to  micro-crack  growth. 

.'/bending  tests  conducted  on  glass/epoxy  beam  and  plate  specimens 
brought  out  a strong  influence  of  the  strain  gradient. ^.Jt  was  found 
that  linear  theory  car.  be  used  for  most  of  the  useful  engineering 
range  of  application.  Nonlinear  theory  based  on  tensile  tests  of 
unidirection  specimens  predicts  considerably  more  reduction  in 
bending  stiffnc  ~ than  observed  experimentally 

Multiple  cycling  effects  un  the  glass/epoxy  composite  were 
found  to  be  more  sensitive  to  the  stress  normal  to  the  fibers  than 
to  the  shear  stress.  The  second-order  Lebesgue  norm  (L2 ) was 
found  to  approximately  characterize  multiple  cycling  effects  such  as 
seen  in  a solid  rocket  motor  case.  The  L?  norm  which,  at  a given 
time,  is  proportional  to  the  root  mean  square  value  of  the  stress  is, 
therefore,  proposed  as  a parameter  for  defining  the  damage  produced 
by  hydrotesti ng . 
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SECTION  I 
INTRODUCTION 

During  the  past  fifteen  years  there  has  been  an  increased 
emphasis  on  the  development  and  structural  application  of  composite 
materials.  In  general,  there  are  two  primary  divisions  of  compos- 
ite materials,  viz.,  particulate  and  fibrous  composites.  These 
classes  and  several  others  are  discussed  in  detail  by  Holliday  [1]; 
however,  this  study  Is  limited  to  fibrous  composites.  Typical 
fibrous  composites  in  current  use  within  the  aerospace  and  commer- 
cial fields  consist  of  composites  made  from  continuous,  parallel 
fibers  embedded  in  a matrix  material.  The  use  of  fiber-reinforced 
plastic  composites  such  as  glass/epoxy  and  graphite/epoxy  for 
structural  components  has  evolved  as  a result  of  increased  interest 
in  their  mechanical  properties  relative  to  conventional  materials. 
Of  particular  interest  are  their  high  strength-to-weight  and  high 
modulus-to-weight  characteristics . This  type  of  composite  is  dis- 
cussed in  detail  by  Ashton,  Halpin  and  Petit  [2],  Lubin  [3],  Tsai, 
Hal  pin  and  Pagano  [4]  and  [5].  This  interest  is  also  mirrored  in 
the  large  number  of  published  works  on  this  class  of  materials  as 
evidenced  in  the  abstracts  of  recent  literature  surveys  conducted 
specifically  on  the  mechanics  of  fiber-reinforced  plastics  by 


The  -'ormat  of  this  dissertation  follows  the  style  of  the 
Journal  of  Composi te  Materials . 


Beckwl th  et  al . [6-8] . 

Although  the  constituent  materials  are  usually  assumed  to  ex- 
hibit linear  elastic  b_eh_ayj_o_r  up  to  failure,  the  overall  fiber- 
reinforced  plastic  composite  exhibits  a significant  amount  of  time 
and  temperature  dependent  mechanical  behavior  in  many  service  envi- 
ronments. The  efficient  and  safe  design  and  utilization  of  these 
composites  demand  a good  understanding  of  their  viscoelastic  be- 
havior. The  surveys  by  Beckwith  et  al . [6-8]  point  out  that  while  a 
significant  amount  of  research  has  been  conducted  on  the  visco- 
elastic behavior  of  polymers  and  particulate  composite  polymers, 
there  have  been  only  a limited  number  of  investigations  of  the 
viscoelastic  behavior  of  plastic  composites. 

It  has  further  been  established  by  Lou  and  Schapery  [9],  Ashton 
[10],  Ha  1 pi n [11]  and  Schapery,  Beckwith  and  Conrad  [12],  that  the 
behavior  of  many  of  these  composites  is  not  linearly  viscoelastic 
except  at  small  strains,  often  well  below  the  design  limits 
normally  imposed  fcr  structural  applications.  More  recent  studies 
of  the  nonlinear  viscoelastic  behavior  of  unidirectional  glass/- 
epoxy  composites  [S]  seem  to  indicate  that  the  nonlinearity  can  be 
attributed  to  internal  crack  growth  within  the  matrix  (or  at  the 
fiber/matrix  interface)  and  viscous  flow  in  the  plastic  matrix. 

Composites  have  typically  found  their  way  into  several  areas 
involving  pressure  vessels  through  the  process  of  filament  winding 
techniques  [13-16].  The  early  emphasis  on  fiber-reinforced 


composites  can  be  attributed  to  their  application  In  solid  rocket 
motor  cases  in  the  late  1950‘s  and  early  1960's.  This  interest  is 
currently  at  a very  high  level  as  evidenced  by  the  use  of  some  of 
the  "advanced"  composites  such  as  graphite/epoxy  [15].  Part  of 
the  motor  case  design  problem  rests  in  the  complete  mechanical 
characterization  of  the  material  and  subsequent  application  after 
the  vessel  has  been  subjected  to  one  or  more  cycles  of  "hydrotesting" 
to  pressure  levels  about  10-25%  above  the  expected  operating  condi- 
tions [17].  This  technique  of  "non-destructive  testing"  (NDT)  has, 
in  face,  b^en  shown  to  cause  considerable  internal  damage  to  the 
composite  and  in  glass/epoxy  composites  creates  a condition  for 
further  damage  due  to  moisture  [14,  15]. 

The  importance  of  the  nonlinearity  due  to  crack  growth  is 
considered  significant,  particularly  during  the  first  few  cycles 
when  a considerable  amount  of  non-recoverable  deformation  occurs. 
There  is  also  an  associated  softening  of  the  composite  from  cycle- 
to-cycle,  which  should  be  accounted  for  in  determining  subsequent 
deformations  during  the  application  of  the  service  conditions. 

Constitutive  theory  for  fibrous  composite  materials  is  re- 
viewed in  Section  II.  In  addition  to  reviewing  the  basis  for  lin- 
ear elastic  constitutive  theory,  the  definition  of  linear  visco- 
elastic behavior  is  presented  and  the  effects  of  temperature  on  the 
overall  composite  response  are  discussed.  Schapery  [18]  has  re- 
cently reviewed  many  of  the  major  areas  of  viscoelastic  behavior 


normally  observed  in  fibrous  composite  materials.  Some  of  the  more 
important  areas  are  reviewed  here,  with  particular  emphasis  on 
permanent  damage  as  a result  of  crack  propagation  [12,  19,  20]. 
Finally,  the  prediction  of  effective  properties  using  some  of  the 
current  micromechanics  theories  [2,  4,  5]  are  discussed. 

The  experimental  program,  which  was  designed  to  provide  the 
mechanical  characterization  data  necessary  for  the  evaluation  of 
the  glass/epoxy  composite  material  currently  used  in  the  third 
stage  Minuteman  III  solid  rocket  motor  case  [17],  is  described  in 
Section  III.  The  program  consists  of  mechanical  characterization 
tests  on  unidirectional  and  laminated  glass/epoxy  composites  as 
well  as  the  matrix  material.  A small  number  of  tests  were  also 
conducted  on  samples  of  the  laminated  glass/epoxy  composite  taken 
from  the  actual  motor  case. 

Reported  in  Section  IV  are  the  experimental  results  from  a 
series  of  uniaxial  creep  and  recovery  tests  and  several  beam  tests 
designed  to  emphasize  the  effects  of  temperature,  stress  level  and 
load  history.  The  linear  viscoelastic  constitutive  properties  are 


established  for  the  various  materials  and  compared  with  several  of 
the  micromechanics  theories  used  to  predict  the  effective  proper- 
ties. The  observed  effects  of  material  nonlinearities  as  a result 
of  crack  propagation  are  discussed. 

The  effects  of  the  observed  nonlinear  material  behavior  of 
the  glass/epoxy  composite  on  the  design  and  analysis  of  typical 
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solid  rocket  motor  cases  are  reviewed  in  Section  V.  The  current 
approach  to  motor  case  design  and  the  effects  of  hydrotesting  and 
temperature  on  the  case  properties  are  discussed.  A method  of 
accounting  for  the  effects  of  multiple  loading  using  second-order 
Lebesgue  norms  is  proposed. 


f. 


SLOT l ON  II 

HL VIEW  OF  CONSTITUTIVE  THEORY  FOR 
FIBROUS  COMPOSITE  MATERIALS 

I introduction 

The  analysis  of  any  structure  or  body  in  terms  of  a resultant 
stress  (strain)  field  requires  the  satisfaction  of  a set  of  equil- 
ibrium and  strain-displacement  equations  [21-23].  This  set  of 
equations  is  independent  of  the  particular  material  makeup  but  the 
general  solution  depends  on  the  relationship  between  the  stress  and 
strain  tensors.  This  relationship  is  known  as  the  constitutive 
equation  and  in  the  theory  of  linear  anisotropic  elasticity  is 
-•el  erred  to  as  generalized  Hooke's  law  [22]  and  simply  as  Hooke's 
law  in  the  case  of  isotropic  materials.  The  application  of  Hooke's 
law  has  been  extensively  studied  during  the  past  century,  partic- 
ularly for  homogeneous,  isotropic  materials. 

Indeed,  several  excellent  treatments  of  the  subject  are  given 
in  Love  [21],  Sokol nikoff  [22]  and  Fung  [23]  to  mention  only  a few 
of  the  more  classic:  references.  These  references,  developed  pri- 
marily along  the  lines  ol  classical  elasticity  theory,  treat  Hooke's 
law  in  generalized  terms  initially  and  then  tend  to  concentrate  on 
the  more  specialized  case  uf  isotropic  materials. 

The  analysis  of  anisotropic  materials  is  obviously  more  com- 
plicated due  to  the  nature  of  the  constitutive  theory  and  therefore 


has  not  received  the  same  degree  of  attention  until  recent  years. 

The  works  of  Lekhnitskii  [24,  25]  and  Ambartsumyan  [26]  concentrate 
on  the  development  of  anisotropic  theory.  With  the  current  interest 
in  the  use  of  fiber-reinforced  composite  materials  at  an  all-time 
high,  their  works  have  served  as  primary  references  for  current 
analytical  developments  in  areas  where  fibrous  composites  are  being 
used  as  structural  materials. 

It  has,  in  fact,  only  been  within  the  last  few  years  that 
these  works  have  been  translated  and  presented  to  the  Western  world. 
Some  of  the  more  notable  presentations  of  their  works  in  the  eng- 
ineering terminology  were  done  by  Ashton  et  al . [2],  Ashton  and 
Whitney  [27],  Dong  et  al.  [28]  and  [5],  The  use  of  fibrous  com- 
posites by  the  aerospace  industry  for  structural  applications  has 
resulted  in  several  analytical  developments  and  a complete  reassess- 
ment of  the  overall  design  philosophy. 

Fibrous  composites  have  been  developed  along  two  different 
paths  classified  by  the  matrix  material  holding  the  load  carrying 
fibers;  namely,  metal  matrix  and  plastic  matrix.  The  metal  matrix 
composites  are  of  interest  primarily  for  their  higher  temperature 
capability,  although  they  often  possess  potentially  higher  strength 
and  stiffness  due  to  the  matrix  properties.  Plastic  matrix  compos- 
ites have  a lower  density  and  can  be  readily  fabricated  into  more 
complex  shapes  than  metal  matrix  composites,  thus  making  them  more 
attractive  for  filament-wound  pressure  vessels  and  similar 
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structures.  We  shall  be  more  concerned  in  this  dissertation  with 
the  plastic  (polymeric)  matrix  composites  which  exhibit  a signif- 
icant amount  of  time-dependent  mechanical  behavior  in  many  service 
envi ronments.  This  behavior,  termed  viscoelasticity,  has  only 
recently  been  investigated  with  regard  to  fibrous  composites 
[9,  11,  12,  18,  29]  although  considerable  attention  has  been  given 
to  particulate  composites  [20,  30,  31]  (e.g.  solid  propellants). 
Glass  fiber/epoxy  and  graphite  fiber/epoxy  composites  are  typical 
of  this  class  of  materials. 

The  application  of  linear,  anisotropic  viscoelastic  theory  to 
fibrous  composites  in  recent  years  is  evident  in  the  works  of 
Schapery  [18,  32],  Halpin  [11,  33]  and  Hashin  [34],  Unfortunately, 
as  a result  of  the  high  volume  fraction  of  fibers  in  the  composite 
and  the  gross  difference  in  both  physical  and  mechanical  properties 
of  the  constituents,  linear  theory  is  not  always  adequate  to  predict 
the  response  to  various  loading  situations.  In  many  cases  linear 
theory  at  least  provides  the  starting  point.  Nonlinear  viscoelastic 
behavior  is  reported  by  Schapery  [35]  and  Ashton  [10]  in  the 
United  States  and  by  Ogibalov  and  Tiuneeva  [36,  37],  Rabotnov  et  al . 
[38]  and  Martirosian  [39-4?]  in  the  Soviet  Union.  The  nonlinear 
behavior  may  be  "reversible"  or  "nonreversibl e" , the  latter  norm- 
ally considered  to  be  due  to  "microcracking".  "Reversible  non- 


linearities"  are  usually  attributed  to  secondary  bond  failure  within 
the  polymer  and  essentially  consist  of  polymer  molecules  sliding 
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past  one  another.  On  the  other  hand,  primary  bond  failure  in  which 
the  molecules  are  actually  torn  apart  leads  to  microstructural 
damage.  During  fabrication  there  are  a number  of  voids  and  flaws 
which  develop  as  a result  of  the  polymer's  curing  process  and  the 
difference  in  the  physical  properties  of  the  various  constituents. 

As  a result,  microstructural  damage  in  fibrous  composites  is 
thought  to  be  due  primarily  to  the  growth  of  these  small  flaws  or 
"microcracks"  [18]. 

Analysis  of  fibrous  composite  structures  requires  a careful 
consideration  of  the  constitutive  theory  of  the  material.  The 
decision  to  use  linear  or  nonlinear  anisotropic  viscoelastic  thee  y 
depends  on  the  particular  application  and  material.  The  develop- 
ment of  viscoelastic  theory  with  regard  to  fibrous  composites  will 
be  reviewed  in  this  section.  Particular  attention  will  be  given  to 
the  contribution  of  the  polymeric  matrix  to  the  overall  response 
and  the  influence  of  microstructural  damage  on  the  composite  con- 
stitutive relationship.  Thermal  effects  will  also  be  discussed 
since  they  strongly  affect  the  behavior  of  the  polymeric  phase. 

Viscoelastic  Behavior  and  Matrix  Constitutive  ’’'heory 

Typical  polymeric  matrix  materials  fall  into  several  classes 
of  epoxies,  phenolics,  polyesters,  etc.  depending  on  the  particular 
strength  characteristics  desired  and  many  other  design  considera- 
tions. All  of  the  materials  exhibit  some  degree  of  viscoelasticity 


and  are  generally  considered  to  be  of  a homogeneous,  isotropic 
nature.  Therefore,  In  addition  to  reviewing  matrix  constitutive 
theory  for  this  class  of  materials,  it  is  also  fitting  that  several 
ground  rules  be  defined  in  terms  of  what  is  meant  by  the  term 
"linear  viscoelasticity". 

The  behavior  of  viscoelastic  materials  falls  into  two  major 
divisions;  linear  and  nonlinear.  Farris  and  Schapery  [20]  recently 
conducted  an  extensive  review  of  the  entire  field  of  linear  and 
nonlinear  viscoelastic  constitutive  theory  and  the  reader  is  re- 
ferred to  this  excellent  review  article.  The  authors  noted  that 
although  an  extensive  amount  of  work,  both  theoretical  and  exper- 
imental, had  been  conducted  during  the  last  decade  toward  the  dev- 
elopment of  a nonlinear  viscoelastic  constitutive  theory,  the 
definition  of  linearity  was  often  incomplete  or  misunderstood.  As 
a starting  point  we  shall  review  linear  viscoelastic  theory  and 
proceed  from  there. 

Linear  Viscoelastic  Theory 

An  understanding  of  linear  viscoelastic  theory  is  important 
for  two  reasons.  First  of  all,  linear  theory  is  at  least  an 
approximation  to  real  behavior  and  one  should  be  familiar  with  its 
range  of  applicability.  Second,  the  exact  definition  of  material 
linearity  should  be  firmly  established  before  one  can  assess  the 
significance  of  nonlinear  behavior.  All  too  often  it  has  been 
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found  that  only  part  of  the  definition  of  linearity  has  been 
applied  [43,  44]  while  in  actuality,  the  material  is  indeed  non- 
linear. A definition  of  linearity  compatible  with  that  used  in 
current  literature  [18,  20]  will  be  given  in  order  to  provide  a 
precise  statement  of  linear  viscoelastic  behavior.  We  shall  neglect 
the  effects  of  large  strains  since  the  strains  in  fibrous  composites 
are  usually  below  3-5  percent. 

Statement  of  linearity.  Following  [18],  the  mechanical  be- 
havior of  a linear  viscoelastic  material  can  be  described  in  the 
form  of  a general  input-response  relationship.  By  applying  the 
definition  of  linearity  along  these  lines  we  treat  the  material  as 
a "black-box"  without  specifying  a priori  its  physical  makeup  and 
the  physical  significance  of  the  input  and  response  quantities. 

The  representation  of  a response  function,  R,  due  to  an  input 
function,  I , is  given  by 


R = R{I } (1) 

in  which  the  brackets  { } indicate  that  the  current  value  of  R is 
dependent  upon  the  entire  history  of  I,  and  not  just  its  current 
val ue. 

The  response  R is  said  to  be  linear,  and  the  material  conse- 
quently, linearly  viscoelastic,  if  and  only  if  it  satisfies  the 
property  of  homogeneity  and  obeys  the  superposition  principle  [18]. 
The  property  of  homogeneity,  or  proportional ity  is  given  by 
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R{cl } = cR{ I } ; c = constant  (2) 

while  the  superposition  principle  is  given  by 

RU  + iKl  = R(  I 1 + RU.  1 (3) 

d D d 0 

where  I and  I,  are  arbitrary  input  histories. 

d t) 

Although  the  use  of  Equation  (2)  alone  is  not  sufficient  to 
prove  linearity,  many  standard  characterization  methods  used  today 
use  only  this  criterion  to  establish  linearity.  However,  if  the 
material  obeys  Equation  (3),  it  can  be  shown  to  automatically  sat- 
isfy homogeneity  for  all  rational  values  of  c (including  c = 0). 
Thus,  the  primary  condition  for  a material  to  be  linear  is  that  it 
satisfies  superposition. 

Knowledge  that  a material  is  linear  is  sufficient  to  establish 
explicit  single-integral  expressions  connecting  responses  to  inputs 
[45].  The  response  R can  be  written  as  an  integral  containing  the 
actual  input  I and  the  response  of  one  preselected  input  history. 

For  use  in  the  characterization  of  viscoelastic  material'  it  is 
customary  and  convenient  to  use  the  Heaviside  unit  step  function, 

H,  as  the  preselected  input  history.  The  step  function  is  giver,  as 

«t  - to , jf;  *;  *:  w 

where  t is  the  current  time  and  t'  the  time  at  which  the  input  I is 
applied.  The  response  to  H is  usually  denoted  by  , and,  for  many 
applications  where  aging  considerations  are  not  important,  the 
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representation  of  becomes 

^ = RH(t  - f)  (5) 

As  a consequence  of  the  definition  of  the  step  function,  RH 

vanished  when  t < t ' In  as  much  as  aging  effects  are  excluded. 

If  the  Input  is  time  varying,  the  response  Is  found  by  considering 

the  input  to  be  the  limit  of  a sum  of  step-inputs.  An  integral 

relationship  between  input  and  response  Immediately  follows  in  the 

form  given  by  Vol terra  [45]  and  Pipkin  [46]  as 

t 

RafR H(t  " t')ffr,dt'  (6) 

— oo 

Equation  (6)  is  called  a hereditary  law  and  is  sometimes  referred 
to  as  the  convolution  Integral,  the  Duhamel  integral  or  the 
Boltzmann  superposition  Integral. 

The  superposition  integral  may  be  generalized  to  represent 
multiple  responses  due  to  several  inputs  when  the  response  of  the 
material  or  body  is  linear  with  respect  to  all  Inputs.  This  gener- 
alized form  of  the  response  Is  given  by 

‘I  W‘  - *'>31^'  <7> 

— ao 

in  which  R^  is  the  cumulative  response  due  to  all  inputs  and 
where  Ru  a is  the  unit  response  of  R when  only  I = H(t)  is 

nCtp  QL  P 

applied.  The  summation  convention  over  the  range  of  indices  a and 
if  will  be  implied  whenever  the  index  is  repeated  in  a term. 
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Isothermal  constitutive  equations.  The  general  linear  rela- 
tions between  the  stresses,  , (i , j * 1,  2,  3)  and  infinitesimal 
strains,  may  be  found  by  replacing  the  responses  and  inputs  in 
Equation  (7)  with  these  variaoles,  respectively,  hence 

t 

°u  - f W*  - <8> 


where  is  the  "effective  relaxation  modulus  tensor.  The 

spatial  dependence  of  o^.,  and  C.  are  implied  although  not 
shown  for  simplicity  in  the  notation  scheme,  and  the  stress,  a.., 

J 

and  strain,  , may  be  any  time-dependent  functions.  The  inverse 
relations,  in  terms  of  the  "cre^p  compliance  tensor",  S.^,  are 
given  by 


t . 

r'i j = J Sijks,(t  " ^ 

— go 

The  compliances  S..^  and  moduli  C.^,.  are  fourth-order 
tensors  which  possess  symmetry  with  respect  to  index  changes  of 
i with  j and  k wi th  i due  to  the  symmetry  of  the  stress  and  strain 
tensors  [23], 


Cijk.W  ■ cjm(t>  ’ CiV"1  - cjjtk(t)  (,0a) 

° Soi£k^''  O0b) 

Additional  symmetry  of  the  tensor  with  respect  to  an  interchange  of 


the  first  two  (1j)  and  second  two  (kt)  Indices,  i.e., 


ijkt^  = Ck*1%1^ 

(11a) 

ijkJi^  s SkJtij^ 

(lib) 

has  been  shown  by  Biot  [47]  based  on  the  thermodynamics  of  stable, 
irreversible  systems.  These  symmetries  considerably  reduce  the 
number  of  material  properties  required  to  completely  characterize 
linear  material  behavior. 

In  order  to  bring  out  the  significance  of  the  material  prop- 
erties in  Equation  (9),  consider  a generalized  creep  test  where  the 
time-dependent  input  stress  is  given  by 


°ij  " r1J  H<t>  02) 

where  all  the  a!,  are  constant.  Substitution  of  Equation  (12) 

into  the  constitutive  relation.  Equation  (9),  leads  to  the  function 
dH 

gx-i  which  is  the  Dirac  delta  function  S(t')  given  by 


6(f) 


0 , t'  f 0 
- , t'  = 0 


03) 


Consequently,  Equation  (9)  becomes 

t 

04. 

Jc 


Using  the  sifting  property  associated  with  the  Dirac  delta  function 
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Equation  (14)  reduces  to 


Lij  “ Sijk£^0ij  • 


05) 


If  the  test  is  a uniaxial  creep  test  then  a..  * and  we  find 

* J 

that 


cn  = Siin(t)o 


1 1 


(16) 


or,  with  single-index  notation  for  stress  and  strain  [22], 


i-  i = 5 n ( t ) a j 


07) 


where  is  used  to  denote  the  strain  in  the  x:  direction,  and  oj 
is  the  constant  value  of  stress  in  this  direction;  Srj(t)  is 
called  the  "uniaxial  creep  compliance".  When  the  material  is  iso- 
tropic it  is  often  common  to  rewrite  Equation  (17)  after  a change 
in  notation  t.o 


e = D(t)o 


(18) 


where  the  corresponding  terms  still  have  the  same  significance  as 
before . 

Noni so thermal  constitutive  equations.  It  has  generally  been 
found  that  viscoelastic  materials  exhibit  significant  temperature 
dependence.  The  influence  of  temperature  on  viscoelastic  behavior 
can  be  divided  into  reversible  and  irreversible  effects.  Irrever- 
sible effects  result  in  permanent  changes  such  as  primary  bond 
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rupture  and  weight  loss,  brought  about  by  thermal  degradation. 

These  effects  will  not  be  considered  within  this  discussion. 

Reversible  effects,  on  the  other  hand,  consist  of  thermal 
expansion  and  temperature  dependence  of  those  mechanical  properties 
which  appear  In  the  constitutive  equations,  l.e..  Equations  (8)  and 
(9).  Thermal  expansion  effects  will  not  be  treated  here;  a de- 
tailed analysis  of  these  effects  is  presented  by  Schapery  [48]  for 
anisotropic  materials  in  the  context  of  irreversible  thermodynamics. 
We  shall  concentrate  on  the  temperature  dependence  of  the  mechan- 
ical properties  in  our  discussion  of  reversible  thermal  effects. 

In  general,  these  effects  further  divide  material  behavior  into 
two  categories:  thermorheologically  simple  behavior  and  thermo- 
rheologically  complex  behavior. 

The  simplest  realistic  representation  of  material  response 
for  viscoelastic  materials  under  transient  temperatures  Is  that 
for  a "thermorhealogically  simple  material"  (TSM).  Morland  and 
Lee  [49]  originally  defined  a TSM  for  an  isotropic  material  and 
Schapery  [50]  extended  the  concept  to  anisotropic  materials  using 
irreversibly  thermodynamic  principles.  By  definition  [50]  for  a 
TSM,  the  constitutive  equation  is  given  by 

W5  ('9a) 


or,  equivalently. 


in nit 


W 


where  >:  is  the  so-called  "reduced  time"  defined  by 


C(t)  £ 


f dt' 

h ■ 


S'  = €'(f)  -3 


The  scalar  function  ay  reflects  the  influence  of  temperature  on 
internal  viscosity  and  is  quite  sensitive  to  temperature  changes. 
The  so-called  "temperature  shift-factor"  is  given  by 

aT  = aT[T(t')]  (21 

and  is  used  to  evaluate  the  reduced  times  r and  r, 1 . For  isothermal 


tests,  Equation  (20)  becomes 


; 


since  aT  is  constant  in  this  case.  By  rewriting  Equation  (22), 
(where  log  log10), 


log  f = log  t - log  ay  (23 

we  find  that  plots  of  isothermal  moduli,  or  compliances,  can  be 
shifted  horizontally  along  the  time  scale  with  the  magnitude  of  the 
shift  equal  to  log  ay.  It  is  interesting  to  note  that  Schwarzl 
and  Staverman  [51]  shifted  isothermal  data  in  this  manner  to  obtain 
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a single,  "master",  curve  as  a function  of  ? and  subsequently 
defined  a TSM  as  one  which  behaved  In  this  manner.  A reference 
temperature  Is  usually  selected  with  a^ defined  as  = 1. 

The  shift-factor  Is  usually  obtained  by  graphically 
shifting  the  data  horizontally  along  the  time  scale  as  noted  by 
Equation  (23).  However,  several  analytical  representations  have 
been  postulated  over  the  years  for  polymers.  When  T < T , an 
Arrhenius  type  of  temperature  dependence  Is  usually  observed 
[11,  52].  This  shift  factor  has  exponential  dependence  with  re- 
spect to  the  inverse  of  absolute  temperature  (1/T), 

1og  aT  = O WR  4 ' fR)  (24) 

where  aH  is  the  activation  energy  (per  mole),  R is  the  universal 
gas  constant,  and  TR  is  an  arbitrary  reference  temperature.  By 
plotting  log  a-j-  against  the  reciprocal  of  absolute  temperature, 

1/T,  a straight  line  results  and  the  activation  energy,  aH,  may  be 
found. 

When  T > Tg  the  so-called  WLF  equation  normally  applies  [52], 

"ci (T  - Tr) 

log  aT  = r2’  + T-"T^  (25) 

where  and  C2  are  constants.  The  form  of  the  equation  is  sim- 
ilar regardless  of  the  choice  of  TR,  only  the  values  of  Cj  and  C2 
change.  When  TR  is  chosen  to  be  approximately  50°C  above  the  T^, 
the  values  of  Cj  and  C2  have  been  found  to  be  universal  constants 


for  many  polymer  systems  [52,  53]  with  the  following  values, 
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C,  = 8.86 

C2  = 101.6  (°K) 


(26) 


Schapery  [32]  has  used  a power  law  form  of  aT  when  T > T , 


<T*  - V 


■T 

P 


(27) 


where  T and  )i  are  material  constants.  Again,  the  value  of  ID  is 

arbitrary  and  it  has  been  found  that  p assumes  a typical  range  of 

12-15.  T,  takes  on  values  several  degrees  below  the  T . This 
a g 

power  law  has  been  found  to  be  a good  approximation  of  experimental 
data  over  a wide  range  of  temperature  for  filled  polymers  [54]  and 
enables  the  reduced  time,  viz.,  Equation  (20),  to  be  evaluated 
analytically  for  constant  rates  of  change  of  temperature  whereas 
the  WLF  equation  does  not. 

The  inverse  representation  of  Equations  (19a)  and  (19b),  in 
terms  of  the  :reep  compliance  tensor,  is  given  by 

£,•  i " I s.~  (£  - (28a) 


C / ”i 


or,  equivalently. 


HI  * / *, 

*A 


3o 


“ ?,)W^  st' 


(28b) 
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In  both  cases,  Equations  (19)  and  (28),  we  have  neglected  the 
effects  of  thermal  expansion,  however,  these  terms  generally 
represent  only  a simple  addition  to  the  strain  tensor  as  shown 
In  [32],  Furthermore,  by  specializing  these  relations  to  isotropic 
materials  and  changing  notation  to  correspond  with  that  used  In 
Equation  (18)  we  find 

E(S  - S'^.dt'  (29) 

and 

D(e  - eOaf.dt*  (30) 

for  the  uniaxial  stress-strain  behavior  of  a TSM.  These  relation- 
ships correspond  to  the  constitutive  theory  proposed  originally  by 
Mori  and  and  Lee  [49]. 

There  Is  a significant  amount  of  published  data  to  verify 
Isothermal  curve  shifting  of  relaxation  moduli  and  creep  compli- 
ances from  several  sources  [52,  55].  Moehlenpah  et  al . [56] 
treated  isothermal  relaxation  moduli  for  an  epoxy  resin  to  form  a 
master  curve  In  terms  of  The  superposition  of  individual  iso- 
thermal responses  to  form  a master  curve  is  not  sufficient  to  estab- 
lish the  behavior  of  a TSM  although  most  characterization  methods 
tend  to  totally  ignore  this  fact.  Transient  temperature  tests  are 
necessary  to  completely  define  the  TSM  behavior.  Only  recently 
tests  of  this  type  have  been  conducted  by  Leeming  [57]  and 


i 


I 
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Farris  [58]  on  solid  propellants,  Johnson  et  al . [59]  and  Watkins 
[60]  on  epoxy  resins,  and  by  Kabelka  and  Vejchar  [61]  and  Schapery 
et  al . [12]  on  fibrous  composites..  With  few  exceptions  [12,  60] 
the  behavior  of  the  materials  characterized  using  transient  temp- 
erature  tests  have  been  TSM. 

The  data  of  Moehlenpah  et  al . [56],  Schapery  et  al . [12], 

Watkins  [60],  and  Sims  and  Hal  pi n [62]  defi ni tely  show  that  hori- 
zontal shifting  along  the  time  scale  is  not  entirely  sufficient  to 
form  a smooth,  well-defined  master  curve.'  However,  by  using  vert- 
ical and  horizontal  translations,  it  is  possible  to  superpose 
the  data  over  a large  time-temperature  range  [18]. 

Viscoelastic  materials  whose  temperature  dependence  cannot  be 
characterized  by  Equations  (29)  and  (30)  for  isotropic  materials  ^ 

and  Equations  (19)  and  (28)  for  general,  multiaxial  behavior  of 
anisotropic  materials  are  defined  as  "thermorheologically  complex 
materials"  (TCM) . Some  materials  are  composed  of  several  constitu- 
ents or  phases,  each  behaving  as  a TSM  with  different  a^  shift 
factors.  This  particular  type  of  TCM  has  been  designated  by 
Schapery  [18]  as  TCM-1.  The  behavior  of  materials  in  this  class 
has  been  studied  only  under  isothermal  conditions  which,  as  we 

[ 

noted  earlier,  is  net  sufficient  to  completely  verify  behavior  | 

h 

under  transient  temperature  conditions  [33,  63]. 

Schapery  [18]  has  also  defined  a second  class,  TCM-2,  composed 
of  materials  which  by  definition  satisfy  the  uniaxial  stress-strai n 
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behavior  under  isothermal  or  transient  temperatures  as  qiven  by 

t 

d 


= ' °0°  */  iD(E  - E'ljf-  (f^df 


(31) 


where  5 and  £'  are  given  by  Equation  (20),  a is  any  arbitrary 
stress  input,  and  DQ  = DQ(T)  is  the  initial  value  of  the  creep 
compliance.  For  the  case  of  an  isothermal  creep  test,  the  creep 
compliance  becomes 


0 ■ Do<T>  + 0f} 


(32) 


where  £ = — . Most  data  reduction  techniques  follow  a normalization 
aT 

condition  by  assuming  aG  = a^  = 1 at  some  arbitrary  reference 
temperature,  IR,  which  results  in 


AD(t)  = D(t,TR)  - Do(Tr) 


(33) 


where  AD(t)  is  the  transient  component  of  the  compliance  at  temp- 
erature Tr. 

In  order  to  determine  aT  and  aG  from  experimental  data  at 
different  temperatures  Schapery  [18]  rewrote  Equation  (32)  in  the 
form 

log  [D  - Dq]  = log  AD  - log  a^,. 


A plot  of  log  [D  - D ] against  log  t for  isothermal  tests  at  various 
temperatures,  T,  will  be  Identical  to  that  at  TR  except  for  a rigid 
horizontal  translation  of  log  aT  and  a similar  vertical  translation 
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of  log  a^.  The  final  result  is,  again,  a master  curve  similar  to 
that  obtained  for  a TSM. 

A special  case  of  Equation  (31)  occurs  when  the  temperature 


dependence  of  DQ  is  given  by 


which  results  in 


, .W 

-o  ■ "s^rrr 


(36a) 


where 


D(C)  Dq(Tr)  + &D(0  . 


(36b) 


As  a result  the  constitutive  relation  given  in  Equation  (31)  becomes 


°(t  - C)g|r  (~-)dt‘ 


which  is  the  same  as  Equation  (30)  except  for  ?G. 

Superposition  procedures  or  normalization  techniques,  in  add- 
ition to  that  of  Schapery’s  based  on  Equation  (34),  have  been 
suggested  for  accounting  for  the  behavior  of  TCM  [64-69],  McCrum 
and  Pogany  [70]  have  reviewed  several  of  these  procedures  and  com- 
pared the  master  curve  predictions  for  an  epoxy  resin  over  a temp- 
erature range  which  Included  the  glass  transition  temperature,  T . 
Four  of  these  different  techniques  are  given  in  the  following 


discussion. 


1.  Tobol sky-Ferry  procedure  [65,  65]: 


D0(T) 


_ Do<V 


aG(T)  9 TT 


- L_£_ 
R PR 


(3 


where  p and  pR  are  the  cities  at  temperatures  T 
and  Tj^,  respectively. 

2.  Ferry-Fi tzgerald  procedure  [68]: 


Do<T>  9 Do<V 


ar(T)  = L-H- 
G Tr  Pr 


3.  Ke  procedure  [66]: 


D_(TR) 

Do(T)  2 -brr 


G' 

Mtd> 


a m ~ ov  ff' 
aG  T = ~dJTT 


4.  McCrum-Morris  procedure  [67]: 


D0(T} 


aQ(T) 


D(0,T) 

D»<TR> 


0„(tr) 


or 


(3 


(4 


(< 


The  Tobol  sky- Ferry  and  Ferry-Fitzgerald  procedures  are  based 
on  the  kinetic  theory  of  rubber  elasticity;  the  Ke  procedure  uses 
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the  initial  compliance  to  define  the  temperature  dependence  of  a^, 
and  the  McCrum-Morrl s procedure  contains  the  first  three  as  special 
cases  by  allowing  for  arbitrary  variations  of  the  initial  and  long- 
term compliances  with  temperature.  The  most  successful  of  these 
procedures  is  the  McCrum-Morris  normalization  although  the  specific 
reduction  technique  given  by  McCrum  and  Pogany  [70]  requires  know- 
ledge of  the  initial  and  long-term  compliances  which  are  not  always 
available  because  of  experimental  limitations.  Another  form  of 
aG  proposed  by  Schapery  and  Martin  [69]  is  based  on  kinetic  theory 
and  is  given  as 

y In(j-) 

a6(T)  - '-Trr5 — <42> 

R _ i 

e 

where  y is  a factor  related  to  pressure,  volume  and  excess  molar 
energy  and  In  = logg. 

Graphical  shift-methods  [18]  appear  to  be  the  easiest  proce- 
dure for  providing  the  best  "average"  material  properties  in  the 
absence  of  the  limiting  values  of  compliance.  The  general  techni- 
que involves  the  use  of  Equation  (32)  by  making  preliminary  esti- 
mates of  Dq  and  then  forming  master  curves  of  AO(^) . By  using  a 
smoothing  technique,  the  best  master  curve  may  be  found  by  adjusting 
the  values  of  the  initial  estimate  of  D . 

■ Analytical  representation  of  time-dependent  properties.  No 

mention  has  been  made  up  to  this  point  about  the  representation  of 

i. 

A. 

f1 
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the  relaxation  and  creep  functions  in  Equations  (8)  and  (9). 
According  to  thermodynamic  theory  [47],  they  have  the  following 
time-dependence, 


Uijk£^^ 


C1  jki,  + 


(s)  -t/ps 

M e 


S1jkt^  = S?jk£  + Si  jk{.L  + XvS 


(s) 
1 jkn 


(1  - e"t/Ts) 


(43a) 

(43b) 


where  the  range  of  the  summation  index,  s,  depends  on  the  particular 

material,  p are  relaxation  times  and  x are  retardation  times.  The 

constants  C”.^  are  elastic  moduli  which  produce  long-term  stress 

o 

response  to  strains  whereas  the  constants  are  elastic  compli- 

ances which  produce  the  short-term  strain  response  to  stresses.  The 
exponential  coefficients  and  sj ^ and  the  corresponding  con- 

stants ps  and  ts  define  the  time  and  rate-dependence  of  the  material. 
The  coefficients  lead  to  steady-flow  under  constant  stress 

such  as  might  be  found  in  uncrossl inked  polymer  systems. 

Thermodynamics  gives  us  the  general  form,  Equations  (43a)  and 
(43b),  of  the  material  properties.  Although  there  may  be  a large 
number  of  time  constants  as  a result  of  the  polymer's  molecular 
configuration,  creep  and  relaxation  functions  are  often  approxi- 
mated by  exponential  series  [71]  consisting  of  relatively  few  terms 
(typically  ten  to  twenty). 

The  modified  power  law  [72] 


o 


(4 


E(t)  = Em  + 


0 +J>" 


v.nere  Em,  Eq,  tq  and  n are  Independent  of  time,  often  provides  an 


excellent  approximation  to  polymer  behavior  above  their  T . Also, 


when  t/t0  >>  1 then  Equation  (44)  reduces  to 


E(t)  ■ E. 


* <E0  ' E, 


)(t  )-n 

i ' T # 

0 


(4 


and  further  to 


E(t)  = Eq(”  )“" 
o 


(4 


when  Eu  <<  Eq.  A similar  form  for  the  compliance  is  [72] 


(D  - Dj  . „ 

D(t)  = D + ■■  , -0- £ )n 


(i  + i )n  To 
To 


where  DQ,  0^,  t0,  and  n are  material  properties.  If  t/TQ  <<  1 
then  Equation  (47)  reduced  to 

D(t)  = D + (D„  - D )(i  )"  (' 

0 

and  further  to 

D(t)  = Dm(~-  )n  (- 

o 

when  D >>  D„. 

CO  0 

The  forms  of  Equations  (48)  and  (49)  have  often  been  found  in 
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the  literature  as 

DC t)  - D0  + Djtn  (50) 

and 

D(t)  = Djtn  (51) 

respectively.  These  particular  forms  have  been  found  to  represent 
the  behavior  of  many  rigid  plastics  with  and  without  reinforcement 
[73-76], 

Nonlinear  Viscoelastic  Behavior 

We  have  already  defined  the  conditions  necessary  for  linearity, 
viz.  Equations  (8)  and  (9),  and  have  further  commented  on  the  use- 
fulness of  linear  viscoelasticity.  However,  the  behavior  of  most 
materials  is  generally  nonlinear  prior  to  complete  fracture  and, 
for  some  materials,  this  nonlinearity  exists  even  at  small  stress 
levels  well  within  the  design  range  of  structural  applications. 

There  are  many  theoretical  works  on  constitutive  theory,  both  lin- 
ear and  nonlinear,  such  as  those  of  Green,  Rivlin  and  Spencer 
[77,  78],  Coleman  and  Noll  [79,  80],  Truesdell  [81,  82],  Wang 
[83,  84],  Volterra  [45],  Pipkin  [85,  86],  Hermann  [87],  Lianis  and 
co-workers  [88,  89],  Schapery  [90,  91],  Huang  and  Lee  [92],  Dong 
[93,  94],  Williams  [72],  Tobolsky  [65,  95],  Alfrey  [96],  Ferry  [52] 
and  others  [97-107],  The  development  of  nonlinear  constitutive 


theory  has  also  resulted  In  several  applications  to  various  mater- 
ials such  as  reported  in  the  works  of  Findley,  Lai  and  Onaran 
[108-113],  Gottenberg  et  al.  [114],  Lockett  [115,  116],  Valanis 
and  Landel  [117],  Lee  and  Huang  [92],  Schapery  [35,  91],  Freudenthal 
[118],  Ward  [119-121]  and  others  [122-127],  It  would  be  a formid- 
able task  to  comment  on  all  of  these  works;  therefore,  in  the 
sections  to  follow,  we  shall  comment  only  on  observed  deviations 
from  linear  behavior  and  discuss  some  of  the  more  general  approaches 
taken  to  describe  nonlinear  viscoelastic  behavior. 

Equations  (2)  and  (3),  homogeneity  and  superposition,  respec- 
tively, completely  define  the  mathematical  rules  of  linear  material 
response.  Linear-  viscoelastic  behavior  is  achieved  through  the  use 
of  these  relationships  to  derive  Equations  [8)  and  (9).  Nonlinear- 
ity is  defined  only  as  the  failure  to  satisfy  one  or  both  of  these 
rules.  It  is  therefore  possible  for  a material  to  satisfy  che 
homogeneity  rule  and  fail  the  superposition  rule  [128,  129]. 

Multi  pie- Integral  theories . Because  of  its  generality,  the 
various  approaches  to  multiple-integral  theory  have  received  wide- 
spread attention  in  recent  years.  These  mathematical  theories  are 
probably  general  enough  to  account  for  almost  all  types  of  observed 
nonlinearities,  but  they  are  very  impractical  with  strong  nonlinear- 
ities.  This  type  of  representation  was  developed  by  Green,  Rivlin 
and  Sperc°r  [77.  130]  and  has  been  applied  Ly  several  investigators 
over  the  years  [110,  119,  124,  131], 
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For  the  case  of  uniaxial  loading,  the  multiple-integral  repre- 
sentation takes  the  following  form. 


e~J  Kl(t-Tl)-gn—  dr1  + 


f C <Mt.)  do(x„) 

j I K^t-T^t-T,)-^ dtldT2  + 


/// 


do(-r  ) do(r  ) d a(x  ) 

K3U-Tpt-T2,t-T3)^17--1^—  + 


t t t t 

/'// 

where  K , K2,  K3,  ...  are  the  kernal  functions  containing  time  (t) 
and  material  constants,  and  tj,  t2,  t3,  ...  are  the  dummy  time 
variables  of  integration.  The  corresponding  creep  relation 
(constant  stress)  takes  the  form 

e - K^tjcr  + K2(t,t)a2  + :<3(t,t,t)o?  + ...  (53) 

For  a linear  viscoelastic  representation  the  higher  order  kernals 
K?  and  K3  do  not  appear;  however  In  published  applications  of  this 
theory,  the  series  is  usually  truncated  at  three  irtegrals  for  the 
representation  of  nonlinear  behavior. 
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In  order  to  define  the  material  behavior  one  may  either  solve 
the  relations  by  numerical  techniques  or  assume  some  form  of  the 
kernal  functions  K j , Kz  and  K3.  Findley  and  Onaran  [110]  used  a 
product  form  suggested  by  Nakada  [107]  where  the  kerrial s take  the 
form 

Ki  = u^t-Tj)  (54) 

K*  = a2(t-T1)1/2  (t-x2)!/2 

«3  = a3(t-Ti)  ( t~T  2 ) ^ 3 ( t-T  3 ) 

where  uj , a2  and  a3  are  material  constants.  Nolte  and  Findley 
[132]  have  also  used  creep  kernal  functions  of  the  form 

Kl  = «l  + Mt-r  j )n  (55) 

K?  - a-/  + B2(t-Tl)n(t-T2)n 

K3  = a3  + e3(t-T1)n(t-r?)n(t-T3)n 

where  a j , a2,  a3,  s a , 83  and  n are  again  material  constants. 

A more  general  form  of  Equation  (55)  was  recently  proposed  by 
Smart  and  Williams  [133]  where  the  kernal  functions  may  be  written 

Ki  = aj ( t-T  i ) ^ (56) 

^2  = a?.(t-Ti)P+nl(t-T?)  P 

K3  = ^{t-r1)qfn(t-r2)r(t-T3)S 
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with  q + r + s a 0 and  where  alt  a2,  <*3  and  1 , m,  p,  q,  r and  s are 
the  material  constants. 

Without  going  further  Into  multiple-integral  representation  of 
nonlinear  viscoelastic  materials,  one  must  be  Impressed  at  this 
point  with  the  complex  nature  of  such  an  approach.  In  fact,  the 
representation  of  nonlinear  behavior  often  requires  more  than 
three  integrals,  particularly  if  a strong  nonlinearity  exists. 

The  solution  of  boundary  value  problems  and  the  inversion  of  even 
the  simplest  formulation  becomes  very  involved  and  complex.  In 
addition,  the  determination  of  the  material  functions  requires  an 
experimental  program  consisting  of  many  multiple  step  loading 
sequences.  Generally  speaking,  the  multiple-integral  approach  has 
not  found  Its  way  into  serious  engineering  consideration  with  regard 
to  material  characterization  and  subsequent  structural  applications 
of  composites.  As  a result  we  shall  not  consider  this  approach  any 
further. 


Single-integral  theories.  Several  methods  exi'.t  which  fall 
into  a class  called  single-integral  theories  of  no  d 1 near  visco- 
elasticity. The  modified  superposition  principle  (MSP)  was  first 
suggested  by  leaderman  [134]  who  observed  that  the  creep  behavior 
of  certain  fibers  and  plastics  could  be  separated  into  time  and 
stress-dependent  parts  so  that  the  creep  response  to  a stress,  0, 
could  be  written  as 


t 

DQa  + J AD(t-T)^dt 


(57) 


where  DQ  is  the  initial  value  of  the  compliance,  aD  is  the  trans- 
ient component  of  the  compliance  and  g is  a nonlinear  function  of 
the  stress,  a.  This  theory  has  met  with  varying  degrees  of  success, 
the  accuracy  depending  on  the  particular  material  as  well  as  on  the 
type  of  loading.  The  MSP  cannot  represent  a material  which  posses  - 
es  permanent  memory  effects  which  exist  as  a result  of  permanent 
internal  changes  such  as  microcracking . 

Schapery  [90,  91,  135]  used  thermodynamic  theory  to  develop  a 
nonlinear  constitutive  theroy  which  has  been  successfully  applied 
to  several  ma'erials  [9,  35].  When  stress  is  treated  as  the  inde- 
pendent state  variable,  then  the  theory  can  be  written  as 

t 


/dg2a 

)-g^-dT 


(58) 


i:! 
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where  DQ  and  aD  are  the  initial  and  transient  components  of  the 
1 inear  viscoelastic  creep  compliance,  t p is  the  so-called  reduced- 
t*;.,e  defined  by 


* - *«>  s f 5JFF 

’o 


IT 


and 


■/> 


dt ' 

RFTT 


(59a) 


(59b) 
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and  the  material  properties  gQ,  gl4  g2  and  ao  are  functions  of 

stress.  When  the  applied  stress  is  sufficiently  small, 

g0  8 gj  B g,  B aff  * 1 , and  Equation  (60}  reduces  to  the  familiar 

Boltzmann  superposition  Integral  for  linear  viscoelastic  behavior, 

t 

e = Dqo  + £ AD(t-r)^-  dt.  (60) 

The  MSP,  l.e.,  Equation  (57),  Is  obtained  from  Equation  (58)  by 
setting  gQ  = 9:  s aa  * 1 and  allowing  all  of  the  nonlinearity  to  be 
contained  in  g2.  The  stress-dependent  properties  have  specific 
thermodynamic  significance;  changes  in  gQ,  g:  and  g2  reflect  third 
and  higher  order  dependence  of  the  Gibb's  free  energy  on  the 
applied  stress,  and  aQ  arises  from  similar  high-order  effects  in 
both  entropy  production  and  free  energy. 

A constitutive  equation  with  strain-dependent  properties  has 
also  been  developed  [90]  and  can  be  put  In  the  form 

f dh  e 

a = heEee  + h1  / AE(P-p')-^-dT  (61) 

0 

with  the  reduced-time,  p,  defined  as 


and 


= p(t)  = 


t 

f dtV 

J ^ 
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= P*(t)  = 


•/ 


dt1 

aeLe C t ’ ) J 


(62a) 


(62b) 
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The  strain-dependent  properties  are  hg,  nl , h2  and  a£;  variations 
of  the  first  three  are  due  to  third  and  higher  order  strain  effects 
in  the  Helmholtz  free  energy,  while  changes  in  ac  arise  from  similar 
strong  strain  influences  in  both  entropy  production  and  free  energy. 

Schapery's  theory  involves  only  a limited  amount  of  experi- 
mental testing  in  addition  to  that  normally  required  by  linear 
theory.  The  determination  of  the  nonlinear  properties  does  not. 
represent  much  in  the  way  of  additional  work,  especially  considering 
the  improvement  over  MSP  when  significant  nonlinearity  exists. 

Consider  the  case  of  a creep  and  recovery  test  as  shown  in 
Figure  1,  where  the  stress  input  is  given  by 

o » oQ[H(t)  - H (t  - f)].  (63) 

Equation  (58)  yields  the  creep  strain 

= 9uVo + WD(r-}v  0 " 1 •'  r <64) 

a 

and  the  recovery  strain 

= 9-[^D(ip)  - - if‘)]a0>  t > t'  (65) 

where  in  Equation  (65) 

r = ~ ; if  = ~ + t - t\  (66) 

a u 

We  have  already  noted,  e.g.  Equation  (50).,  that  many  metals 
and  plastics  obey  a power  law  in  time  where  the  transient  component 
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may  be  represented  as 


AD(<p)  D^n  (67 

where  Dj  and  n are  independent  of  stress  level  and  time.  Rewritlny 
Equation  (64)  using  Equation  (67)  we  find 


n + 9l92  n 

e - goDo°o  + 7n  Dit  ao 


(68 


and  also  by  a similar  substitution  into  Equation  (65), 


Ac 


= --  CO  + a x)n  - (a  x)n] 
r g a a 


(69 


where 


t - t' 


(70 


and 


Aej  = e(t')  - eQ  = gig2D1K»,no0  (71 

is  the  transient  component  of  strain  existing  immediately  before 
the  stress  is  removed. 

Findley  and  co-workers  [76]  have  found  that  several  of  the 
material  properties  obey  a hyperbolic  sine  function  as  suggested 
by  Eyring's  rate  theory, 


sinh  a/a 


where  the  constants  ag  and  have  values  which  depend  on  the  part- 
icular material.  The  function  gQ  Is  generally  close  to  g * 1 and 
consequently,  c?e  » om  for  many  materials.  Schapery  [35]  has 
shown  that  ao  contributes  most  of  the  stress  dependence  in  Equation 
(72b)  with  gx  e g2  B 1 for  a glass  fiber-reinforced  phenolic 
composite. 

Lamina  Constitutive  Theory 

In  order  to  study  the  behavior  of  laminated  fiber-reinforced 
composites.  It  Is  first  necessary  to  establish  the  constitutive 
theory  of  the  lamina  (or  laminae)  using  linear,  elastic,  aniso- 
tropic theory  [24,  25].  We  shall  first  review  the  basis  for  elas- 
tic theory  before  proceeding  on  to  viscoelastic  behavior  and  sub- 
sequent discussions  on  structure-property  relationships.  In  des- 
cribing the  behavior  of  fibrous  composites.  It  Is  often  helpful  to 
keep  In  mind  the  general  nature  of  the  material  at  hand.  Figure  2 
shows  the  surface  of  a typical  glass  fiber/epoxy  lamina  which  has 
had  the  epoxy  matrix  removed  by  a resin-burnoff  process.  Evidence 
of  some  degree  of  fiber  misalignment  is  clearly  visible.  The  fiber 
ends  are  shown  In  Figure  3,  and  It  can  be  seen  that  there  are  act- 
ually groups  or  bundles  of  fibers.  This  Is  common  in  glass/epoxy 
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composites  because  of  the  fabrication  technique,  but  generally  does 
not  exist  with  boron/epoxy  and  graphite/epoxy  composites. 

Elastic  Anisotropic  Materials 

Governing  equations  of  anisotropic  elasticity.  SI  nee  the 
principles  of  anisotropic  elasticity  represent  the  foundation  for 
the  mathematical  description  of  the  elastic  response  of  composite 
materials,  the  governing  equations  for  this  theory  will  be  presented 
here. 

Using  the  notation  of  Timoshenko  and  Goodier  [136]  we  can  de- 
note the  stresses  acting  on  a small  cubic  element  of  the  material 

In  the  following  manner;  the  three  normal  stresses  are  a , a and 

x y 

a and  the  shearing  stresses  are  t , t and  t where  we  have 
2 xy  Xt  yL 

assumed  symmetry  with  respect  to  = t ^ . The  basic  equations  of 
static  equilibrium  for  an  anisotropic  body  are  the  same  as  those 
given  for  an  Isotropic  body, 

3cx  Stxv  3txz 

+ "1.  + + Y = n 

3X  3y  3z  * u 


9t  3a  3t 

XV  V V7 

—JSZ.  + + —is.  + Y = 0 

3x  3y  3Z 


9txz  3tvz  3oz 

— ^ + —is.  + + i = o 

3 x 3y  3 z 


where  X,  Y and  1 designate  body  forces  referred  to  a unit  volume 
in  the  directions  x,  y and  z respectively. 
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The  infinitesimal  linear  strain-displacement  relations  are 
defined  in  the  same  way  as  in  isotropic  elasticity,  namely. 


30  . „ jw  . 3w 

Ex  “ ax  ' cy  “ ay  ’ Gz  “ az 


(74) 


= M + M + H.  . aw  au 

Yxy  “ay  ax  ’ Yyz  “ay  az  * Yxz  ax  az  * 


In  the  conventional  engineering  definition  of  strain,  the  components 
are  the  same  as  the  general  tensor  notation,  e.g.  Sokol nikoff  [22], 
except  for  the  shear  strain  components  which  are  multiplied  by  a 
factor  of  2.  Equations  (73)  and  (74)  are  tensors  of  order  two  and, 
along  with  the  boundary  conditions  and  constitutive  equations,  spec- 
ify the  state  of  stress  and  strain  at  any  point  within  the  body. 

The  transformation  relations  are  found  in  the  reference  by  Ashton 
et  al . [2]  and  will  not  be  repeated  here  for  the  sake  of  space. 

Generalized  Hooke's  law.  Equations  (73)  and  (74)  are  insuffic- 
ient to  solve  problems  of  equilibrium,  motion  or  stability  of  an 
elastic  body  without  defining  the  relations  between  the  components 
of  stress,  o..,  and  the  strain,  e • • . As  in  Equations  (8)  and  (9), 
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we  can  formulate  the  relations  for  an  elastic  body  to  obtain 
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where  the  are  the  elastic  moduli  and  the  are  the 
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elastic  compliances.  Each  set  of  material  constants  defines  a 
fourth-order  tensor  and  possesses  the  symmetry  discussed  earlier, 
namely,  Equations  (10)  and  (11).  As  a result,  for  a general  aniso- 
tropic material  there  are  21  Independent  elastic  constants. 

A contracted  notation  has  frequently  appeared  In  the  literature. 
In  dealing  with  fourth-order  tensors,  this  notation  reduces  the  num- 
ber of  indices  from  four  to  two  but  expands  the  range  from  three  to 
six.  Iii  contracted  notation,  engineering  strain  is  used  instead  of 
tensorlal  strain.  Table  1 shows  the  correlation  between  normal  and 
contracted  notation  where  e. . represents  tensoriai  strain  given  by 

* J 

+ (76) 

and  the  comma  denotes  differentiation.  A mixed  notation  used  by 
Ashton  et  al . [2]  is  also  shown. 

Table  1.  Conversion  Between  Tensor, Contracted 
and  Mixed  Contracted  Notations 

Tensor  Contracted  Mixed  [2] 


Oil 

ell 

oj 

El 

01 

El 

022 

e22 

02 

E2 

02 

E2 

0 33 

e3  3 

°3 

e3 

03 

E3 

°23 

2e2  3 

04 

E4 

*23 

Y23 

031 

2e3i 

05 

E5 

T 3 1 

Y31 

0 1 2 

2ei2 

06 

E6 

*12 

Y12 

After  changing  to  the  mixed  contracted  notation,  the  constitu- 
tive relation  Is  given 
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where  use  has  been  made  of  the  symmetry  conditions.  This  represents 
the  most  general  case  of  an  anisotropic  material. 

Various  kinds  of  geometric  symmetry  are  often  present,  which, 
in  turn,  lead  to  elastic  symmetry.  Elastic  symmetry  is  expressed 
by  the  property  that  under  certain  coordinate  rotations  the  elastic 
moduli  or  compliances  remained  unchanged.  Some  of  the  constants 
become  equal  to  zero  and  dependencies  appear  between  other  constants 

A monoclinic  material  has  a single  plane  of  symmetry,  which  by 
selecting  a suitable  reference  system  yields 
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with  13  Independent  elastic  constants. 

If  an  anisotropic  material  possesses  two  orthogonal  planes  of 
symmetry,  the  material  Is  called  orthotropic  and  has  the  following 
constitutive  relation 


with  9 independent  elastic  constants.  This  form  occurs  in  many 
structural  materials  such  as  wood,  plywood,  fiber-reinforced  rubber 
and  plastics,  etc. 

When  a material  has  a plane  In  which  the  coefficient  matrix  is 
isotropic,  it  is  called  a transversely  isotropic  material.  If  we 
assume  that  the  y - z (or  2-3)  plane  is  isotropic,  then  there  is 
no  preferred  orientation  In  this  plane  and  the  constitutive  relation 
becomes 
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The  number  of  independent  plastic  constants  for  this  material  is  5 
if  we  assume  that  Si2  = S2i  or  6 if  complete  symmetry  is  not  assum-  ■ 
ed,  e.g., Schapery  [18].  The  Inverse  constitutive  relation  for  the 
transversely  isotropic  material  is  given  by  [5] 
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In  the  case  of  complete  material  symmetry,  the  material  is 
isotropic  and  we  obtain 


ei 

Sll  Sl? 

S 1? 

0 

0 

0 

1 

£? 

Si, 

si? 

0 

0 

0 

'•>? 

c 3 

Si  1 

0 

0 

0 

u ■> 
' 

Y?3 

\ 

2(Sn-S,2) 

0 

0 

T ? 3 

r 3 1 

Symmetri c 

2 ( Si  1 -S 1 2 ) 

0 

TJ1 

Yl? 

2(Sj  1 -S]  2 

*12 

82) 


wKh  2 independent  elastic  constants. 

Engineering  constants  usually  refer  to  Young's  moduli, 
Poisson's  ratio  and  shear  moduli  which  can  be  measured  from  simple 
tests.  The  following  relations  between  the  components  of  S- . and 
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( i , j = 1,  2,  6)  and  the  engineering  constants  can  be  estab- 
lished immediately  from  the  nature  of  uniaxial  and  simple  shear 
tests  for  the  case  of  a transversely  Isotropic  material: 
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where 


V = [(1  + V23) (1  - V23  - 2vi2V2l)] 
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Two-dimensional  composites.  Fibrous  composites  are  generally 
used  in  a manner  such  that  the  stress  state  is  essentially  two- 
dimensional.  As  a result,  the  assumption:,  of  plane  strain  or  plane 
stress  are  often  Invoked.  We  shall  review  the  constitutive  rela- 
tions given  In  Equations  (80)  and  (81)  for  the  transversely  iso- 
tropic material  In  terms  of  plane  stress  assumptions.  This  state 
exists  more  often  than  most  since  the  composite  generally  is  a 
thin  laminated  structure  used  in  plates  or  shells. 

The  plane  stress  problem  can  be  formulated  by  assuming 

G3=T23=t31=0  (84) 

and  we  note  that 

Y23  = Y31  = 0 (8$) 


as  a result  of  Equation  (84). 


Equation  (80)  now  becomes 
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Thus,  for  plane  stress,  the  remains  the  same  as  a three-dimen- 
sional material,  whereas  the  must  be  replaced  by  a "reduced 
stiffness"  denoted  as  the  Q^.  The  Q^.  (1,  j * 1,  2,  6)  are  given 
by 


where 
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This  is  the  constitutive  relationship  for  a specially  orthotropic 
material  in  a plane  stress  state.  It  is  called  specially  ortho- 
tropic when  the  lamina  principal  axes  (1,  2)  coincide  with  the 
reference  axes.  When  the  material  axes  are  referred  to  any  other 
direction  (x,  y)  as  in  Figure  4,  the  constitutive  relations  must  be 
transformed  accordingly,  e.g.  Ashton  et  al.  [2], 
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where  the  and  qu  matrices  are  now  fully  populated  and  the  S.!j 
and  Qy  are  the  transformed  compliances  and  stiffnesses,  respec- 
tively. S,! j and  Qjj  are  related  to  the  orthotropic  lamina  proper- 
ties and  «u  and  the  angle  of  orientation  of  the  lamina,  9. 

The  transformation  equations  are  found  in  the  reference  by  Ashton 
et  al . [2]  and  are  repeated  here  for  convenience, 


S(i  = Sncos4o  + (2S12+S66)sin2ecos2e  + S22s1n4e 

S22  " Sns1n4e  + (2Sx2+S66)s1n2ecos2e  + S22cos4e 

Sl2  - S12(sin48+cos4e)  + (Sn+S22-S66)sin29cos2e  (89a) 

Se6  = 2(2S11+2S22-4Sl2-S65)sin2ecos20  + S66(sin4e+cos4e) 

S 1 q - 2 ( 2S 1 1 -2S i2“Sgg )s i necos 39  - 2(2S22-2Si2-Sgg)sin30cos0 
S2g  - 2 ( 2S j 2 -2S 1 2-Sgg )s 1 n2 9cos0  - 2(2S22-2Si2-Sgg)sm0cos30 

Q11  = Qncos4e  + 2(Q12+2Q66)sin2ecos2e  + Q22sin4e 

Q22  = Qnsin40  + 2(Q22+2Q6s)s1n20cos20  + Q22cos4e 

Qi2  = (Qn+Q22-4Q66)sin29cos2e  + Q12(sin4e+cos4e)  (89b) 

Qee  = (Qn+Q22“2Qi2-2Q66)sin2ecos2e  + Q66(s1n4e+cos4e) 

Qie  “ (Qi 1 -Qi2-2Q66 ) si necos 3e  + (Q12-Q22+2Q66)s1n39cose 
Q26  a (Qn-Qi2-2Qg6)sin39cos0  + {q12-Q22+2Q66)sinecos30 
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Viscoelastic  Composite  Material  Behavior 

The  fibers  and  matrix  materials  used  to  form  fiber-reinforced 
composites  are  usually  assumed  by  most  investigators  to  exhibit 
linear  elastic  behavior  until  they  fail.  In  many  cases,  the  matrix 
material  is  actually  a viscoelastic  material  which  will  exhibit  a 
significant  amount  of  time  dependence  within  certain  ranges  of 
stress  and  temperature,  However,  it  has  only  been  within  the  pest 
decade  that  serious  attention  has  been  given  to  the  viscoelastic 
behavior  of  fibrous  composites. 

The  works  of  Hashin  [34,  137],  Schapery  [9,  18,  32]  and 
Hal  pin  [11,  33]  have  done  much  to  bring  out  the  significance  of 
anisotropic  viscoelasticity.  In  several  recent  literature  surveys 
on  fiber-reinforced  plastic  composites  by  Beckwith  et  al . [6-8], 
it  is  evident  that  these  effects  cannot  be  neglected  in  many  design 
situations.  Experimental  studies  on  the  mechanical  behavior  of 
fibrous  composites  have  been  reported  by  Hanson  [138],  Bott  and 
Barker  [139],  Del  monte  [140],  Zvonar  and  Tamchyna  [141],  Martirosian 
[39-42]  and  Ogibalov  and  Tiuneeva  [36,  37].  Some  of  the  work  done 
in  the  Soviet  Union,  namely  that  of  Antans  and  Skudra  [142,  143], 
Bulavs  and  Skudra  [144]  and  Smotrin  and  Chebanov  [145],  have 
attempted  to  model  the  viscoelastic  behavior  of  glass/epoxy  compos- 
ites using  simple  spring  and  dashpot  models  consisting  of  only  a 
few  elements  typical  of  a Maxwell  or  Voigt  model  [147],  Perhaps 
the  best  documented  work  on  reinforced  polymers  done  in  the  Soviet 
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Union  Is  contained  In  Rabinovich  Cl 46]  and  Tarnopolskiy  and  Skudra 
[147].  Both  references  contain  a vast  amount  of  experimental  and 
theoretical  research  on  glass  and  cotton  fiber-reinforced  composites. 
Time  and  temperature  dependence  are  treated  by  using  spring  and 
dashpot  models.  In  most  cases  it  Is  difficult  to  assess  the 
accuracy  of  these  models  since  their  data  do  not  cover  a very 
long  time  scale  and  is  almost  always  plotted  against  real  time 
rather  than  log  time.  Bryzgalin  [148,  149]  used  the  modified 
superposition  principle  to  predict  the  creep  behavior  of  glass/- 
epoxy  plates,  however,  his  work  is  not  typical  of  the  main  approach 
used  by  Soviet  researchers,  namely,  mechanical  modeling. 

Kaye  and  Saunders  [150]  investigated  the  creep  behavior  of  a 
glass/epoxy  laminate  in  the  linear  viscoelastic  range  and  over  a 
small  temperature  range.  They  were  not  able  to  ascertain  the 
symmetry  of  the  creep  compliance,  i.e.,  S12  = S2l,  as  a result  of 
experimental  limitations.  Cessna  [151],  who  studied  a glass/- 
polypropylene  composite  and  Findley  and  Worley  [152,  153],  who 
studied  a glass  fabric/resin  composite  both  used  the  activation 
energy  theory  of  Eyring  [154]  and  the  hyperbolic  sine  law  to  des- 
cribe the  creep  behavior,  e.g.  Equation  (74). 

Schapery  [18,  32]  has  recently  reviewed  the  use  of  viscoelastic 
analysis  of  fiber-reinforced  composite  materials  and  comments  on  the 
techniques  required  to  characterize  the  material  in  terms  of  linear 
and  nonlinear  viscoelastic  behavior.  Sims  [29]  and  McQuillen  [155] 
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have  also  recently  applied  viscoelastic  theory  to  composite  mater- 
ials using  several  approximate  methods  developed  earlier  by 
Schapery  [18,  32]. 

Linear  constitutive  theory.  Throughout  this  discussion  we 
shall  limit  ourselves  to  the  viscoelastic  behavior  of  an  orthotropic 
material  in  a state  of  plane  stress,  viz.  Equations  (86)  and  (88). 
However,  it  should  be  kepy  in  mind  that  the  general,  three-dimen- 
sional state  of  stress  (strain)  can  be  described  in  a similar 
manner  for  a linear  viscoelastic  anisotropic  material  by  Equations 
(8)  and  (9), 
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In  an  analogous  fashion  we  can  rewrite  Equation  (86a)  as 
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If  the  stress  Inputs  are  of  the  type  = o^hKt),  l.e.,  creep 
test,  then  we  obtain 

^1  - Si i(t)ai  + S12(t)a2 
E2  s Sj2(t)cr1  + S22(t)a2  (91) 

Y12  = Sgs(t)xi2 

where  all  of  the  stresses  are  time-wise  constant.  If  the  lamina 
coordinate  system  is  referred  to  an  arbitrary  set  of  coordinates 
(x,  y)  which  are  not  aligned  with  the  principal  material  directions 
(1,  2)  then  Equation  (88a)  becomes 


I 

where  the  S.  . are  related  to  the  S . . compliances  in  the  same  manner 

I J * J 

as  given  by  Equation  (89a)  except  that  they  are  time-dependent.  As 
in  Equation  (91),  when  the  stress  inputs  are  of  the  type 
= o^H(t),  then  we  have 
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Gx  = Sn°x  + S^ay  + S^Txy 

<y  - s)>ax  + $r/ay  + s?otXy  (93) 

yxy  = S^ax  + S^'ay  + s66Txy 

The  solution  of  the  field  equations,  e.g.,  Equations  (73)  and 
(74),  along  with  any  boundary  conditions  and  either  Equation  (90) 
or  (92)  for  isothermal  conditions  can  be  accomplished  with  the  use 
of  Laplace  transforms  in  order  to  reduce  the  problem  to  an  equival- 
ent elasticity  problem.  This  analogy  is  called  the  "correspondence 
principle"  and  is  limited  to  problems  in  which: 

1.  The  boundaries  do  not  move  except  by  infinitesimal 
displacements, 

2.  the  stress  boundary  conditions  do  not  change  to  dis- 
placement boundary  conditions  with  time,  or  vice-versa, 
and 

3.  the  differential  equations  relating  stress  and  strain 
have  time  independent  coefficients.  If  the  stress-strain 
relations  are  characterized  by  integral  relations,  then 
the  integral  relations  must  be  of  the  convolution  type. 

This  principle  was  shown  by  Lee  [156]  for  isotropic  media,  and  by 
Biot  [157]  for  anisotropic  materials.  Unfortunately  conventional 
Laplace  transform  inversion  is  often  not  feasible  because  the  elas- 
tic solution  is  known  only  numerically,  or  is  so  analytically 


complex  that  standard  methods  are  not  adequate  [18].  Schapery 
[158]  proposed  two  methods  of  transform  inversion,  the  collocation 
method  and  the  direct  method. 

Another  procedure  whlrh  has  been  used  very  extensively  Is  the 
"quasi-elastic"  method  of  analysis  proposed  by  Schapery  [159]. 

The  method  is  easy  to  use  in  that  transform  Inversion  is  avoided. 

In  the  most  general  form,  the  method  Is  equivalent  to  approximating 
the  constitutive  Fquatians  (8)  and  (9)  by 

(94> 

and 

a Si jkjt^cla^ 

where  we  have  assumed 

This  procedure  has  been  applied  successfully  by  Schapery  [32]  and 
Sims  [29]  to  fibrous  composites. 

In  terms  of  general  stress  and  strain  inputs.  Equations  (90) 
and  (92)  may  be  rewritten  as 


ej(t)  - Si](t)o](t)  + Si2(t)o2(t) 
C2(t)  - Si2(t)o,(t)  + S22(t)o2(t) 

Y l j ( t ) - Sg6(t) r i 2 ( t ) 


(97) 
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ex(t)  * Sn(t)ox(t)  + s;?(t)cy(t)  + Si6(t)Txy(t) 

ey(t)  = s;2(t)cx(t)  + s;2(t)0y(t)  + S26(t)Txy(t)  (98) 

Yxy(t)  = Si6(t)ax(t)  + S26(t)oy(t)  + S 66  (t),xy(t) 

Three-dimensional,  nonisothermal  constitutive  equations  for  an 
anisotropic  TSM  have  already  been  given  as 
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where  t;  and  c, ' are  defined  by  Equation  (20).  The  correspondence 
principle  exists  for  a TSM  if  the  temperature  is  spacewise  constant; 
however,  it  may  be  transient  as  long  as  the  first  condition  is 
satisfied.  When  the  temperature  is  both  transient  and  spacewise 
nonuniform,  then  the  correspondence  principle  does  not  exist.  When 
the  material  is  a TCM  of  the  form  of  Equation  (37),  then  the  corres- 
pondence principle  is  similar  to  that  for  a TSM  where  all  the 
stresses  are  divided  by  a^.  The  same  conditions  required  by  a TSM 
must  be  met  by  the  TCM  [18]. 


If  one  (or  more)  phases  of  the  composite  is  of  the  type  TCM-2, 


then  a correspondence  principle  does  not  exist  under  transient 


59 


temperatures.  The  overall  composite  behavior  will  generally  be 
more  complex  than  TCM-2  under  these  conditions.  Schapery  [18] 
discusses  several  other  aspects  of  nonisothermal  behavior  of  com- 
posite materials  in  terms  of  predicting  effective  properties. 

Nonlinear  constitutive  theory.  We  shall  follow  the  previous 
definition  of  material  nonlinearity,  l.e.,  failure  to  satisfy 
homogeneity  and/or  superposition.  In  some  fibrous  composites  one 
important  source  of  nonlinearity  appears  to  be  the  growth  of 
cracks  as  part  of  the  phenomena  we  shall  call  mi  ^structural 
damage.  These  cracks  are  generally  a result  of  the  fabrication 
process  and  the  mismatch  in  physical  and  mechanical  properties  of 
the  matrix  and  fiber.  The  latter  problem  leads  to  significant 
thermal,  or  residual,  stress  in  the  matrix  material  surrounding 
the  fibers  [160]  and  is  further  aggrevated  by  the  high  volume 
fraction  of  fibers  which  creates  tri axial  stresses  between  the 
fibers.  This  damage  and  consequent  nonlinearity  may  be  quite  sig- 
nificant at  relatively  small  strains  (stresses)  compared  to  ultimate 
values.  In  particulate  composites,  e.g.,  solid  propellant,  an 
advanced  stage  of  crack  growth  causes  a large  amount  of  measurable 
dilatation  [161,  162].  Studies  on  several  glass/epoxy  and  graphite/- 
epoxy  composites  by  Beckwith  [163]  at  room  temperature  on  both 
laminae  and  laminates  show  no  detectable  (greater  than  0.1%  volume 
change)  dilatation  prior  to  ultimate  failure. 

Lou  and  Schapery  [9]  and  Schapery  [35]  have  reported  reversible 
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nonlinearities  In  glass/epoxy  composites  due  primarily  to  highly 
stress-dependent  viscosity  behavior.  Their  experiments  were  con- 
ducted after  mechanically  conditioning  the  specimens  until  crack 
growth  effects  reached  a steady  state,  i.e.,  damage  reached  .1  con- 
stant. No  attempt  has  been  made  to  date  to  explicitly  predict  time- 
dependent  micro structural  damage  using  experimental  data  from  fib- 
rous composite  materials. 

Figure  5a  shows  the  forward  dome  of  a typical  solid  rocket 
motor  case  made  from  a glass/epoxy  composite.  A closeup  view 
(Figure  5b)  shows  the  winding  pattern  and  the  light  areas  represent 
regions  of  visible  damage  after  the  pressure  vessel  has  been  loaded 
just  short  of  failure.  The  light  areas  represent  severe  damage  as 
evidenced  by  very  fine  clacks  in  the  epoxy  resin,  limited  fiber 
fracture,  and  delamination.  The  fiber  fracture,  which  is  presumed 
to  occur  only  in  areas  where  the  local  stress  is  close  to  the  ult- 
imate fiber  strength,  is  shown  in  Figure  6 for  a graphite/epoxy 
composite  [12].  Matrix  failure  and  interfacial  adhesive  failure  is 
clearly  visible  in  another  region  of  the  same  composite  as  shown 
in  Figure  7.  Microstructrual  damage  of  this  type  as  not  always  as 
severe  as  shown  here;  however,  the  same  mechanisms  are  thought  to 
exist  at  very  low  stress  (strain)  levels  within  most  structural 
composites  in  use  today. 

Perhaps  at  this  point  we  should  recognize  that  there  are 
additional  classifications  of  viscoelastic  materials  and  divide 
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them  Into  categories  of  fading  memory  and  nonfading  or  permanent 
memory. 

1.  Fading  memory  viscoelastic  theories  are  valid  for 
viscoelastic  materials  whose  time  effects  are 
reversible,  such  as  those  caused  by  Internal  vis- 
cosity effects. 

2.  Permanent  memory  viscoelastic  theories  are  required 
for  viscoelastic  materials  possessing  additional 
time  effects  or  memory  phenomena  not  caused  by  In- 
ternal viscosity  and  are  more  permanent  in  nature, 
e.g.,  micros true tural  damage. 

A single- integral , fading  memory  representation  was  used  by 
Lou  and  Schapery  [9]  to  characterize  a glass/epoxy  composite  In 
the  absence  of  significant  crack  growth  and  similarly  by  Schapery 
[35]  on  a glass/phenolic  fabric  composite.  Equation  (60),  along 
with  Equations  (61),  (69)  and  (74b),  were  used  to  describe  the 
material  behavior.  It  was  found  that  the  nonlinear,  uniaxial  creep 
compliances  obeyed  a power  law  in  time  with  the  exponent  being 
independent  of  fiber  orientation  and  stress  level.  The  nonlinear 
material  properties,  e.g.,  gQ,  g1 , g2  and  a^,  depended  mainly  on 
the  average  octahedral  shear  stress  in  the  matrix.  They  were  able 
to  determine  the  principal  creep  compliances,  e.g.,  S (t),  S22(t) 
and  S12(t),  directly  from  experiments  and  calculated  S66 ( t)  using 
fourth-order  tensor  transformations  and  additional  creep 
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compliance  data.  In  [164],  some  preliminary  work  was  done  to 
relate  the  viscoelastic  behavior  at  various  temperatures  using 
the  Arrhenius  form  of  ay  given  by  Equation  (24).  Additional  work 
done  on  a graphite/epoxy  in  [12]  gave  similar  results,  indicating 
good  agreement  when  T < T . The  main  conclusion  appears  to  be  that 
the  single-integral,  fading  memory  representation  used  in  [9] 
describes  the  material  behavior  only  in  the  absence  of  time-depend- 
ent microstructural  damage.  A constitutive  theory  which  contains 
permanent  memory  effects  is  needed  to  describe  the  behavior  cf 
fibrous  composites  which  generally  exhibit  time-dependent  crack 
growth.  There  are  currently  two  approaches  to  modeling  this  type 
of  behavior;  the  theories  put  forth  by  Farris  [44,  58,  128,  129, 

161,  162,  165,  166]  and  those  by  Schapery  [12,  19,  20]. 

Fitzgerald  and  Farris  [44]  proposed  the  theory  originally 
for  filled  elastomers  such  as  solid  propellant.  Farris  successfully 
applied  the  theory  under  isothermal  [44,  165]  and,  subsequently, 
noni sothermal  [58]  conditions.  Farris  based  his  theory  on  a linear 
cumulative  damage  model  (viz.  Miner's  law)  for  the  failing  micro- 


structure using  the  form 


D'(t)  <tk/tfk) 


where  D‘(t)  is  the  accumulating  damage,  is  the  time  the  sample 
spent  at  strain  level  and  t^  is  the  time-to-fai lure  for  the 
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constant  history  of  strain  level  ek.  An  important  parameter  in 
this  theory  which  results  from  Equation  (99)  is  the  so-called  pth 
order  Lebesgue  norm  of  strain, 


[/i.(f>  iV],/p 


(100) 


where  |e|  is  the  absolute  value  of  strain.  By  judicious  selection 
of  the  Lebesgue  norm  functions,  one  is  able  to  satisfy  the  propor- 
tionality (homogeneity)  condition,  but  not  superposition.  This 
behavior  is  similar  to  that  seen  in  fiber- reinforced  composites 
except  that  stress,  rather  than  strain,  is  the  predominant  factor. 
Motivated  by  these  results  and  his  own  earlier  work,  Schapery  et  al . 
[12,  20]  used  fracture  mechanics,  which  resulted  In  Lebesgue  norms 
of  stress,  to  model  the  rr.icrostructural  damage  in  composite 
materials. 

The  basic  mathematical  concepts  which  describe  the  crack 
growth  behavior  in  viscoelastic  media  are  presented  in  [19].  Using 
these  concepts,  Schapery  [20]  developed  the  constitutive  equations 
for  particulate  composites,  viz.  solid  propellant,  with  micro- 
structural  damage.  The  theory  was  developed  to  predict  micro- 
structural  damage  effects  during  the  initial  stages  of  crack  growth, 
prior  to  extensive  dilatation  within  the  propellant.  It  was 
assumed  that  the  cracks  propagated  in  an  "opening  mode"  and  that  an 
initial  distribution  of  cracks  existed  within  the  matrix,  or 
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possibly  between  the  matrix  and  filler  particles.  It  was  further 
assumed  that  during  most  of  the  time  required  for  local  fracture, 
the  creep  compliance  is  given  by  the  power  law, 

vW  !'0,) 


where  D]  and  n are  material  constants.  Fracture  mechanics  pre- 
dicted Lebesgue  norms  of  stress,  ||o||  , rather  than  strain,  appear 
in  the  general  constitutive  theory. 

If  linearity  for  a fixed  crack  length  is  assumed,  then  the 
total  strain  is  the  sum  of  the  strain  due  to  the  externally  applied 
stress  acting  on  the  undamaged  composite  material  plus  the  strain 
due  to  microcracking.  The  resultant  constitutive  relation  is 
given  by 


(102) 


where  (.  and  c'  take  on  the  normal  definitions  for  reduced-time,  viz. 
Equation  (20),  G(g)  is  a distribution  function  which  reflects  the 
distribution  of  the  crack  lengths  and  stress  concentrations  in  the 
matrix,  and 


g‘  gU')  = 
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where 
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The  exponent  q is  given 


g(C)  « — 

IMIlfc, 


q h 2(1  + 1/n) 


and  ||o||n)q  is  a "weighted"  Lebesgue  norm  where  M = M(c)  accounts 
for  aging  and  rehealing  effects.  In  the  absence  of  these  effects 
M is  unity. 

Equation  (102)  may  be  rewritten  as 


Bmu  - «'>arLdr 


where  oef  is  the  "effective  stress"  seen  by  the  material  and  Is 
given  by 


[oo 

1+/ 


G(g)dg 


°ef  = a[1  + fMal  l]  • 


The  effective  stress  exhibits  both  time  and  temperature  dependence 
in  the  usual  manner  but  also  is  stress-dependent  through  the 
Lesbegue  norm  function  represented  by  f||o||.  It  may  also  be  con- 
venient to  work  with  the  inverse  relationship  given  by 
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°ef  ^ Em^  ~ df'  (108) 

The  constitutive  theory  represented  by  Equations  (106)  and  (108) 
was  developed  under  the  assumption  that  the  particles  are  rigid 
relative  to  the  matrix.  Tnerefore,  the  crack  growth  is  governed  by 
the  stress  intensity  factor  and  fracture  energy  associated  with 
the  matrix  material.  The  theory  represents  the  condition  of  cracks 
which  are  initially  isolated  from  each  other  and  do  not  interact. 

It  has  been  found  that  the  cracks  actually  interact  in  some 
materials  depending  on  the  constituent  material  properties. 

Studies  using  a scanning  electron  microscope  on  qraphite/eooxy  and 
glass/epoxy  composites  show  that  a crack  arrest  mechanism  is  present 
under  certain  loading  conditions  [12],  Figures  8 and  9 show  the 
crack  arrest  mechanism  and  the  degree  of  fiber  bending  present  in 
an  E-glass/epoxy  laminae  which  has  been  subject  to  a stress  normal 
to  the  fibers. 

Motivated  by  this  possibility,  Schapery  et  al . [12]  studied  the 
effects  of  fiber  bending  on  the  overall  strain  energy  release  rate 
as  a function  of  the  loading  condition  and  composite  properties. 

The  energy  associated  with  matrix  stretching  was  ultimately  neglect- 
ed and  the  effects  of  fiber  bending  and  shear  stress  were  used  to 
model  the  crack  growth.  Equation  (101)  was  again  used  to  model  the 
matrix  creep  compliance  and  the  distribution  function  representing 
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crack  lengths  and  stress  concentrations  was  assumed,  as  an  example, 
to  be 


Fi  = Kir. 


-P 


(108) 


where  Kx  and  p are  positive  constants  with  p>l . The  model 
allows  for  the  possibility  of  crack  arrest,  a phenomenon  which  is 
clearly  evident  in  the  behavior  of  multiple  cycle  creep  and  recov- 
ery tests  on  fibrous  composites  [12,  167]  in  which  most  of  the 
damage  is  done  during  the  first  loading  cycle. 

The  principle  compliance  parallel  to  the  fibers,  Si  i , and  the 
Poisson's  ratio  for  loading  in  the  fiber  direction  will  not  be 
affected  by  cracks  parallel  to  the  fibers.  Therefore,  as  long  as 
appreciable  fiber  failure  does  not  occur,  the  principal  stress- 
strain  equations  are  identical  to  those  for  linear  viscoelasticity 
except  for  additional  terms  of  Ay12  and  ac2.  Hence,  the  principal 
relations  are 
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As  an  example  [12],  we  can  rewrite  the  shear  strain  as 


S*6(t“t)-pp  dT 


where  S66  is  the  principal  nonlinear  compliance,  and  is  given  by 


n 1/  / % 4-2p 

c*  _ c x Vi  4Asin4e  + Bsin22e  n n /t  \2-n 
S“  • S«  * (p-D'(2-p)  j 4T j °o  <5^ 


for  the  case  of  a uniaxial  creep  test  where  a = aQH(t).  In 
Equation  (111),  r is  the  fracture  energy  of  the  matrix,  and  A and  B 
are  constants  which  depend  on  the  constituent  properties  and  frac- 
ture properties.  Their  explicit  definition  is  given  by  Equation 
(39)  of  [12]. 

Neither  theory,  Equations  (102-108)  or  Equations  (109-111), 
has  been  compared  with  a specific  fibrous  composite  material  al- 
though they  appear  to  be  general  enough  to  predict  the  types  of 
nonlinearity  presently  observed. 

Structure-Property  Rel ationships 

Micromechanics  deals  with  the  mechanical  interaction  between 
the  constituent  materials  of  a composite.  The  study  of  micro- 
mechanics  of  composite  materials  has  received  considerable  atten- 
tion since  the  early  1 960 ' s in  an  attempt  to  relate  constituent 
properties  to  the  overall,  macromechanical  response  of  a single 
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ply,  or  1 ami nd.  We  have  already  shown,  vi/.,  liquation  ('{(.),  l.hat 
the  constitutive  rt;  lotion  lor  an  orthotropic  iiinteri.il  in  a state 
oF  plane  stress  depends  on  the  determination  of  four  compliance:;, 
naniely,  S M , 5,.,,  Sp:,,  and  S , , , , . It  has  also  been  rioted  in  l.nuolion 
(8.3)  that  these  material  properties  can  he  related  to  functions  of 
the  ernj  incori  nq  constants  I.  ■ j , l;v,  R,;.  and  ci  tlier  v!:.  or 
The  latter  two  are  related  through  the  symmetry  condition  that 
requ i res 

'■I  I y 1".  ;j  ;>  = Vp  ! h . . ( i I r’  ) 

so  that  there  are  again  only  four  independent  n:a  tori  a I cons  tan  !.■ 

I.o  be  determined.  As  wo  noted  earlier,  the  viscoelastic  iiatmo  of 
the  polymeric  matrix  also  makes  these  properties  time  and  l.emoer- 
al.ure  doperidcnt,  with  G];)  (or  SM()  generally  exhibiting  the  pre- 
dominant viscoelastic  effects. 

Schapery  [18]  discusses  the  measurement  of  the  principal  com- 
pliances experimentally  using  off-angle  unidirectional  composites 
and  fourth-order  tensor  transformations.  The  direct  measurement 
of  S;|,  S;.:.,  S | cjnd  S,.,,  f a 1 1 s i nto  the  classification  of  macro- 
mechanics  where  the  ply  or  lamina  is  treated  as  a single  material 
with  anisotropic  properties,  neglecting  the  overall  constitution  of 
the  lamina.  Several  experimental  techniques  have  been  proposed  lor 
the  measurement  of  these  properties  directly  or  indirectly  [2, 

18,  33,  168], 


In  the  subsections  which  follow  we  shall  discuss  some  of  the 
aspects  of  current  micromechanics  theory,  paying  ^articular  atten- 
tion to  the  relationship  between  the  constituent  properties  and  the 
overall  response  observed  experimentally. 

Basic  assumptions.  There  are  several  basic  assumptions  which 
are  common  to  almost  all  of  the  theories  presently  in  use  today 
with  regard  to  elastic  behavior,  namely: 

1.  The  ply  is  macroscopically  homogeneous,  linearly 
elastic  and  orthotropic. 

?..  The  fibers  are  linearly  elastic  and  homogeneous. 

3.  The  matrix  is  linearly  elastic  and  homogeneous. 

4.  Both  the  fiber  and  matrix  are  free  of  voids. 

5.  There  is  complete  bonding  at  the  interface  of  the 
constituents  and  there  is  no  transitional  region 
between  them. 

6.  The  fibers  are  (a)  regularly  spaced  and  (b)  aligned. 

7.  In-si tu  mechanical  properties  are  the  same  as  the 
properties  of  the  constituents  when  not  in  the 
composite. 

In  general , it  is  realized  that  the  lamina  often  exhibits 
behavior  which  is  nonlinearly  viscoelastic,  particularly  at  high 
stress  levels  and  temperatures.  There  are,  of  course,  a great 
number  of  possible  deviations  from  these  assumptions  and  the  net 
result  is  an  abundance  of  theories  based  on  relaxing  one  or  more 
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of  the  basic  assumptions. 

Prediction  o_f  effective  properties.  We  shall  not  attempt  to 
critique  the  various  micromechanics  theories  since  time  and  space 
do  not  permit  even  a minimal  treatment.  The  interested  reader  is 
referred  to  several  excellent  critiques  of  the  area  given  by 
Chamis  and  Sendeckyj  [169],  Ashton  et  al.  [2],  Pagano  and  Tsai 
[170]  and  others  [171,  172].  The  various  theories  which  have  been 
proposed  to  date  can  be  classified  as  follows:  netting  analysis, 
mechanics  of  materials,  sel f-consistent  models,  variational  methods, 
exact  (elasticity)  methods,  statistical  methods,  discrete  element 
methods  and  semi -empirical  approaches.  An  exhaustive  bibliography 
of  many  of  these  theories  can  be  found  in  [6-8,  169]  and  will, 
therefore,  not  be  repeated  here. 

The  prediction  of  the  longitudinal  composite  properties, 
viz.,  En  and  v12»  has  generally  proceeded  along  the  classical  lines 
of  parallel  element  models  [137,  173-176].  The  longitudinal  com- 
posite modulus,  En,  is  given  by  [137,  174] 


v.E.  + vE  + vv.Gm 
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(113) 


where  E is  the  modulus  of  elasticity,  G is  the  shear  modulus, 

K is  the  plane  strain  bulk  modulus,  v is  the  volume  fraction,  and 
the  subscripts  f and  m signify  the  fiber  and  matrix  respectively. 
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Hill  [174]  demonstrated  that  En  is  actually  bounded  from  below  by 
the  "rule  of  mixtures"  and  that  the  last  term  can  be  neglected  such 
that 


Ei i = VrEr  + v_E  . 
11  f f mm 


(114) 


Tsai  [172]  allowed  for  the  possibility  for  fiber  misalignment  by 
applying  a correction  factor  to  Equation  (114) 


Eli  = k(v,E,  + vmEj 
11  ' f f mm' 


(115) 


however,  the  value  of  k is  often  very  close  to  unity  for  present 
composites  considering  the  fabrication  techniques  used. 

The  major  Poisson's  ratio,  v12,  was  also  derived  in  a form 
similar  to  Equation  (113),  namely,  [137,  174] 
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(116) 


where  v is  Poisson's  ratio  and  the  other  symbols  are  the  same  as 
previously  defined.  Again,  neglecting  the  last  term  we  have 

»i2  “ Vf  + V,n  . (117) 

Most  experimental  evidence  indicates  that  the  "rule  of  mixtures" 
representation,  viz  , Equations  (114)  and  (117),  are  auite  accurate 
for  the  calculation  of  En  and  v12  respecti vely . Extension  to 
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include  viscoelastic  behavior  can  be  readily  accomplished  using 
the  quasi-elastic  approach, 

El  i(t)  = vfEf (t)  + vmEm(t)  (118a) 

and 

vj2(t)  * vfvf(t)  + vniv.n(t)  (118b) 

where,  for  most  advanced  composites,  the  properties  of  the  fiber 
are  taken  to  be  independent  of  time,  except  possibly  at  very  high 
temperatures . 

The  representation  of  the  composite  transverse  modulus,  E22> 
and  shear  modulus,  Gi2,  is  not  quite  so  elementary;  also,  these 
properties  are  much  more  sensitive  to  voids  and  fiber  geometric 
arrangement.  Many  investigators  have  looked  at  square  arrays,  rec- 
tangular arrays,  hexagonal  arrays  and  random  fiber  geometries  in 
various  attempts  to  model  the  actual  behavior.  Hal  pin  [171]  and 
Ashton  et  al.  [2]  developed  some  approximate  forms  based  on 
Herman's  work  [177]  who  originally  used  a method  developed  b,y 
Kerner  [178],  These  equations  are  based  on  the  application  of 
semi -empirical  adjustment  factors  to  the  theoretical  representa- 
tions of  the  orthotropic  engineering  parameters.  These  adjustment 
factors  depend  on  the  fiber  geometry  and  spacing  and  can  possibly 
account  for  void  and  microstructural  damage  effects.  The  "Halpin- 
Tsai  equations"  are  given  by 
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(119b) 


where  .j£  and  are  the  adjustment  factors.  Reliable  estimates  for 
the  c-factor  can  be  obtained  by  Equations  (119a)  and  (119b)  with 
various  numerical  micromechanics  solutions.  Typical  values  for 
elastic  composites  are  = 2 and  cG  = 1 [2].  Deviations  from 
these  values  account  for  differences  in  the  assumed  mi crostructural 
arrangement  of  the  fibers  and  mi crostructural  damage.  As  in 
Equation  (118)  a quasi-elastic  viscoelastic  solution  may  be  ob- 
tained as 


E22(t) 


P i.x  r Ef(t)(i+cEvf> + i 
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(120a) 


Sl2(t)  = Gm(t) 


Gf(tm+5Gvf)  + cGyt)d-Vf)  * 
G fTt)  (l-vf)  + Ta^tJTVg+v'fT" 


(120b) 


where  we  can  normally  assume  that  and  are  independent  of  time. 

The  principal  elastic  compliances  Sn,  S22*  S12,  and  S66  can 
be  found  by  applying  Equation  (83)  to  the  engineering  properties. 
Corresponding  viscoelastic  compliances  are  derived  by  using  Equation 
(83)  and  the  quasi-elastic  method  to  obtain  Sn(t),  S22(t),  S;2(t) 
and  i' (.,-,( t ) . We  can  also  include  temperature  dependence  by  using  the 
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reduced-time,  rlt  provided  we  do  not  violate  the  correspondence 
principle  discussed  earlier.  Halpin  [33]  has  used  the  Halpin-Tsai 
equations  to  predict  and  for  a viscoelastic  rubber  rein- 
forced wi th  unidi rectional , nylon  fibers  over  a temperature  range 
using  isothermal  data.  Sims  [29]  has  also  employed  these  relations 
analytically  to  calculate  the  relaxation  moduli  for  unidirectional 
graphite  and  boron  fiber-reinforced  epoxy  composites. 

Laminate  Constitutive  Theory 

The  behavior  of  laminated  composite  materials  is  directly 
related  to  the  response  of  the  various  laminae  after  taking  into 
account  the  geometry  and  individual  properties  of  the  layers. 
Laminated  composite  theory  is  given  in  Ashton  et  al . [2],  Ashton 
and  Whitney  [27],  Ambartsumyan  [26],  and  Lekhnitskii  [25]  for 
elastic  materials  and  may  be  extended  to  linear  viscoelastic  be- 
havior by  using  the  correspondence  principle  or  the  quasi -el astic 
approach . 

A typical  laminate  consisting  of  n layers,  or  lamina,  is 

shown  in  Figure  10.  The  strain-displacement  relationships  are 

generally  derived  for  small  deformations  of  the  laminate  and  are 

written  in  terms  of  the  midplane  strains  e* , e° , 

curvatures  k , k and  k as  follows, 

* y 


and  the  plate 
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^X  " ~3X 


cy  " 3y 
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v°  = —2.  + 2. 

'xy  3y  3x 


where  uo  and  vo  are  the  midplane  displacements  and  w is  the  trans- 
verse displacement.  The  strains  in  any  given  layer  are  related  to 
the  midplane  strains  and  curvatures  in  Equation  (121)  by  the 
relationships 


. = c°  + zk 

XX  X 


ey  = ey  + zky 


v = y + ZK 

rxy  rxy  xy 


where  z is  the  coordinate  distance  measured  perpendicular  to  the 
laminate  reference  surface  as  shown  in  Figure  10. 


Elastic  Orthotropic  Materials 


Under  the  assumption  that  plane  sections  remain  plane  the 
constitutive  equations  for  an  elastic,  orthotropic  laminate  con- 
sisting of  n laminae  having  the  constitutive  relationship  which 
obeys  Equations  (88b)  and  (89b)  are  given  as  (see,  for  example, 
Ashton  et  al . [2]) , 
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where  N N and  N are  the  in-plane  forces  (stress  resultants) 
x . y xy 

and  M I , M and  M are  the  moments  (moment  resultants)  given  by 
x y xy 
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N y = / °ydz  (124a)  My  = / °yzdz  (1241 

-h/2  -h/2 


N = / t dz 

xy  / xy 

-h/o 


Mxy  = 7 Vdz 

-h/2 


The  quantities  A..,  B-.  and  D.  . are  the  laminate  extensional  stiff 
1 J I J J 

nesses,  coupling  stiffnesses  and  bending  stiffnesses,  respectively 
These  laminate  stiffnesses  can  be  defined  in  terms  of  the  lamina 
stiffnesses  given  in  Equation  (89b)  as 
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n 

Aij  “ Xv  (QijMhk  “ hk-l^ 
k=l 

n 

BiWZ(QijVhk-hLi>  <l2! 

k=l 

n 

Dij  = 3 S (Qi jJk(hk  " hk-l^ 
k=l 

where  h^  is  the  coordinate  distance  to  a lamina  from  the  reference 
surface  as  shown  in  Figure  10.  The  Q^'s  are  the  transformed 
stiffnesses  which  depend  on  the  orthotropic  lamina  stiffness  and 
the  angle  of  orientation  of  the  lamina  as  shown  in  Equation  (89b). 

The  most  important  feature  of  Equation  (123)  is  the  coupling 
phenomena  which  exists  between  stretching  and  bending  through  the 
B . matrix.  If  the  (q!.),  is  an  even  function  of  z (symmetric 

1 J ' J K 

layup  of  the  laminate),  then  B.  ■ = 0 and  coupling  is  eliminated. 

* 

Many  laminates  are  constructed  in  this  manner  and,  therefore,  the 
governing  constitutive  equations  are  considerably  simplified. 

Viscoelastic  Orthotropic  Materials 

The  constitutive  equations  for  a linear  viscoelastic  laminate 
can  be  formulated  by  replacing  the  products  of  the  time-dependent 
quantities  in  Equation  (123)  by  superposition  integrals  of  the 


form 
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and,  similarly, 
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where  the  Q , 1 s are  the  time-dependent  transformed  stiffnesses 

I vJ 

defined  by  Equation  (89b). 

In  the  application  of  Equations  (123),  (126)  and  (127)  for  the 
solution  of  many  viscoelastic  problems,  the  procedure  often  becomes 
time  consuming.  Under  certain  conditions  when  the  time-dependent 
material  functions  or  input  quantities  are  slowly  varying  functions 
of  time,  the  convolution  integral  may  be  omitted.  The  quasi-elastic 
approach  used  earlier  to  derive  the  lamina  constitutive  equations 
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may  also  be  used  here.  The  procedure  involves  direct  substitution 
of  the  time-dependent  properties  of  the  A'.,,  8 ! . and  D*.  . matrices. 
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Structure-Property  Relationshi ps 


In  order  to  apply  the  inicromechanics  theory  derived  earlier 
to  the  behavior  of  the  laminate,  it  is  necessary  to  establish  the 
relation  between  the  S,  . and  Q.  ..  These  relations  are  given  in 
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where  SM,  So;;,  Si?  and  SGf,  may  either  be  measured  directly  or 
predicted  by  the  various  inicromechanics  theories  already  discussed. 
The  Halpin-Tsai  relations,  Equations  (114),  (117)  and  (119),  are 
often  used  because  of  their  relative  simplicity.  We  may  also 
relate  the  Q. . to  the  common  engineering  constants  and  apply  the 

1 o 

Halpin-Tsai  relations  directly,  [2] 
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El 1 
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Q?2  = rv  vf2V'21  (87b) 

Q = V21E’l  1 _ VI2E22 

^12  1 - Ul2V21  1 - V12V21 

Q6  e s ®i 2 

Iri  either  case,  the  application  of  the  quasi-elastic  approach  to 
the  solution  of  viscoelastic  problems  proceeds  by  inputting  time- 
dependent  properties  in  Equations  (128)  and  (87b)  by  using,  for 
example,  time-dependent  Halpin-Tsai  relations.  Equations  (118) 
and  (120). 
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SECTION  III 

MECHANICAL  CHARACTERIZATION  TESTS 

The  experimental  program  was  designed  to  provide  the  necessary 
isothermal  data  required  to  evaluate  the  effective  (overall)  lam- 
inate stiffness  of  a typical  glass/epoxy  composite  material.  The 
glass/epoxy  fibrous  composite  studied  is  currently  being  used  in 
the  construction  of  the  third  stage  Minuteman  III  solid  rocket 
motor  case.  State-of-the-art  filament  winding  techniques  are  used 
to  fabricate  the  actual  case.  The  basic  barrel  (cylindrical  sec- 
tion) consists  of  14  layers  of  hoop-oriented  S-901  glass  wraps  and 
8 layers  of  glass  wraps  oriented  14.63°  to  the  longitudinal  motor 
axis.  The  S-901  glass  filament  rovings  are  used  in  a tape  prepreg 
system  consisting  of  a Shell  58-68R  epoxy  resin. 

Glass  cloth  is  interwoven  between  hoop  layers  at  the  barrel 
ends  to  provide  stub  skirts  for  attachment  purposes.  In  the  for- 
ward and  aft  dome  regions  the  fiber  or  wrap  angle  is  variable  be- 
cause of  the  shape  of  the  end  closures.  The  geometry  Is  further 
complicated  with  the  provision  for  six  thrust  termination  ports  in 
the  forward  dome.  These  areas  are  reinforced  locally  with  a glass 
tape/epoxy  during  the  winding  process. 

As  a result,  the  effects  of  both  variable  thickness  and  wrap 
angle  must  be  taken  into  account  in  the  forward  dome  area  in  order 

’Formerly  designated  as  S-HTS  or  S-994-HTS  glass.  Owens- 
Corning  Fiberglas  Corporation  trademark. 
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to  determine  the  effective  stiffness  as  a function  of  the  fiber 
angle.  In  the  discussion  of  effective  stiffness  in  Section  II  it 
was  noted  that  the  controlling  factors  are  primarily  constituent 
properties,  volumetric  content  and  stacking  sequence.  The  fiber 
content  may  be  somewhat  variable  in  the  dome  regions  due  to  the 
fabrication  technique  which  causes  both  fiber  slippage  ("gaps") 
and  shingling  ("overlapping").  Current  motor  case  analysis  con- 
ducted by  the  Aerojet  Solid  Propulsion  Company  ( ASPC ) 2 assumes  a 
nominal  65.7%  fiber  content  (volume)  as  representative  of  the 
actual  motor  case  [17]. 

The  mechanical  characterization  tests  described  in  this  section 
were  conducted  in  a manner  that  would  permit  the  evaluation  of  the 
effective  stiffness  in  the  glass/epoxy  composite  as  a function  of 
fiber  angle,  stacking  sequence  and  load  history.  The  effects  of 
multiple  loading  cycles,  typical  of  the  proof testing  of  glass/epoxy 
rocket  motor  cases,  and  temperature  were  considered  in  the  latter 
area. 


Materials  and  Equipment 


Materials 

The  prediction  of  effective  properties  through  the  use  of  the 
Halpin-Tsai  Equation  (119)  and  other  associated  micromechanics 
theories  rely  very  heavily  upon  accurate  measurements  of  the 


;'Aerojet  Solid  Propulsion  Company,  Sacramento,  California. 
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constituent  properties.  In  the  glass/epoxy  system*  the  glass  fiber 
properties  are  generally  assumed  to  be  linearly  elastic  up  to  fail- 
ure. Epoxy  matrix  properties,  on  the  other  hand,  exhibit  consider- 
able time  and  temperature  dependence. 

The  materials  tested  were  representative  of  both  the  consti- 
tuent materials  as  well  as  the  actual  rocket  motor  case.  Four 
materials  were  used  in  the  experimental  program: 

1.  Shell  58-68R  epoxy  resin. 

2.  S-901  glass/Shell  58-68R  epoxy  resin  unidirectional 
composite  laminae. 

3.  S-901  glass/Shell  58-68R  epoxy  resin  composite 
laminates. 

4.  S-901  glass/Shell  58-68R  epoxy  resin  laminate 
sections  removed  from  a third  stage  Minuteman  III 
solid  rocket  motor  case. 

The  first  three  materials  were  fabricated  by  Structural  Composites 
Industries  (SCI)3  for  the  experimental  program  and  furnished  under 
a subcontract  with  the  ASPC.  All  of  these  materials  were  made 
according  to  the  current  ASPC  specifications  pertaining  to  the 
materials  preparation  and  method  of  construction  as  actually  used 
in  the  third  stage  Minuteman  III  motor  case. 

Epoxy  resin.  The  Shell  58-68R  epoxy  resin  system  was  orig- 
inally developed  by  SCI  in  the  early  Polaris  solid  rocket  motor 


■^Structural  Composites  Industries,  Azusa,  California. 
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cases  and  has  been  used  by  NASA  In  several  filament-wound  tank 
applications  for  cryogenic  fuels  [179,180],  The  epoxy  system  used  In 
our  tests  employs  the  standard  formulation  given  In  Table  2. 


Table  2.  Ingredients  of  Shell  58-68R  Epoxy  Resin  [181] 


Material 

Formulation: 
Parts  by  Weight 

Epon  828 

50 

Epon  1031 

50 

Nadic  Methyl  Anhydride  (NMA) 

90  ± 5 

Benzyl dimethyl  amine  (BDMA) 

0.55  ± 0.05 

The  epoxy  mixture  was  prepared  by  mixing  the  Epon  828  and 
Epon  1031  components  In  an  oven  at  150°F  and  then  blending  them 
together  once  the  Epon  1031  was  completely  melted.  While  the 
components  were  still  warm,  the  NMA  (curing  agent)  and  BDMA 
(accelerator)  were  added  separately  and  mixed  thoroughly  after  each 
addition.  The  resulting  mixture  was  placed  In  a vacuum  bell  jar 
in  order  to  remove  any  entrapped  air  bubbles. 

After  the  Initial  mixing  and  degassing,  the  epoxy  mixture  was 
poured  Into  flat  casting  molds  with  8”  x 8"  x 1/8"  dimensions. 
Sufficient  quantity  of  the  epoxy  was  prepared  to  cast  three  sheets 
from  the  same  batch.  The  molds  were  then  placed  in  an  oven  in  a 
vertical  position  and  cured  for  2 hours  at  200°F  and  then  2 hours  at 
350°F.  Upon  removal  from  the  molds,  the  plates  were  visually  estab- 
lished to  be  void  free  and  found  to  have  a uniform  thickness  of 


0.125  ± 0.003  Inches.  The  density  of  the  cured  plates,  measured 
by  water  displacement  method,  was  found  by  SCI  to  be  1.24  grams/cc 
[181],  This  Is  consistent  with  other  studies  using  the  same  resin 
system  with  slight  variations  of  cure  times  and  temperatures 
[179,  180]. 

Glass/epoxy  composites.  Both  the  unidirectional  laminae  and 
the  composite  laminates  were  prepared  by  using  the  Shell  58-68R 
epoxy  resin  and  Owens-Corning  S -901  Glass  roving  (20-end)4.  The 
20-end  roving  is  used  in  the  actual  motor  case  because  experience 
has  shown  it  to  have  a higher  strength  than  the  12-end  roving 
which  is  also  available  [179]. 

All  prepreg  broadgoods  materials  required  for  panel  fabrica- 
tion were  prepared  by  a drum-winding  process.  This  process  involved 
winding  the  S -901  Glass  roving  onto  a Myl ar-film-1 ined,  24-inch 
diameter  cylinder  at  a controlled  winding  speed  and  load,  and  then 
applying  a quantity  of  the  Shell  58-58R  epoxy  resin  uniformly  over 
the  entire  surface.  The  amount  of  glass  roving  and  epoxy  required 
for  the  prepregs  was  precalculated  to  yield  a desired  ply  thickness 
and  glass  volume  percent  close  to  65%  in  the  molded  panels-  To 
facilitate  processing,  a 50%  resin  solution  was  made  in  methyl - 
ethyl-ketone  (MEK).  Residual  solvents  in  the  prepregs  were  removed 
by  heating  in  an  oven  at  150°F  prior  to  the  final  layup. 

4An  "end"  is  defined  as  the  smallest  commercially  available 
bundle  of  glass  filaments;  each  "end"  generally  consists  of  204 
monofilaments  in  a single,  continuous,  untwisted  strand. 


The  prepreg  broadgoods  were  cut  into  13-inch  strips  with  the 
fibers  running  at  the  required  angle  of  the  molded  panel.  Single 
ply  patterns  approximately  13"  x 13"  were  cut  and  fifteen  (15) 
plies  were  laid  up  In  a single  stack.  All  panels,  unidirectional 
and  laminated,  had  fifteen  plies  such  that  the  outer  layers  were 
oriented  in  the  same  direction  in  order  to  reduce  warpage  during 
the  curing  process.  The  plies  were  also  stacked  symmetrically 
such  that  the  bending-stretching  (coupling)  stiffness,  B.  .,  is  zero. 

* J 

The  prepreg  plies  were  laid  up  on  polished  aluminum  plates,  covered 
with  a Teflon-impregnated  release  fabric,  followed  by  another 
13"  x 13"  x 0.065"  aluminum  cover  plate  over  the  outer  layers,  and 
then  subsequently  covered  with  2 plies  of  glass  bleeder  cloth. 

The  entire  assembly  was  bagged  with  a Nylon  film  and  placed  in  an 
autoclave  for  cure.  Curing  was  executed  for  2 hours  at  200"F  and 
then  2 hours  at  350°F  in  a 75  psig  autoclave  pressure  in  addition 
to  the  vacuum-bag  pressure. 

The  cured  panels  were  inspected  for  visible  defects.  The 
panel  thickness  was  measured  at  25  locations  and  then  trimmed  to  a 
net  delivery  size  of  12"  x 12"  x 1/8".  Unidirectional  glass/epoxy 
panels  were  fabricated  at  fiber  angles  of  0,  20,  45  and  90°  while 
the  laminated  glass/epoxy  panels  were  made  in  symmetric  layups  of 
0/90,  ±10,  ±30,  ±45,  ±60  and  ±80°. 

Filament-wound  case  materials.  The  third  stage  Minuteman  III 
solid  rocket  motor  case  is  shown  in  Figure  11.  The  motor  shown  is 


instrumented  and  being  prooftested  by  ASPC  prior  to  casting  the 
solid  propellant.  The  barrel  section  along  with  the  strain  gage 
instrumentation  can  be  clearly  seen.  Figure  12  shows  the  section 
of  the  forward  dome  which  has  been  removed  from  the  barrel  section. 
The  thrust  termination  ports  and  igniter  adapter  are  shown  in  more 
detail.  In  Figure  11  they  are  hidden  by  the  interstage  assembly 
in  the  forward  end  of  the  motr  asa.  A closeup  of  the  area  be- 
tween two  of  the  thrust  termination  ports  is  shown  in  Figure  13. 

In  order  to  make  a comparison  between  the  effective  stiffness 
of  the  composite  plates  and  the  actual  motor  case,  the  ASPC  furn- 
ished both  the  forward  dome  (Figure  12)  and  the  aft  end  of  a third 
stage  motor  case  which  had  undergone  prooftesting  ("hydrotesting"). 
Although  the  effects  of  the  hydrotesting  damage  may  be  significant 
with  regard  to  the  first  loading  cycle,  the  results  of  subsequent 

A,  L 

loading  cycles  (second  through  nin  cycle)  could  be  compared  directly 
with  the  glass/epoxy  plates  fabricated  by  SCI.  For  these  series 
of  tests  the  forward  dome  was  chosen  since  the  aft  end  has  an 
outer  layer  of  cork  insulation  which  could  not  easily  be  removed 
without  damaging  the  case  material. 

Samples  were  obtained  from  an  area  between  two  of  the  thrust 
termination  ports  similar  to  that  shown  in  Figure  13.  In  this 
region  the  motor  case  does  not  possess  mid-plane  symmetry  since 
there  are  eight  (8)  alternating  layers  and  hence  f 0.  The 
fiber  angles  of  samples  taken  from  the  meridional  and 
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circumferential  directions  were  approximately  ±20  and  ±70° 
respectively. 

Aluminum.  A 6061 -T6  aluminum  plate  was  also  used  as  a refer- 
ence material  for  comparison  with  the  glass/epoxy  plate  tests 
which  are  discussed  in  a subsequent  section. 

Specimen  Preparation 

Uniaxial  tensile  coupons.  The  majority  of  the  tests  conducted 
during  the  course  of  the  experimental  program  utilized  uniaxial 
tensile  coupons  cut  from  the  plates  fabricated  by  SCI.  Sample 
preparation  technique  for  both  the  Shell  58-68R  epoxy  resin  and 
S-901  glass/Shell  58-68R  epoxy  resin  materials  was  the  same. 

The  surfaces  of  the  plates  which  were  to  be  used  were  first 
taped  completely  with  masking  tape.  A sample  layout  corresponding 
to  the  desired  fiber  angle  and  coupon  dimensions  was  drawn  directly 
on  one  taped  surface.  All  sample  dimensions  were  drawn  slightly 
oversize  to  allow  for  final  trimming.  Using  a tungsten  carbide 
band  saw  blade  and  a slow  band  speed,  the  samples  were  rough  cut 
from  each  plate  used  in  the  tests.  After  carefully  removing  the 
masking  tape,  several  samples  were  mounted  in  a vertical  mill  and 
machined  to  the  final  dimensions  using  a fly  cutting  attachment. 
Visual  Inspection  of  the  sample  edges  showed  little  damage  from 
the  cutting  operation  provided  the  cutting  tool  was  resharpened 
after  each  surface  was  completed. 
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Several  sets  of  samples  were  prepared  from  the  SCI  plates  at 
different  fiber  angles,  o,  with  the  emphasis  placed  on  the  off-angle 
specimens  (0°  •-  o • 90°).  These  latter  specimens  exhibit  significant 
creep  and  nonlinearity.  All  of  the  samples  were  cut  to  the  nominal 
dimensions  of  6"  x 1/2"  x 1/8". 

Specimens  for  the  failure  tests,  constant  crosshead  rate  tests 
and  the  creep  and  recovery  tests  required  end  tabs  to  be  bonded  to 
coupons  in  order  to  transfer  the  applied  load  to  the  specimen. 
Aluminum  end  tabs  were  sandblasted  and  bonded  to  the  coupon  ends 
using  Micro-Measurements'  M-Bond  AE-155  adhesive  cured  under  slight 
pressure  for  2 hours  at  170°F.  The  end  tabs  for  the  constant  cross- 
head rate  tests  were  basically  5°  wedges  measuring  1.5"  x 1"  x 1/0" 
and  the  tabs  for  the  failure  and  creep  and  recovery  tests  measured 
1 11  x 5/8"  x 1/8"  with  the  edge  tapered  to  about  20-30°  (Figure  14). 
Different  types  of  tabs  were  used  because  of  the  difference  in 
grips. 

Tensile  coupons  used  for  the  creep  and  recovery  tests  were 
loaded  by  means  of  pin-connected  grips  while  the  constant  crosshead 
rate  coupons  were  held  by  a set  of  wedge-action  grips.  In  order 
to  reduce  the  possibility  of  load- induced  bending  due  to  any  pin 
misalignment,  an  alignment  fixture  was  used  to  clamp  and  align  the 
ends  while  a 1/8-  or  3/16-inch  hole  was  drilled  through  the  tabs. 
These  holes  were  actually  drilled  slightly  oversize  to  accommodate 
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the  pin  easily.  The  larger  pin  size  was  used  for  the  higher  loads 
to  prevent  pin  failure. 

All  of  tne  samples  were  strain-gaged  using  standard  techniques. 
The  same  adhesive  system  and  cure  conditions  used  to  bond  the  end 
tabs  was  employed  during  the  strain-gaging.  Axial  strain  was  meas- 
ured by  two  Mi  cro-Measurements ' strain  gages,  type  F.A-06-125BZ-35Q, 
They  were  bonded  along  the  longitudinal  direction  of  the  sample 
(Figure  14),  one  on  each  side,  and  connected  to  the  signal  condi- 
tioning system  individually  so  that  the  strain  on  each  gage  was 
known  at  all  times  in  order  to  assess  the  effects  of  bending,  if 
present.  Some  of  the  epoxy  specimens,  as  well  as  the  8=0°  and 
o = 90°  glass/epoxy  specimens,  were  strain-gaged  with  the  type 
EA-06-1 OOVA-350  gages.  These  gages  had  an  additional  transverse 
gage  for  measuring  Poisson's  ratio,  \>,  in  the  epoxy  resin,  and 
u1;,  and  in  the  glass/epoxy  for  the  two  respective  fiber  angles. 

Beam  bendi ng  specimens . Several  specimens  of  the  glass/epoxy 
laminates  were  also  prepared  for  four-point  beam  bending  tests. 

These  specimens  were  all  cut  and  prepared  in  a manner  similar  to  the 
tensile  coupons  with  several  exceptions.  The  sample  configuration 
was  slightly  larger,  namely,  7"  x 1"  x 1/8"  and  no  end  tabs  were 
necessary.  For  comparison  purposes  only  the  0/90  and  ±45°  fiber 
angle  layups  were  prepared.  Only  the  bottom  (tension),  central 
section  of  the  beam  was  strain-gaged  for  these  tests  using  type 
EA-06-1 25BZ-350  gages. 
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Filament-wound  case  material  specimens.  Preparation  of  the 
specimens  from  the  forward  dome  required  considerably  more  care 
because  of  the  case  flexibility,  curvature  and  lack  of  support 
during  the  cutting  operation.  The  case  thickness  varies  from  a 
minimum  of  0.08  inches  (eight  plies)  near  the  barrel  section  and 
thrust  termination  (TT)  ports  to  approximately  0.375  Inches  near 
the  igniter  adapter  and  other  areas  of  local  reinforcement 
(Figure  13).  An  area  between  two  of  the  thrust  termination  ports 
was  selected  because  cf  the  uniform  thickness  (0.08  inches)  and 
relatively  small  curvature.  The  interior  of  the  forward  dome 
had  a 1/2-inch  layer  of  rubber  Insulation  bonded  over  the  complete 
dome  area  which  had  to  be  removed  from  the  case  material. 

Two  areas  measuring  10"  x 10"  were  marked  to  be  cut  out.  One 
section  was  to  be  used  for  samples  with  a fiber  angle  of  ±20° 
(meridional  direction)  while  the  other  was  oriented  at  ±70°  (cir- 
cumferential direction).  Each  section  was  cut  out  by  using  a saber 
saw  with  a special  tungsten  carbide  blade.  Before  further  section- 
ing into  useable  sample  sizes  could  be  accomplished,  the  rubber 
insulation  had  to  be  removed.  Removal  was  done  with  a high  speed, 
electric  hand  grinder  in  such  a manner  that  only  a thin  layer  of 
the  insulation  was  left  to  protect  the  fiber  surface.  The  sections 
were  then  marked  and  further  cut  such  that  at  least  four  beam 
specimens,  of  nominal  dimensions  7"  x 1"  x 0.08",  and  two  tensile 
specimens,  measuring  approximately  6"  x 3/4"  x 0.08",  were  obtained 
from  each  section.  Final  trimming  and  smoothing  of  the  edges  was 
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done  with  a high  speed  router. 

All  of  the  samples  were  strain-gaged  by  following  the  tech- 
niques established  previously  for  the  tensile  coupons  and  beams, 
respectively,  with  one  exception.  The  surface  of  the  filament- 
wound  case  material  was  very  irregular  compared  to  the  SCI  plate 
materials,  and  therefore,  the  surface  was  first  coated  with  an 
epoxy  filler  material  and  sanded  smooth  prior  to  gaging.  Several 
of  the  completed  samples  are  shown  in  Figure  15. 

Plate  twist  specimens.  Two  samples  were  used  for  these  tests, 
a ±45°  glass/epoxy  and  6061 -T6  aluminum  plate.  Both  plates  were 
12"  x 12"  x 1/8".  The  glass/epoxy  plate  was  one  of  the  original  SCI 
plates  described  earlier  and  therefore  possessed  midplane  symmetry, 
i.e.,  B . . = 0.  The  plates  were  strain-gaged  on  both  surfaces  at 
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the  intersection  of  their  diagonals  using  a Micro-Measurements' 
EA-06-1 25RD-350  rectangular  rosette  gage  configuration.  The  gages 
were  mounted  with  the  AE-15  adhesive  with  the  two  principal  gages 
oriented  in  the  direction  of  the  diagonals  (Figure  16).  A room 
temperature  cure  for  24  hours  followed  by  postcuring  at  140°F  for 
4 hours  was  used  for  these  samples  in  order  to  reduce  plate  warping. 


Equipment  and  Procedures 


All  of  the  tests  involved  the  use  of  strain-gaged  specimens. 


I 
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Whenever  a test  was  conducted,  a dummy  sample  and  identical  gage 
was  mounted  in  the  same  environment  in  order  to  compensate  for  both 
thermal  expansion  and  humidity  changes.  Strain  gage  measurements 
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were  recorded  using  a B & F Instruments'  Strain  Gage  Acquisition 
System6  and  a modified  Hewlett-Packard  561  Digital  Printer7.  The 
printer  was  modified  to  Include  an  Input  from  a timing  clock.  The 
system  operated  at  a rate  of  two  channels  per  second,  accommodating 
up  to  ten  channels  of  strain  gage  signal  conditioning,  and  provid- 
ing reliable  data  within  a few  seconds  of  the  actual  loading. 

Creep  and  recovery  testers.  The  variation  of  fiber  angle  and 
stress  level  studied  during  the  program  required  the  use  of  several 
creep  testers  in  order  to  provide  the  maximum  amount  of  data  in  a 
reasonable  time  period.  Three  sets  of  creep  testers  were  designed 
and  fabricated  for  the  test  program: 

1.  Multi-station,  dead-weight  test  unit  (500  pounds 
maximum) . 

2.  Single  unit,  lever  arm  system  for  medium  loads 
(2000  pounds  maximum). 

3.  Single  unit,  lever  arm  system  for  heavy  loads 
(6000  pounds  maximum). 

The  multi-station,  dead-weight  creep  tester  consists  of  five 
test  stations,  each  capable  of  sample  loads  up  to  about  500  pounds. 
This  unit,  shown  in  Figure  17,  was  used  for  the  majority  of  the 
testing  where  the  low  load  (stress)  levels  were  of  interest.  For 
the  tensile  coupon  configuration  shown  in  Figure  14  this 
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corresponded  to  a stress  level  of  almost  8000  psl. 

The  dead  weights  were  connected  to  the  lower  grips  which  in 
turn  were  lowered  (creep-load  application)  and  raised  (recovery- 
load  removal)  by  means  of  a scissors  jack  which  controlled  the 
motion  of  the  base  plate  directly  under  the  grips.  Figure  18  shows 
a typical  tensile  coupon  mounted  in  the  tester  without  the  heating 
chamber.  The  specimen  was  pin-connected,  both  at  the  top  and 
bottom,  and  a universal  joint  located  at  the  top  was  used  to  min- 
imize bending.  The  lower  grips  contained  an  oversized  slot  so  that 
the  sample  was  completely  free  of  any  loads  during  the  recovery 
period.  Details  of  the  lower  grip  and  pin-connection  are  also 
shown  in  Figure  18. 

Several  tests  required  load  (stress)  levels  higher  than  500 
pounds  which  in  turn  lead  to  the  design  of  the  lever  arm  creep 
testers.  Five  single  unit,  lever  arm  systems  were  built  such  that 
two  load  ranges  could  be  obtained  by  simply  changing  the  lever 
advantage.  The  basic  unit,  shown  in  Figure  19,  consisted  of  a 
lever  arm,  counter-weight  and  essentially  the  same  grips  as  used  in 
the  dead-weight  creep  tester.  The  unit  shown  in  Figure  19  is  set 
up  for  the  intermediate  load  range  of  2000  pounds  and  has  a 9:1 
lever  advantage.  By  making  a suitable  change  in  lever  arms  and 
connections,  these  same  units  have  a 30:1  lever  advantage  and  can 
be  used  with  wedge-action  grips  to  achieve  6000  pounds.  The  exact 
lever  arm  advantage,  necessary  for  accurate  determination  of  the 
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stress,  was  obtained  through  an  extensive  series  of  calibration 
tests  using  dead  weights  on  the  loading  pan  assembly  and  strain- 
gaged  aluminum  and  steel  specimens  mounted  on  the  grips. 

Cven  with  the  care  taken  to  assure  perfectly  aligned  loading 
pin  holes,  It  was  not  always  possible  to  prevent  some  degree  of 
hole  misalignment.  In  order  to  further  reduce  the  bending  effect, 
a small  load  was  applied  to  the  sample  prior  to  the  actual  test 
(usually  enough  to  Induce  about  50pe).  The  strains  on  both  sides 
of  the  sample  were  recorded  and  then  the  lower  grip  pin  alignment 
was  changed  until  both  strain  readings  were  approximately  equal. 
This  was  accomplished  by  adjusting  two  small  Allen  screws  which 
engaged  the  pin  vertically  in  the  upper  end  of  the  slot  until  the 
strains  were  equivalent.  The  entire  procedure  yielded  strains 
which  were  usually  within  5-10%  of  each  other. 

Temperature  conditioning.  The  creep  and  recovery  tests  as 
well  as  some  of  the  constant  crosshead  rate  tests  were  performed 
at  several  temperatures.  Most  of  the  tests  at  temperatures  between 
-20°F  to  140°F  were  conducted  in  the  large  walk-in  environmental 
rooms  located  In  Texas  A & M University's  McNew  Laboratory.  These 
rooms  had  controls  for  the  humidity  as  well  as  temperature,  there- 
fore providing  a stable  environment  for  long  tern  tests.  Almost 
all  of  the  creep  and  recovery  tests  on  the  glass/epoxy  composites 
were  conducted  in  these  rooms. 

High  temperature  creep  and  recovery  tests  on  the  Shell  58-68R 
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epoxy  resin  were  run  by  enclosing  each  specimen  within  individual 
heaters  such  as  the  ones  shown  in  Figures  20  and  21.  The  cylin- 
drical heating  elements  were  constructed  with  electrical  heating 
tape  wrapped  around  a pre-Formed  wire  cylinder  serving  as  the 
inner  wall.  The  exterior  was  insulated  with  two  layers  ol  asbestos 
tape  and  two  layers  ol  cardboard.  During  testing  tne  ends  of  the 
chambers  were  covered  with  asbestos  plates  to  miriiiiiuo  heat  losses. 
Temperature  control  was  accomplished  with  the  use  of  a small  ther- 
mistor sensor  located  rear  the  specimen  center  and  an  1(1 1.  Industries 
Model  70  Proportional  Controller”.  Temperature  measurement  was 
achieved  by  inserting  a glass  thermometer  through  a small  hole  in 
the  side  of  the  heat  chamber.  These  heating  chambers  were  capable 
of  temperatures  up  to  and  exceeding  260'T  [12]  with  control  to  !2"'F 
or  better.  Previous  calibration  tests  showed  that  the  axial  temp- 
erature gradient  at  the  center  of  the  chambers  was  small. 

Ileam  tests,  lour-point  bonding  tests  were  conducted  in  a 
manner  which  essentially  followed  that  given  in  the  A5TM  1)790-/1 
procedure  for  determining  the  flexural  properties  of  plastics  [102]. 
These  tests  were  conducted  for  both  creep  and  recovery  and  constant 
crosshead  rate  bending  using  the  same  basic  loading  fixture  shown 
in  Figure  22.  The  support  was  a section  of  5-inch  heavy  aluminum 
channel  with  the  edges  sanded  smooth  and  rounded.  The  loading 
points  were  the  Lips  of  a small  section  of  3-inch  steel  channel 
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Figure  20.  Interior  view  of  heat  Figure  21.  Heat  chamber  mounted 

chamber.  around  sample  in  lever 

arm  creep  tester. 


which  had  also  been  sanded.  AH  of  the  loading  and  support  points 
were  covered  with  a 1-inch  wide  strip  of  Teflon  tape  to  minimize 
the  friction. 

The  loading  channel  was  fastened  to  the  bottom  of  one  of  the 
load  pans  in  the  multi-station,  dead  weight  creep  tester  described 
earlier  and  then  lowered  onto  the  upper  surface  of  the  sample.  In 
the  constant  crosshead  rate  tests  the  channel  was  bolted  to  the 
lower  moveable  crosshead  of  an  Instron9  Testing  Machine.  A cross- 
head rate  of  0.2  inches  per  minute  was  used  for  these  tests. 

Plate  twist  tests.  Some  of  the  previous  background  and  int- 
erest in  the  plate  twist  tests  has  already  been  discussed  in  Section 
II  Experimentally,  the  test  is  relatively  simple  to  perform, 
provided  that  the  theoretical  limitations  concerning  large  deflec- 
tions and  localized  loading  conditions  are  considered.  In  the  plate 
twist  test,  a pure  twisting  moment  is  imposed  on  a square  plate  by 
loading  all  four  corners  with  equal  forces.  The  forces  are  perpen- 
dicular to  the  plate  with  those  forces  at  the  first  and  third 
(diagonal)  corners  being  upward  and  the  other  two  forces  downward. 
The  corner  loads  cause  the  square  plate  to  assume  a hyperbolic 
paraboloid  or  saddle-shaped  surface  '183,  184]. 

The  vertical  upward  forces  were  applied  through  the  base 
support.  The  support  consisted  of  a heavy  steel  plate  about 
1/2-inch  thick  with  two  large  steel  blocks  welded  to  it  at  the  ends 
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of  one  diagonal  (see  Figure  23  for  detail).  A hole  was  countersunk 
into  each  block,  making  sure  that  the  distance  between  hole  centers 
was  14  Inches.  A 3/4-inch  diameter  steel  ball  bearing  was  then 
welded  into  the  countersunk  hole  to  provide  a rounded  loading 
support  which  would  reduce  the  effect  of  a concentrated  load. 

The  vertical  downward  forces  were  applied  in  a similar  manner. 

The  ball  oearings  and  blocks  were  welded  14  inches  apart  to  the 
ends  ot  a section  of  3-inch  steel  channel.  A 1/2-inch  diameter 
steel  rod  was  welded  to  the  inside  of  the  channel  , midway  between 
the  ball  bearings.  This  rod  was  used  to  bolt  the  channel  to  the 
upper  crosshead  of  the  Instron  /constant  crosshead  rate  cests)  and 
served  as  a guide  for  a load  bucket  (creep  tests).  All  of  the  ball 
bearing  surfaces  were  covered  with  Teflon  tape  to  minimize  the 
frictional  effects. 

The  displacement  of  the  corners  under  the  downward  forces  was 
monitored  by  two  linear  variable  differential  transformers  (LVOT). 
Placement  of  the  spring-loaded  LVDT's  was  very  critical,  as  was  the 
alignment  of  the  channel  loading  fixture.  They  were  mounted  14 
inches  apart,  directly  under  the  applied  downward  forces.  Cross- 
head deflections  from  the  Instron  tests  were  also  used  and  found  to 
be  in  excellent  agreement  with  the  LVDT  deflections. 

The  constant  crosshead  rate  tests  on  the  Instron  testing 


machine  were  performed  by  mounting  the  base  support  plate  on  a 
compression  load  cell  and  the  loading  channel  to  the  upper  crosshead. 
Dead-weight  creep  tests  were  run  in  a similar  manner  by  carefully 
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weighing  the  entire  loading  channel  assembly  (channel,  load  bucket, 
weights,  etc.)  before  the  test.  The  load  was  lowered  onto  the 
plate  diagonal  as  rapidly  as  possible,  using  several  guide  clamps 
to  assure  the  alignment  of  the  applied  load. 

Experimental  Considerations 
Bending  and  Grip  Effects 

The  use  of  off-angle  tensile  specimens  (0°  < o < 90°)  under 
uniaxial  loads  creates  bending  moment  and  shear  stress  distribu- 
tions along  the  length  of  the  specimen  due  to  the  existence  of 
rigid  clamps  which  are  prevented  from  rotating.  The  clamping  con- 
dition was  Investigated  in  [167,185]  and  it  was  shown  in  [167]  that 
the  difference  between  the  actual  compliance  and  the  compliance 
defined  by  the  centerline  strain  and  the  axial  stress  (axial  force 
divided  by  the  cross-sectional  area)  was  less  than  1%  for  an 
E-glass/epoxy  composite.  It  was  pointed  out  that  the  effect  is  a 
strong  function  of  the  geometry  of  the  tensile  coupon  as  well  as 
its  anisotropy.  For  the  materials  and  specimens  employed  in  this 
program  this  condition  was  found  to  be  negligible. 

Strain  Gage  Heating 

The  use  of  bonded  foil  resistance  strain  gages  has  been  found 
to  lead  to  significant  gage-specimen  interaction  if  the  material 
is  near  the  softening  or  transition  temperatures . High  resistance 
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gages  of  350  ohms  were  used  to  minimize  the  local  heating.  An 
earlier  study  [60]  conducted  on  the  Shell  58-68R  epoxy  resin  showed 
that  local  heating  for  these  same  gages  was  on  the  order  of  3-5°F 
for  a 2.5  volt  bridge  excitation  voltage.  The  gages  used  had  a 
lower  power  output  than  some  other  available  gages  such  as  the  more 
common  120  ohm  foil  gages. 


Humidity  Effects 

Humidity  has  been  shown  to  affect  polymeric  materials  In  much 
the  same  way  as  temperature  [11].  An  increase  in  relative  humidity, 
or  water  content,  is  equivalent  to  an  Increase  in  temperature. 
Several  experiments  were  conducted  on  an  E-glass/epoxy  composite 
[186]wh1ch  showed  that  humidity  control  should  be  considered  during 
any  characterization  test  program.  In  view  of  this  fact,  the 
samples  were  always  allowed  to  reach  an  equilibrium  in  the  test 
environment  for  a period  of  at  least  24  hours.  The  humidity  in  the 
environmental  chambers  was  held  constant  during  the  experiments. 
However,  it  was  not  possible  to  maintain  a given  relative  humidity 
level  at  all  of  the  temperatures  tested  because  of  equipment 
limitations. 


Softening  and  Transition  Temperatures 


The  softening  temperature  of  a material  has  generally  been 
defined  as  the  approximate  value  at  which  a rapid  increase  In  creep 
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rate  begins  as  the  specimen  Is  slowly  heated.  Above  the  softening 
temperature  the  composite  becomes  extremely  soft  and  loses  Its 
practical  structural  value.  One  way  to  obtain  an  estimate  of  this 
temperature  has  been  to  conduct  a low  stress  level  creep  test  on 
an  off-angle  specimen  while  slowly  Increasing  the  temperature 
[9,  12].  The  point  or  region  at  which  the  creep  rate  rapidly  in- 
creases is  then  defined  as  the  "softening  temperature". 

Another  temperature,  the  glass  transition  temperature,  T , 

, y 

represents  the  point  at  which  the  physical  mechanism  of  deformation 
within  the  polymer  changes  due  to  a significant  change  in  free 
volume  [52].  As  an  epoxy  resin  or  composite  system  is  heated 
through  the  transition  temperature,  the  expansion  coefficient  has 
been  found  to  increase  by  a factor  of  two  or  three.  The  Tg  repre- 
sents a more  physical  change  in  the  polymeric  system  and  is  always 
below  the  softening  temperature.  It  was  necessary  to  determine  the 
region  where  the  transition  took  place  in  order  to  define  the  max- 
imum temperature  for  the  characterization  tests. 

The  thermal  expansion  behavior  of  the  test  materials  (Shell 
58-68K  epoxy  resin,  S -901  glass/epoxy  resin  and  aluminum)  was  deter- 
mined with  a DuPont  Thermal -Mechanical  Analyzer  (TMA)10.  The  TMA 
consists  of  a fused  quartz  holding  tube  and  an  LVDT  probe  to  measure 
the  expansion  characteristics  of  small  samples  under  variable  heat- 
ing and  cooling  rates.  A sample  of  each  material  measuring 
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approximately  1/4"  x 1/4"  x 1/8"  was  cut  from  a representative 
section  of  the  plates.  The  specimen  was  a unidirectional  laminae 
with  the  expansion  measured  In  the  direction  of  the  plate  thickness. 
The  aluminum  sample  was  used  primarily  as  a reference  standard. 

Each  sample  was  thermally  cycled  (heating  followed  by  cooling) 
several  times  at  a rate  of  10°C/minute.  The  results  of  these  tests 
are  shown  In  Figures  24  through  26. 

All  of  the  materials  exhibited  a thermal  hysteresis  loop  which 
was  believed  to  have  resulted  from  the  high  heating  rate  used. 

Even  with  the  small  sample  dimensions  employed,  the  temperature  of 
the  specimen  was  not  expected  to  be  completely  uniform  at  that  rate. 
However,  the  tests  showed  several  distinct  trends  which  were  con- 
sidered Important,  The  hysteresis  loop  in  the  epoxy  resin  and 
glass/epoxy  composite  can  be  attributed  to  some  degree  to  the 
shrinkage  and  weight  loss  associated  with  the  Initial  moisture 
content  of  the  materials.  Freeman  and  Campbell  [187] reported  the 
same  effects  on  epoxy  resin  and  graphite/epoxy  composites  using  a 
quartz  tube  dilatometer  and  slow  heating  (cooling)  rates.  After 
several  cycles  above  the  boiling  point  of  water,  212°F,  the  material 
appeared  to  reach  a stable  state.  It  is  not  known  at  this  time  if 
the  process  Is  reversible  or  If  any  permanent  damage  has  been  done. 

In  the  glass/epoxy  composite  (Figure  26),  the  effect  was  more 
pronounced.  The  relaxation  of  thermal  stresses  (originally  devel- 
oped during  curing)  and  drying  out  could  possibly  explain  the 
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TEMPERATURE  (°F) 

Thermal  expansion  behavior  of  aluminum  reference  material. 


Figure  25-  Theraa'i  expansion  of  Shell  58-68R  epoxy  resin. 


relatively  large  changes  in  the  first  and  second  thermal  cycles. 
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The  third  cycle  appeared  to  be  very  stable.  The  T^  for  these  mat- 
erials appears  to  be  around  180°F  or  so,  while  the  softening  tem- 
erature  was  at  a higher  temperature.  As  a result,  all  of  the 
characterization  tests  were  conducted  below  160°F  after  drying  the 
samples  for  at  least  24  hours.  These  effects  are  discussed  in 
Section  V with  regard  to  composite  design  and  analysis. 

Preliminary  Characterization  Tests 

Several  preliminary  tests  were  conducted  to  determine  some  of 
the  basic  properties  of  the  S-901  glass  and  glass/epoxy  composite 
system.  These  properties  were  used  as  a guide  in  establishing 
various  test  parameters. 

Fiber  Properties 

The  determination  of  the  fiber  properties  by  SCI  was  considered 
necessary  to  Insure  that  the  S-901  glass  used  in  the  laminates  was 
within  the  normal  specifications.  Consequently,  SCI  conducted 
several  strength  tests  on  the  20-end,  S-901  glass  roving  taken 
from  each  of  the  three  spools  used  in  the  fabrication.  These  strand 
tensile  tests  indicated  that  the  tensile  strength  was  nominally 
525,000  psi,  well  above  the  minimum  specification  of  500,000  psi 
[181]. 

The  specific  gravity,  tensile  modulus  and  Poisson's  ratio  of 
the  glass  fibers  were  not  determined  by  SCI  since  this  data  has  been 


reported  previously  [179,  188] 


- 2.465 


Ef  = 12.4  x 106  psi 


vf  = 0.22 


The  value  used  for  the  specific  gravity  represents  an  average 
of  the  reported  values.  These  properties  are  used  in  subsequent 
analysis  of  the  overall  composite  behavior  as  discussed  in 
Section  IV. 


Fiber  Volumetric  Content 

The  fiber  content  (by  volume),  vf,  of  a composite  material  , 
strongly  influences  the  overall  mechanical  and  physical  behavior  of 
the  composite  as  shown  by  the  Halpin-Tsai  Equation  (119).  As  a 
result  of  its  importance  in  the  analysis  of  the  data,  glass  volu- 
metric tests  were  run  at  SCI  and  TAMU.  SCI  conducted  gravimetric 
tests  on  the  edge  trimmings  taken  from  each  panel  and  determined 
the  resin  content,  voids  and  glass  volume.  Void  determinations 
were  generally  found  to  be  less  than  0.1%  and  were  discontinued 
after  the  first  four  panels  [181],  The  SCI  data  is  shown  in  Table  3 
along  with  the  values  determined  by  TAMU.  Glass  fiber  content 
values  determined  by  SCI  are  considered  questionable  with  regard  to 
being  representative  of  the  panel  since  the  samples  were  edge 
trimmings. 
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Fiber  content  values  more  representative  of  the  actual  test 
samples  were  determined  from  samples  cut  from  the  interior  of  the 
panels  (i.e.,  away  from  the  edges).  These  tests  were  conducted  on 
samples  which  had  been  tested  and  were  not  needed  for  further  test- 
ing. Specimens  were  generally  cut  from  the  center  sections  of  the 
tensile  coupons  after  removing  the  strain  gage.  E.-ch  sample  was 
carefully  weighed  and  then  heated  in  a muffle  furnace  at  600°F  for 
4 hours  and  then  1200°F  for  24  hours  [179].  After  this  sequence  the 
resin  was  completely  burned  off.  The  sample  was  removed  to  a 
dessicator  where  It  cooled  to  ambient  temperature.  The  sample  was 
weighed  again  and  the  fiber  content  determined  from  the  values  of 
specific  gravity  and  measured  weights.  These  values  are  shown  in 
Table  3 and  were  considered  to  be  more  representative  of  the 
actual  glass/epoxy  composite  use  in  the  program. 

The  fiber  content  of  samples  taken  from  the  forward  dome  of 
the  Minuteman  III  case  was  determined  after  removing  the  rubber 
insulation  material  and  the  adjacent  glass/epoxy  lamina.  Using  the 
same  procedure  as  before,  the  fiber  content  shown  in  Table  3 was 
found  to  be  very  close  to  the  SCI  panels.  This  makes  the  compar- 
ison of  plate  and  case  properties  more  valid  since  the  volumetric 
contents  of  the  constituents  are  almost  exactly  the  same.  .A  value 
of  0.616,  representing  the  average  of  the  TAMU  fiber  volume  content 
on  the  panel  specimens,  was  used  for  analytical  purposes. 
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Tensile  Coupon  Stress-Strain  Behavior 

In  order  to  establish  the  creep  stress  levels  and  approximate 
strains  for  the  test  program,  constant  strain  rate  tests  were  con- 
ducted on  the  various  materials  which  were  used  during  the  program. 
These  tests  were  run  using  an  Instron  tensile  tester  at  a crosshead 
rate  of  0.02  inches  per  minute  at  temperatures  of  75°F  and  140°F. 
These  temperatures  essentially  represent  the  range  investigated 
curing  the  program. 

The  rectangular  tensile  coupons  described  in  an  earlier  sub- 
section were  used  for  the  tests.  Wedge-action  grips  were  used 
rather  than  the  pin-connected  creep  grips  since  the  samples  were 
taken  to  complete  failure.  It  should  be  noted  that  this  particular 
specimen  design  provides  excellent  mechanical  characterization 
data  for  a uniaxial  stress  field:  however,  because  the  sppclmen 
does  not  have  a reduced  test  section,  the  specimen  will  usually 
fail  prematurely  near  the  grips  or  end  tabs.  Consequently,  these 
specimens  were  not  used  to  generate  failure  stress-strain  informa- 
tion and  the  constant  rate  tests  employed  here  served  only  as  a 
guide  for  establishing  test  limit  conditions. 

The  stress-strain  behavior  of  the  Shell  58-68R  epoxy  resin  and 
S-90T  glass/epoxy  resin  is  shown  in  Figures  27  through  37.  The  in- 
fluence of  temperature  is  evident,  and  the  possible  presence  of 
nonlinear,  viscoelastic  material  behavior  is  particularly  notice- 
able for  the  off-angle  glass/epoxy  orientations.  Strain  data  past 
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Figure  36.  Stress-strain  behavior  of  laminated  glass/epoxy  (e 
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2%  was  not  available  for  the  0°/90°  and  ±45°  fiber  angles  because 
of  equipment  limitations,  but  the  extrapolated  behavior  is  shown 
as  a dotted  line  in  the  respective  figures. 

Constant  Crosshead  Rate  Tests 

Uniaxial  Tension 

Using  the  tensile  coupons  with  the  5°  wedge-shaped  end  tabs, 
a series  of  constant  crosshead  rate  tests  were  conducted  on  several 
45°  and  ±45°  glass/epoxy  samples.  These  tests  were  carried  out  at 
temperatures  of  75°F  and  140°F  using  the  Instron  and  the  environ- 
mental rooms.  Each  sample  was  tested  at  a crosshead  rate  of  0.02 
inches  per  minute  to  a stress  level  of  6000  psi  and  then  unloaded 
to  zero  stress  at  a crosshead  rate  of  0.2  inches  per  minute.  Three 
cycles  were  run  on  each  sample  after  waiting  approximately  5 min- 
utes between  cycles.  All  of  the  tests  were  conducted  with  the 
Instron  wedge-action  grips.  A second  set  of  5°  wedges,  placed  in 
the  grips  with  a piece  of  Teflon  tape  facing  the  sample  end  tabs, 
permitted  the  sample  to  hang  free  of  obstructions  upon  unloading. 
This  technique  gave  more  consistent  results  and  allowed  the  spec- 
imen to  recover  unhindered  (zero  stress). 

Four-Point  Beam  Bending 

All  of  the  beam  bending  tests  were  performed  at  75° F using  the 


Instron  tensile  tester  and  a crosshead  rate  of  0.2  inches 
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per  minute.  Glass/epoxy  specimens  with  fiber  angles  of  ±45,J  and 
0°/90°  and  filament-wound  specimens  with  nominal  fiber  angles  of 
±20°  and  +70°  were  used  in  the  tests.  Since  only  one  side  of  the 
specimen  was  strain-gaged,  the  side  with  the  gage  was  mounted  face 
downward  in  order  to  record  the  tension  strain.  In  the  case  of  the 
filament-wound  specimen  this  corresponded  to  the  surface  without 
the  rubber  insulation  material.  The  temperature  compensating 
specimens  were  mounted  on  the  end  of  the  support  fixture  in  the 
same  position  so  as  to  compensate  for  creep  effects  due  to  the 
sample  weight.  This  procedure  was  necessary  primarily  for  the 
thinner  filament-wound  case  materials  which  were  more  flexible. 

The  beam  samples  were  loaded  to  a preselected  maximum  stress 
(load)  and  then  unloaded  completely,  always  at  the  same  crosshead 
rate.  Strain  and  load  information  were  recorded  continuously 
:ring  this  period  and  at  specified  time  intervals  following  the 
complete  unloading.  A recovery  period  of  approximately  three  times 
the  complete  cycle  time  was  allowed  tefore  starting  the  next  load 
cycle.  Each  maximum  stress  (load)  cycle  was  repeated  in  the  same 
manner  until  two  or  three  cycles  were  completed.  The  beam  was  then 
subjected  to  another  preselected,  higher,  maximum  stress  (load) 
cycle  following  the  previous  procedure. 

The  maximum  stress  (load)  levels  for  each  of  the  beams  tested 
are  given  in  Table  4 in  terms  of  the  bending  moment,  Mx  (in-lbs), 
which  existed  at  the  center  section  of  the  beam. 
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Table  4.  Summary  of  Maximum  Bending  Moment, 
Mx  (In-lbs),  for  Four-Point  Bending 
Tests  (Constant  Crosshead  Rate) 


Fiber 

Angle 

±45° 

0°/90° 

±20° 

±70° 

(case) 

(case) 

10 

10 

5 

3 

25 

50 

10 

6 

m¥ 

X 

40 

90 

25 

70 

150 

Plate  Twist  Test 


Constant  crosshead  rate  tests  conducted  at  75°F  and  0.1  inches 
per  minute  were  carried  out  on  a 6061 -T 6 aluminum  and  a i45° 
glass/epoxy  plate.  During  the  tests  the  outputs  from  ti.  Instron 
load  cell,  the  surface  strain  gages  and  the  corner  LVDT's  were  mon- 
itored continuously.  The  aluminum  plate  was  used  primarily  as  an 
isotropic  reference  material  and  was  loaded  only  within  its  linear 
elastic  range. 

Each  test  consisted  of  loading  and  unloading  the  plate  :o  a 
preselected  maximum  load  level  at  the  same  crosshead  rate.  After 
unloading,  the  plate  was  allowed  to  recover  for  a period  of  approx- 
imately three  times  the  complete  cycle  time  before  starting  the 
next  cycle.  In  the  case  of  the  aluminum,  the  plate  was  cycled 
twice  to  a load  level  of  about  50  pounds.  This  corresponded  to  a 
twisting  moment,  MXy,  of  almost  25  in-lbs.  The  ±45° 
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glass/epoxy  plate  was  cycled  in  a similar  manner.  Each  cycle  was 

repeated  three  times  before  increasing  the  maximum  load  level. 

Three  load  levels  were  studied,  corresponding  to  twisting  moments, 

M , of  25,  42  and' 53  in-lbs  in  succession, 
xy 

Single  and  Multiple  Cycle  Creep  and  Recovery  Tests 
Uniaxial  Tension 


In  order  to  characterize  the  epoxy  resin  and  glass/epoxy 
composites,  both  single  and  multiple  cycle  creep  and  recovery  tests 
were  conducted.  The  creep  test  consisted  of  rapidly  applying  a 
constant  uniaxial  load  to  the  tensile  coupon  and  measuring  the 
resultant  strain  for  one  hour.  This  test  was  immediately  followed 
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After  completion  of  one  cycle  (creep  and  recovery),  the 
total  strain  was  recorded,  the  strain  electrically  rezeroed  and  the 
next  cycle  was  Initiated.  The  strain  was  therefore  referred  to  the 
specimen  length  which  existed  at  the  end  of  each  cycle  and  little 
time  was  permitted  for  bond  reformation  between  flaw  surfaces. 

This  last  point  was  considered  important  since  the  earlier  work  on 
the  E-glass/epoxy  Indicated  that  long  recovery  periods  created  a 
necessity  for  repeating  the  mechanical  conditioning  due  to  re- 
healing. 

Tests  on  the  Shell  58-68R  epoxy  resin  matrix  material  were 
conducted  over  a temperature  range  from  20° F to  160°F  for  two  cycles 
of  creep  and  recovery.  Several  stress  levels  ranging  between 
200  psi  and  2000  psi  were  investigated  during  this  phase  of  the 
program.  Most  of  the  creep  and  recovery  tests  carried  out  on  the 
S-901  glass/Shell  58-68R  epoxy  resin  were  performed  at  75°F  and 
140°F  with  the  exception  of  a few  tests  conducted  on  45°  and  ±45° 
glass/epoxy  specimens  at  20°F  to  establish  the  validity  of  time- 
temperature  superposition  (Tables  5 through  7).  Three  cycles  of 
creep  and  recovery  were  completed  for  each  sample.  As  many  as 
five  different  stress  levels  were  studied  for  each  fiber  angle 
studied,  depending  on  the  degree  of  nonlinearity  observed  in  the 
constant  crosshead  rate  stress-strain  behavior.  Most  of  the  work 
centered  around  the  20°,  45°,  ±30°  and  ±45°  glass/epoxy  composite 
layups . 
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Table  5.  Summary  of  Uniaxial  Creep  and 
Recovery  Tests  (T  = 20°F) 


Fiber  Angle 

(e) 

Stress 

(psi) 

Number  of 
Cycl  es 

Epoxy 

600 

2 

1200 

2 

0 

7000 

2 

45 

520 

3 

3010 

O 

sJ 

90 

500 

3 

+45 

520 

3 

3270 

3 
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Table  6.  Summary  of  Uniaxial  Creep  and 
Recovery  Tests  (T  = 75°F) 


Fiber  Angle 

Stress 

Number  of 

(0) 

(psl) 

Cycles 

Epoxy 

200 

2 

400 

2 

600 

2 

1200 

2 

2000 

3 

0 

7000 

3 

42,000 

3 

20 

2000 

2 

4800 

3 

8700 

3 

10,500 

3 

12,500 

3 

45 

500 

3 

3000 

3 

5000 

2 

6400 

3 

8500 

1 

90 

550 

3 

2000 

3 

4000 

3 
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Table  7.  Summary  of  Uniaxial  Creep  and 
Recovery  Tests  (T  e 140°F) 


Fiber  Angle 

Stress 

Number  of 

(e) 

(psi) 

Cycles 

Epoxy 

200 

2 

(145°F) 

400 

2 

600 

2 

1200 

2 

2200 

2 

0 

42,600 

3 

20 

2000 

3 

4800 

3 

8500 

2 

10,500 

2 

45 

500 

3 

3000 

2 

5000 

3 

6500 

1 

90 

500 

2 

1500 

2 

2400 

2 

I 

I 


I 


I 

1 
t * 


u 


147 


Creep  and  recovery  tests  were  also  conducted  on  specimens  of 
the  filament-wound  case  material  at  75°F.  A fiber  angle  of  ±70° 
was  used  for  this  series  of  tests  because  of  the  lower  specimen 
curvature  in  the  circumferential  direction.  Stress  levels  of 
300  psl  and  1000  psi  were  used  for  this  series  of  tests.  Three 
creep  and  recovery  cycles  were  run  on  each  specimen  at  the  selected 
stress  level. 

Four-Point  Beam  Bending 

Using  basically  the  same  procedures  applied  in  the  constant 
crosshead  rate  tests,  three  cycles  of  creep  and  recovery  were  run 
on  samples  of  the  same  fiber  orientation  at  75°F.  The  maximum 
creep  stress  (load)  levels  for  each  of  the  beams  tested  are  given 
in  Table  8 In  terms  of  the  bending  moment,  M (in-lbs),  which 

A 

existed  at  the  center  section  of  the  beam. 


Table  8.  Summary  of  Maximum  Bending  Moment, 
Mx  (in-lbs),  for  Four-Point  Bending 
Tests  (Creep  and  Recovery) 


Fiber 

Angle 

±45° 

0790° 

±20° 

±70° 

(case) 

(case) 

M 

1.25 

1.25 

3.0 

1.5 

X 

55.5 

85.0 

10.5 

5.5 
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Plate  Twist  Test 


Creep  and  recovery  tests  on  the  +45°  glass/epoxy  plate  were 

conducted  at  load  levels  corresponding  to  a maximum  twisting  moment, 

M , of  3.0  and  16.7  in-lbs  successively.  Each  cycle  was  re- 
xy 

peated  three  times  before  increasing  the  maximum  load  level.  The 
strain  output  from  the  surface  gages  and  the  corner  deflections 
from  the  LVDT's  were  monitored  on  a logarithmic  time  scale  during 
the  creep  and  recovery  phases.  The  loads  were  set  in  place  on  the 
surface  of  the  plate  manually,  using  the  guide  fixtures  in  Figure 
to  assure  alignment.  Recorded  data  was  taken  only  after  the  loads 
were  completely  aligned;  usually  on  the  order  of  20-30  seconds 
after  initial  application. 
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SECTION  IV 

DATA  REDUCTION  AND  ANALYSIS  OF  RESULTS 
Uniaxial  Creep  and  Recovery  of  Epoxy  Matrix 


Determination  of  Material  Constants 


In  Section  II  several  analytical  forms  used  to  represent  the 
compliance  of  viscoelastic  materials  were  discussed.  The  power 
law  form  given  by  Equation  (50)  has  been  found  to  represent  the 
behavior  of  many  rigid  plastics  with  and  without  reinforcement 
[73-76].  The  compliance  given  by  Equation  (50),  repeated  here  for 
convenience,  is 

D ( t ) = Dq  + Djtn  (50) 


where  DQ,  Dj  and  n are  positive  constants  which  are  independent 
of  time.  The  values  of  the  material  constants  may  reflect  temper- 
ature dependence  if  the  material  is  a TCM,  as  was  noted  earlier. 

In  the  case  of  fiber-reinforced  plastic  composites  where  the 
fiber  is  considerably  stiffer  than  the  plastic  matrix,  the  net 
compliance  is  relatively  small  and  the  initial  compliance,  DQ,  is 
not  easily  discernable  from  experimental  data.  This  also  is  true 
of  the  epoxy  matrix  without  reinforcement  when  the  epoxy  is  rela- 
tively stiff,  i.e.,  Eq  = 0.5(106)  psi . As  a result,  the  constants 
defining  the  power  law  cannot  be  found  accurately  from  short-term 
data. 
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Lou  and  Schapery  [9]  developed  an  analytical  procedure  for 
evaluating  the  constants  by  employing  the  power  law  form  of  the 
compliance  with  experimental  data  from  both  creep  and  recovery 
tests.  The  procedure  uses  the  superposition  principle  in  order  to 
represent  the  recovery  compliance  defined  by  (see  Figure  1) 

D (t)  : (130) 

o 

which  results  in 

Dr(t)  = Dj  t' n [(1+A)n  - An]  (131) 

where 

A h r - • (70) 

t 

Equation  (131)  is  of  the  same  form  as  Equation  (69)  when 
9j  = g2  = aQ  = 1 for  the  linear  viscoelastic  region  and  t1  is  the 
time  at  which  the  recovery  initiates.  When  Equation  (131)  is 
plotted  on  log-log  scales,  the  result  is  a standardized  curve 
whose  shape  is  dependent  only  on  the  value  of  n through  the  function 
in  the  brackets.  The  exponent  n may  be  determined  by  plotting 

Dref  = (1+A)n  - An  (132) 

for  several  values  of  n (typically  0 < n < 0.5)  and  then  overlaying 
these  curves  on  the  experimental  recovery  data  curves  to  select  an 


n value  as  shown  In  Figure  38.  The  theoretical  curve  that  fits  all 
of  the  data  most  accurately  was  found  to  be  n = 0.19.  This  value 
appeared  to  be  independent  of  both  stress  and  temperature. 

There  are  two  methods  which  may  be  used  to  evaluate  Dq  and  Dj 
once  the  value  of  n is  known.  The  first  uses  recovery  data  plotted 
on  log-log  scales  and  Equation  (131).  The  vertical  shift,  along  the 
compliance,  is  equal  to  log  (D^’0).  With  the  values  of  n and  t* 
already  known,  the  shift  value  may  be  used  to  determine  D . DQ  may 
be  found  by  choosing  any  point  on  the  creep  curve  and  using  Equation 
(50).  The  second  method  for  evaluating  Dc  and  involves  the  use 
of  Equation  (50)  and  any  two  points  from  the  creep  data,  which  are 
then  used  to  yield  two  simultaneous  equations  with  DQ  and  D1  as 
the  unknowns. 

There  are  several  points  which  must  be  made  with  regard  to 
these  procedures,  all  of  which  are  designed  to  provide  accurate 
results  in  the  analysis  of  the  creep  and  recovery  data.  First  of 
all,  it  should  be  noted  that,  at  least  for  cases  where  the  amount 
of  net  creep  strain  is  relatively  small,  one  should  always  work 
with  data  points  taken  from  a smooth  curve  drawn  through  the  exper- 
imental data  rather  than  the  Individual  data  points.  This  is  part- 
icularly true  when  the  recording  equipment  can  only  record  strain  in 
units  of  1 , 2 or  5 microstrain,  which  at  low  stress  levels  is  in- 
adequate to  show  the  changes  continuously  over  the  time  scale.  As 
a result,  the  data  often  appears  in  quantum,  or  discrete,  jumps. 


38.  Evaluation  of  n usinc  the  recovery  compliance  for  the 
Shell  58-68R  epoxy. 


Determination  by  the  second  method,  solution  ;f  two  simultaneous 
equations,  was  found  to  be  easier  since  the  epoxy  behaved  in  a 
linear  viscoelastic  manner.  In  order  to  accurately  determine  Dq 
and  Dj , creep  data  were  taken  at  21  time  points,  approximately 
evenly  spaced  on  a log  time  scale.  A computer  program  was  developed 
to  solve  pairs  of  simultaneous  equations  taken  at  several  time 
points.  This  procedure  was  easier  and  less  time  consuming  than  the 
method  which  involved  recovery  data. 

The  values  of  Dq1  Di  and  n may  be  checked  by  rewritting 
Equation  (50)  in  the  form 

D(t)  - DQ  « Djtn  (133) 

which,  if  plotted  on  log-log  scales  should  result  in  a straight 
line  with  a slope  of  n and  a value  of  Dj  at  t = 1. 

Effect  of  Stress  and  Temperature 

Creep  and  recovery  tests  were  conducted  over  a temperature 
range  from  20°F  to  160°F  at  stress  levels  up  to  approximately 
3000  psi.  In  all  cases  the  epoxy  satisfied  both  the  homogeneity 
and  superposition  requirements  for  linearity  up  to  the  maximum 
stress  level  studied.  The  ultimate  stress  for  the  epoxy  has  gen- 
erally been  reported  to  be  in  the  9,000-11,000  psi  range  [179,  180] 
although  the  present  tensile  coupon  design  did  not  allow  us  to 
rpach  these  levels  prior  to  failure  near  the  grips.  Therefore, 
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it  is  possible  that  nonlinearity  is  present  at  the  higher  stress 
levels.  The  exponent  in  .he  power  law,  n,  was  found  to  be  equal 
to  0.19  under  all  test  conditions  investigated. 

The  creep  compliances  for  several  temperatures  are  shown  in 
Figure  39  for  the  epoxy  resin  using  data  taken  at  various  stress 
levels.  The  DQ  values  determined  from  the  solution  of  simultaneous 
equations,  i.e..  Equation  (50),  are  given  in  Table  9.  Also  given 
are  the  values  of  Dj  obtained  in  the  same  manner  and  the  values 
of  Dj  obtained  from  the  recovery  data,  i.e.,  Equation  (131). 

Table  9.  Power  Law  Constants  for  Shell  58-68R  Epoxy  Resin 


Temperature 

(°F) 

Dq  x 10“c 
(psi-1) 

0,  x 10~6 
(from  creep  data) 

D,  x 10"& 

(from  recovery  data) 

20 

1.726 

.019 

.025 

75 

1,883 

.069 

.069 

no 

1.934 

.124 

.142 

130 

2.022 

.161 

.171 

145 

2.012 

.198 

.201 

160 

2.070 

.232 

i 

-C* 

The  temperature  dependence  of  DQ  can  be  seen  in  Figure  40  to 
increase  linearly  with  temperature.  As  a result  the  epoxy  must  be 
treated  as  a TCM-2  where  D0  = DQ(T).  This  is  particularly  impor- 
tant when  conditions  of  transient  temperature  exist  in  a structure. 


68R  epoxy 
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In  all  future  analysis,  the  value  of  DQ  for  a specific  temperature 
was  obtained  from  the  curve  drawn  In  Figure  40  and  not  the  individ- 
ual data  points.  It  should  also  be  noted  that  Watkins  [60]  shows  a 
somewhat  different  relationship  for  D (T)  in  some  earlier  work. 

The  difference  appears  to  be  the  result  of  only  a limited  amount 
of  preliminary  data  at  the  time  the  earlier  work  was  done.  The 
preliminary  work  was  done  on  several  samples  which  were  not  ade- 
quately temperature  compensated;  however,  the  general  trend  is  the 
same  in  both  cases. 

After  subtracting  the  initial  compliance,  D , the  net  creep 
compliance,  aD,  given  by  Equation  (133)  was  determined  and  is 
shown  in  Figure  41.  The  data  appears  to  validate  the  n value  of 
0.19  found  from  recovery  data.  The  recovery  compliances  given  by 
Equation  (130)  are  also  shown  in  Figure  42  for  the  epoxy  as  a 
function  of  the  nondimenslonal  time,  X. 

The  values  of  the  shift  factors,  a^  and  ay,  were  also  found 
from  the  isothermal  creep  data.  Actually  it  is  not  possible  to 
separately  determine  a,,  and  ay  for  the  power  law  material  using 
isothermal  creep  data  alone.  However,  Watkins  [60]  found  that  the 
value  of  aQ  = 1 fit  the  data  taken  from  transient  temperature  tests 
on  the  same  epoxy  over  a large  portion  (70°F  <_ T <^160°F)  of  the 
temperature  range  investigated  in  this  study.  As  a result,  the 
evaluation  of  the  time-temperature  shift  factor,  ay,  can  be  easily 
determined  by  employing  a graphical  shift  along  the  time  axis  of 


LOG  AD  (psi 


O 160  F 
• 145°F 
□ 130°F 
■ 110°F 


LOG  TIME  (minutes) 

Figure  41.  Net  creep  compliance,  AD,  for  Shell  58-68R 
epoxy  at  different  temperatures. 
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the  net  creep  compliance,  aD  (Figure  41),  as  well  as  evaluating 
the  ratios  of  Dj  at  different  temperatures. 

The  graphical  shifting  of  the  isothermal  net  creep  compliances 
to  obtain  a smooth  master  curve  is  shown  by  Figure  43  where  only 
a few  of  the  data  points  have  been  shown  for  illustration  purposes. 
The  second  procedure  involves  the  use  of  the  Dj  values  shown  in 
Table  9 for  the  various  temperatures.  Since  the  values  of  Dj  for 
both  creep  and  recovery  are  close,  only  one  set,  the  recovery  Dj 
values,  was  used  for  reduction  purposes.  The  value  of  D: , as  used 
in  Equation  (130)  in  this  instance,  reflects  temperature  dependence. 
In  order  to  place  the  temperature  dependence  where  it  should  be, 
in  the  reduced  time,  it  should  be  realized  that  the  values  given 
in  Table  9 are  really 

V 

D,  : — (134) 

aT 

where  is  a constant.  Thus,  by  choosing  D /at  75°F  as  a refer- 
ence value  (ij  si),  the  ratios  of  at  other  temperatures  to 
at  75°F  will  yield  the  value  of  a-j.n(T)  where  aTn(75)  has  been  given 
the  value  of  unity. 

The  shift  factors  found  by  these  procedures  are  shown  in 
Table  10  and  Figures  44  and  45.  Good  agreement  between  the  two 
methods  is  indicated.  The  shift  factor  is  plotted  against  the 
inverse  temperature  (Figure  44)  in  order  to  determine  the  activa- 
tion energy  using  Equation  (24).  The  activation  energy  was  found 


1 


TEMPERATURE  (c^) 

Figure  45.  Temperature  dependence  of  the  shift  factor,  3j,  for  Shell  58-68R  ep< 
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to  be  30  K cal/g-mole  which  Is  typical  of  that  generally  reported 
for  epoxy  resin  materials  [11]. 


Table  10.  fihift  Factors,  ay,  for  Shell  58-68R  Epoxy  Resin 


Temperature  (°F) 

log  aT 

(from  aD  shift) 

log  aT 

(from  Dj  ratios) 

20 

2.95 

2.32 

75 

0 

0 

no 

-1  .34 

-1.65 

130 

-1 .94 

-2.07 

145 

-2.41 

-2.44 

160 

-2.77 

-2.92 

» t 


The  Poisson's  ratio,  v,  of  the  epoxy  was  measured  at  20,  75 
and  140°F  at  several  stress  levels  below  600  psi.  No  appreciable 
time  dependence  was  noted  although  some  of  the  tests  lasted  as 
long  as  60  minutes.  A value  of  v = 0.393,  which  represents  the 
average  of  tests  at  all  three  temperatures , will  be  used  for  all 
future  analysis.  There  was  no  strong  temperature  dependence  of 
Poisson's  ratio  and  the  value  quoted  represents  a ±5%  variation  of 
the  experimental  values. 
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Multiple  Cycling  Effects 


All  of  the  creep  and  recovery  tests  on  the  epoxy  resin  were 
conducted  for  at  least  two  cycles  in  order  to  assess  any  effects 
from  multiple  cycling.  The  data  from  two  cycles  are  shown  in 
Figures  46,  47  and  48  for  the  creep  compliance,  net  creep  com- 
pliance and  recovery  compliance,  respectively.  Two  different 
stress  levels  are  shown  for  comparison;  however,  as  noted  before, 
there  does  not  appear  to  be  any  strong  nonlinearity  due  to  stress 
present  at  the  levels  Investigated.  The  main  difference  in  the 
two  creep  compliance  curves  shown  in  Figure  46  for  the  different 
stress  levels  represents  less  than  0.6%  which  is  probably  sample 
variation.  When  the  strain  during  the  creep  portion  of  the  test 
is  referenced  to  the  sample  length  existing  just  prior  to  the  start 
of  the  cycle,  the  epoxy  appears  to  have  become  stiffer  after  the 
first  cycle.  However,  by  always  referencing  the  creep  strain  to  the 
initial  sample  length,  the  material  softens  slightly.  This  becomes 
more  evident  in  looking  at  the  net  creep  compliance,  AD,  in  Figure 
47  where  it  is  evident  that  the  value  of  Dj  has  decreased.  In 
general,  the  value  of  was  found  to  decrease  by  about  15%  from 
the  first  to  the  second  cycle,  regardless  of  the  temperature  range. 
This  behavior  is  not  clearly  understood  at  this  time. 

The  previously  quoted  values  for  DQ  and  Dj  are  all  based  on 
the  second  cycle  data,  which  should  be  more  representative  of  the 
actual  material  in  the  glass/epoxy  composite  as  a result  of  the 
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fabrication  process  which  induces  stresses  during  cure. 

Matrix  Creep  Compliance 

The  creep  compliance  for  the  Shell  58-68R  epoxy  can  be  deter- 
mined for  analytical  purposes  by  using  the  DQ  values  taken  directly 
from  the  solid  line  in  Figure  40  and  the  values  of  D determined 
using  D1  ( 75)  - 0.u69(10-6)  and  the  specific  a-j-  values  defined  by 
the  solid  line  in  Figure  44  (or  Figure  45).  The  total  creep  com- 
pliance is  then  defined  by  Equation  (50)  with  n = 0.19. 

Later  on  we  shall  need  to  know  the  creep  compliances  for 
20,  75  and  140°F  in  order  to  evaluate  the  various  micromechanics 
theories.  Following  the  approach  outlined  above,  we  find  the 
values  of  the  epoxy  creep  compliances  as  shown  in  Table  11. 

Table  11.  Shell  58-68R  Epoxy  Creep  Compliances 

Temperature  Creep  Compl iance 

(°F)  D(t)  x 10'6  (psi-1) 

20  1 .732  + .021t‘ 19 

75  1 .865  + .069f19 

19 


140 


2.0?'.  + . 184t‘ 
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Uniaxial  Creep  and  Recovery  of  Unidirectional 
Glass/Epoxy  Laminae  and  Laminated  Composites 

Determination  of  Material  Constants  in  Power  Law  Representation 

The  power  law  representation,  Equation  (50),  has  been  found 
to  describe  the  behavior  of  several  glass/epoxy  anu  graphite  com- 
posites [9,  12]  and  will  be  used  here  to  model  the  S- 901  glass/- 
Shell  58-60R  epoxy  composite  system.  Moreover,  the  use  of  visco- 
elastic micromechanics  theory  implies  that  the  composite  obeys  the 
same  power  law,  i.e.,  t0,19,  at  least  in  the  absence  of  any  signif- 
icant crack  growth  and  fiber  deformation.  When  the  fiber  angle, 

9,  is  not  close  to  zero  degrees  the  glass/epoxy  exhibits  an  appre- 
ciable amount  of  time  dependence,  particularly  at  the  higher  temp- 
eratures and  stress  levels.  Figures  49  through  52  represent  the 
typical  strain  response  for  e = 45°  and  e = ±45°  at  two  of  the 
temperatures  investigated.  Since  the  transient  part  of  the  creep 
strain. rbeys  a po.».  law  in  time,  the  strain  does  not  actually 
level  off  in  time,  as  it  appears  to  do  in  these  figures. 

In  the  case  of  the  off-angle  specimens,  0°<  0 < 90°,  thj 
creep  compliance  measured  was  sjj.  For  6=0°  and  e = 90°,  Sn 
and  S22  were  measured  directly  for  the  respective  angles.  The  use 
of  a transversely  mounted  strain  gage  on  the  latter  fiber  angles 
also  gave  us  the  values  of  the  compliances  S12  and  S , respec- 
tively. The  power  law  creep  compliance  can  be  rewritten  as 


recovery  of  45°  glass/ 


recovery  of  *45°  elass/eDOxv  at  75JF 
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S(t)  - SQ  + s/  (135) 

where  the  notation  more  closely  represents  that  commonly  used  in 
the  composite  area,  and  where 

SQ  = S^O)  (136a) 

and 

AS  » S:tn  (136b) 

The  Sj  coefficient  is  normally  called  the  "creep  coefficient". 

By  solving  pairs  of  simultaneous  equations  using  the  creep  data, 
in  the  same  manner  used  for  the  epoxy  creep  tests,  it  is  possible 
to  determine  both  SQ  and  S.  with  the  assumption  that  n = 0.19. 

The  computer  program  was  modified  in  order  to  account  for  values 
of  n greater  than  0.19  when  crack  growth  appeared  to  be  signifi- 
cant. All  of  the  strain-time  data  was  first  plotted  on  semi-log 
paper  and  a smooth  curve  was  drawn  through  each  data  set.  The  new 
data  points,  taken  from  the  smooth  curve,  were  used  in  all  sub- 
sequent analysis.  The  computer  program  solved  the  simultaneous 
equations  for  each  data  set  and  various  values  of  n from  0.19  to 
0.59.  In  order  to  determine  the  best  fit  for  SQ,  and  n the 
program  output  was  set  up  to  give  the  correlation  coefficient  be- 
tween the  experimental  AS  defined  by 


AS 


exp 


037) 


and  the  theoretical  AS  given  by  Equation  (136b).  The  ratio  of 
ASexp/AS  was  also  printed  out  to  assist  in  the  evaluation  of  the 
fit. 

Linear  Viscoelastic  Creep  Compliances 

In  Section  III  it  was  noted  that  the  residual,  or  permanent, 
strain  was  rezeroed  just  prior  to  the  start  of  the  next  loading 
cycle,  and,  therefore,  the  resulting  creep  compliance  was  referred 
to  the  sample  length  existing  at  the  start  of  the  cycle.  For  the 
purpose  of  the  analysis  to  follow,  the  creep  compliance  will  always 
be  referred  to  the  total  strain  referenced  to  the  initial,  first 
cycle  length.  This  procedure  is  more  desireable  for  engineering 
applications  and  comparisons  with  nonlinear  theory. 

At  sufficiently  low  stress  levels  the  behavior  of  the  glass/- 
epoxy  composite  can  be  considered  linearly  viscoelastic.  The  stress 
range  of  linearity  depends  on  the  temperature  as  well  as  the  fiber 
angle.  The  linear  viscoelastic  creep  compliances  were  determined 
at  stress  levels  between  5C0  and  2500  psi , depending  on  the  fiber 
angle.  In  all  cases  the  value  of  n was  found  to  be  0.19  with  good 
agreement  in  the  correlation  coefficient  in  the  third  or  fourth 
decimal  point,  e.g.,  0.9995.  The  linear  viscoelastic  creep  compli- 
ances found  experimental ly  are  tabulated  in  Table  12. 

Examination  of  the  experimental  data  taken  from  various  ten- 
sile coupons  showed  that  there  is  some  variation  between  samples 
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Table  12.  Linear  Viscoelastic  Creep  Compliances 
for  S -901  Glass/Shell  58-68R  Epoxy 


Temperature 

Fiber  Angle 

Creep  Compliance, 

(°F) 

(e) 

(x  10“6 

psi"1) 

20 

0 

.1210  + 

,0003t* 19 

45 

.3021  + 

.0046t> 19 

90 

.3115  + 

. 0025t* 1 9 

±45 

.2876  + 

.0023t* 19 

75 

0 

.1213  + 

o 

o 

o 

c+ 

ID 

20 

.1990  + 

. 0038t* 1 9 

45 

.3260  + 

.01 24t* 1 9 

90 

.3315  + 

. 0081 t* 19 

0/90 

.1793  + 

. 0C1 5t* 1 9 

±30 

.2037  + 

.0037t*19 

±45 

.3061  + 

.0082t-19 

±60 

.3534  + 

,0103t*19 

±80 

.3416  + 

.0152t-19 

140 

0 

.1265  + 

. 0003t • 1 9 

20 

.1957  + 

. 0439t* 19 

45 

.3524  + 

.0503t*19 

90 

.3490  + 

. 1 468t  • 1 9 

0/90 

.1829  + 

.0096t-19 

±30 

.2042  + 

,0345t* 19 

±45 

.1881  + 

. 3032t • 19 

±60 

.3181  + 

. 1 786 1 * 1 Q 

i r<\ 


within  the  same  plate  as  well  as  plate-to-pl ate  variability. 

This  variability  Is  reflected  somewhat  by  the  differences  in  the 
fiber  contents  shown  in  Table  3 and  can  account  for  the  approxi- 
mately 3-7%  variation  in  initial  compliances.  Strain  gage  mis- 
alignment, small  deviations  from  the  desired  fiber  angle  and 
loading  misalignment  all  are  thought  to  contribute  to  the  differ- 
ences in  initial  compliance,  as  well  as  effecting  the  transient 
compliance  if  the  errors  are  high  enough. 

The  creep  behavior  of  the  0°  fiber  angle,  although  very  small 
(less  than  0.5%),  can  probably  be  attributed  to  the  straightening 
out  of  any  misaligned  fibers.  The  effect  is  larger  than  would 
normally  be  predicted  by  the  rule  of  mixtures,  viz.,  Equation  (114). 
Antans  and  Skudra  [142]  also  noted  this  behavior  in  a glass/epoxy 
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Figure  54.  Creep  compliance  for  different  stress 
levels  (6  - 20° , T = 75°F). 
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Figure  57.  Creep  compliance  for  different  stress 
levels  (e  = 0790°,  T = 75° F). 
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Figure  67.  Creep  compliance  for  different  stress 
levels  (e  = ±60°,  T - 140°F). 
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4500  ps 1 the  strains  within  the  90°  layers  are  considerably  below 
that  which  would  cause  failure  in  a 90°  unidirectional  sample 
(see  Figure  31).  In  Figure  32  this  is  the  region  below  the  “knee'' 
in  the  stress-strain  curve.  At  the  higher  stress  levels,  the 
90°  layers  have  failed  and  the  0°  layers  are  carrying  the  load. 

This  is  typical  of  the  behavior  in  the  barrel  region  of  the  motor 
case  where  the  layers  are  predominantly  0°/90°. 

Crack  arrest  is  not  viuibily  evident  during  the  first  loading 
cycles  for  any  of  the  unidirectional  composites  shown  in  Figures 
53  through  56.  However,  the  behavior  of  the  ±45°  laminated  com- 
posite (Figure  59)  at  the  two  highest  stress  levels  shows  a definite 
decrease  in  the  relative  rate  of  creep  after  a period  of  time. 

It  is  hypothesized  that  some  of  the  cracks  arrest  under  conditions 
where  they  (i)  move  toward  the  fiber  where  the  crack  tip  is  arrested 
and/or  (ii)  the  crack  moves  away  from  the  fiber/matrix  interface 
toward  an  adjacent  fiber,  with  arrest  occuring  when  it  reaches  this 
fiber.  The  laminated  composite  provides  additional  crack  arrest 
boundaries  in  terms  of  the  individual  layers  oriented  at  an  angle  of 
2e  with  respect  to  the  adjacent  layer.  This  layering  may  tend  to 
cause  the  cracks  to  develop  and  propagate  between  layers,  causing 
considerable  delamination  at  the  higher  stress  levels.  This  delam- 
ination is  clearly  visible  as  an  edge  effect  at  stresses  near  fail- 
ure when  the  top  layers  begin  to  curl  upwards  at  the  outer  edge  of 
the  tensile  coupons. 
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The  creep  compliances  which  were  measured  at  the  lowest  stress 
levels  are  summarized  in  Figures  68  through  71  for  the  unidirec- 
tional and  laminated  composites.  These  compliances  were  denoted  as 
the  linear  viscoelastic  creep  compliances  previously  listed  in 
Table  12;  however,  at  140°F  even  the  lowest  stress  level  tests 
conducted  on  the  45,  90,  +45  and  ±60°  samples  indicated  some  degree 
of  irreversible  damage.  Tests  at  lower  stress  levels  were  not 
possible  as  a result  of  the  lack  of  strain  gage  sensitivity  coupled 
with  the  effect  of  small  temperature  changes  during  the  test. 

With  the  small  changes  in  creep  strain  at  lower  stress  levels, 
even  minor  changes  in  temperature  (±1°F)  create  cyclic  readings  of 
a few  microstrain.  A stress  level  of  300  psi  was  considered  to  be 
the  minimum. 

The  effect  of  temperature  can  be  seen  more  readily  in 
Figures  72  and  73  which  show  the  creep  compliance  measured  at  a 
stress  level  of  3000  psi  for  both  the  45°  unidirectional  and  ±45° 
laminated  glass/epoxy  composite.  Both  figures  show  the  strong 
influence  of  temperature  on  the  material's  internal  viscosity  which 
affects  the  creep  coefficient  through  the  reduced  time,  or  aT  shift 
function.  The  vertical  offset  indicates  at  first  glance  that  the 
initial  compliance  may  be  temperature  dependent.  The  initial  com- 
pliances, S^O),  for  the  45,  90  and  ±45°  fiber  angles  are  shown  in 
Figure  74.  Since  the  behavior  of  the  initial  matrix  compliance  is 


that  of  a TCM,  it  is  not  unexpected  that  the  glass/epoxy  composite 
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Figure  69.  Creep  compliances  at  low  stress  levels  for 
verto^f  ber  angles  of  laming  ^posUes 
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Figure  72.  Effect  of  temperature  on  creep  compliance  of  glass/epoxy  compo; 
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Figure  73.  Effect  of  temperature  on  creep  compliance  of  class/epoxy  comnosi 
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exhibits  a similar  behavior.  The  glass/epoxy  composite,  therefore, 
is  a TCM-2  by  definition  that  the  initial  compliance  is  temperature 
dependent.  Similar  behavior  was  noted  for  a graphite/epoxy  [12] 
and  an  E-glass/epoxy  [9].  Unfortunately , many  investigators  do  not 
recognize  this  temperature  dependence  when  trying  to  separate  out 
the  elastic  and  transient  components.  This  results  in  a series  of 
discontinuities  in  the  experimental  data  when  one  attempts  to  form 
a master  curve  using  isothermal  data  from  several  temperatures  [62]. 

The  time-temperature  shift  factor,  a^,  was  determined  for  the 
glass/epoxy  by  using  the  method  of  the  creep  coefficient  ratios  for 
the  45°  unidirectional  ar.d  ±45°  laminated  composite.  These  part- 
icular fiber  angles  were  selected  primarily  because  it  was  felt 
that  the  crack  growth  had  essentially  ceased  by  the  third  cycle  at 
low  stress  levels.  Lesser  fiber  angles  were  expected  to  give  more 
scatter  since  the  creep  coefficient  was  smaller.  The  a^  factor  is 
compared  in  Figure  75  with  the  epoxy  resin  shift  factor.  An  act- 
ivation energy  of  37  Kcal/g-mole  was  found  for  the  glass/epoxy. 

This  corresponds  to  the  values  of  38  Kcal/g-mole  and  35  Kcal/g-mole 
found  for  an  E-glass/epoxy  [9]  and  a graphite/epoxy  [12]  earlier. 

The  value  is  physically  consistent  in  terms  of  the  fiber  stiffnesses 
in  the  three  systems.  It  should  be  noted  that  the  values  of 
determined  for  all  of  the  laminated  composites  and  the  90°  unidir- 
ectional composite  at  140°F  were  considerably  below  the  log 
aT-temperature  curve  for  the  epoxy. 
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Studies  conducted  on  the  .'.70°  motor  case  materials  were 
carried  out  only  at  /5°F;  however,  the  same  behavior  is  expected 
in  terms  of  material  response  at  the  higher  temperatures . The 
linear  creep  compliance  could  not  be  determined  for  this  material 
since  the  motor  case  had  previously  been  pressurized  to  a high 
level.  From  Figures  5 and  13  one  can  see  that  the  case  material 
has  undergone  some  degree  of  permanent  damage  already.  However, 
the  results  of  some  multiple  cycling  tests  carried  out  at  low  and 
high  stress  levels  will  be  discussed  in  the  next  subsection. 

Multiple  Cycling  Effects 

All  of  the  glass/epoxy  composite  samples  were  tested  for  sev- 
eral creep  and  recovery  cycles.  In  order  to  simulate  the  actual 
conditions  in  the  motor  case,  at  least  two,  and  as  many  as  three 
cycles  were  run  on  the  tensile  coupons  at  all  stress  levels  and 
temperatures.  Since  the  total  number  of  tests  was  well  into  the 
hundreds,  it  would  not  be  possible  to  present  the  results  from  all 
of  them.  The  principal  results  from  the  tests  conducted  at  75°F  and 
140°F  on  the  20,  45,  90  and  2-45°  glass/epoxy  composites  are  shown 
in  Figures  76  through  95.  These  results  are  also  limited  to  only  a 
few  representati ve  stress  levels. 

Several  general  comments  appear  to  be  in  order  upon  inspection 
of  the  results.  In  all  of  the  tests  the  initial  compliance  contin- 
uously increases  with  each  cycle,  with  a disproporti onate  change 
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Figure  80.  Effect  of  multiple  loading  on  the  creep  compliance 
of  the  glass/epoxy  unidirectional  laminae 
(8  * 20°,  T « 140°F,  a = 10,500  psi). 
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Figure  82.  Effect  of  multiple  loading  on  the  creep  compliance 
of  the  glass/epoxy  unidirectional  laminae 
(e  = 45°,  T = 75°F,  o = 6350  psi). 
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Figure  86-  Effect  of  multiple  loading  on  creep  compliance 
of  the  glass/epoxy  unidirectional  laminae 
(6  - 90® , T = 75°F,  a = 2015  psl). 


iklL..t-«Ii  uK.li  Jiutf  Mi  MlMlUl 


f rww  wuppHu  'lartni 


iUiitadii4lliit.u~.lia,.  jutkililMtJillUil  Ulrin  iillull  s 


Figure  89.  Effect  of  multiple  loading  on  creep  comp 
of  the  glass/epoxy  unidirectional  lamina 
(e  * 90°,  T = 1 40°F,  c * 1480  psi). 
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Figure  95.  Effect  of  multiple  loading  on  creep  compliance 
of  the  glass/epoxy  laminated  composite 
(e  = ±45°,  T = 140°F,  a = 7000  psi). 
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occuring  between  the  first  and  second  cycle.  The  Initial  compli- 
ance, when  plotted  as  a function  of  the  loading  cycle,  appeals  to 
asymptotically  level  off  after  a number  of  cycles.  At  the  high 
stress  levels  three  cycles  does  not  represent  an  equilibrium.  Even 
at  moderately  low  levels  it  appears  that  three  cycles  are  not 
enough  to  generate  data  from  which  one  can  truly  measure  the  effects 
of  stress  level  in  the  absence  of  crack  growth. 

Also  evident  is  the  change  in  the  net  creep  compliance  with 
each  cycle.  During  the  first  cycle  a considerable  amount  of  crack 
propagation  and  possible  fiber  realignment  occurs.  The  latter  is 
believed  to  occur  more  predominantly  at  the  elevated  temperatures 
where  the  matrix  viscosity  Is  considerably  reduced  and  fiber  re- 
arrangement is  feasible.  All  of  the  75°F  tests  obeyed  a power  law 
with  n = 0.19  regardless  of  fiber  angle  or  stress  level  for  all 
loading  cycles.  However,  the  value  of  the  creep  coefficient  gener- 
ally decreased  very  rapidly  from  the  first  cycle!  The  creep  coeffi- 
cient for  the  second  cycle  was  usually  about  50-80%  c ' the  value 
during  the  first  cycle,  depending  on  the  stress  level.  At  140°F 
the  glass/epoxy  obeyed  a power  law  in  time  where  the  exponent  n 
varied  from  cycle  to  cycle.  The  value  of  n for  tests  at  the  high- 
est stress  level  sometimes  reached  as  high  as  n = 0.59.  As  an 
example  of  this  behavior,  the  creep  compliances  determined  by  ana- 
lytically fitting  the  data  to  a power  law  of  the  form  in  Equation 
(135)  are  shown  in  Table  13  for  several  fiber  angles  and 


Table  13.  Creep  Compliances  for  $-901  Glass/- 
Shell  58-68R  Epoxy  Composite  at  1 40° F 


Fiber  Angle 

Stress 

Cycle 

Creep  Compliance,  sji 

(e) 

(psi) 

(x  10-6  psi"1) 

20 

4800 

1 

.2185  + . 0097t* 4 1 

2 

.2406  + . 01 05t • 39 

3 

.2531  + .0l23t-3> 

20 

8501 

1 

.2633  + , 0582t*  59 

2 

.7196  + . 0944t* 1,7 

20 

10,501 

1 

.2642  + . 0555t*  'J,} 

2 

.7422  + . 08731- 

45 

3000 

1 

.4102  + . 1 077 1 ♦ ^ 1 

2 

.6947  + . 0879 1 * 3 3 

45 

5000 

1 

.4419  + .0221 t- 113 

2 

.4997  + .0602t- 27 

3 

.5349  + . 0722 1 • 3 

45 

6501 

1 

.5264  + . 1 441 1*  57 

90 

1480 

1 

.4357  + . 0498t*  '* 5 

2 

.5318  + . 0649t*  3-c' 

stress  levels. 

It  can  be  seen  that  n decreases  rapidly  with  each  succeeding 
cycle,  with  the  creep  coefficient  generally  increasing.  The  over- 
all result  is  a decrease  in  the  net  creep  compliance  with  each 
subsequent  cycle.  Given  enough  loading  cycles  the  crack,  growth 
should  stabilize  such  that  one  could  measure  the  effects  of  stress 


and  temperature  alone  as  did  Lou  and  Schapery  [9].  It  was  noted 
in  their  work  that  as  many  as  ten  cycles  were  necessary  to  get 
repeatable  results,  l.e.,  the  same  strain  output  under  the  same 
load  for  every  creep  test  thereafter. 

Lou  and  Schapery  [9]  were  able  to  express  their  data  in  terms 
of  Equation  (68)  and  a normalized  octahedral  shear  stress  given  by 

tQ  s a sine[y  - (y-1  )sin2e]1/,z  (138) 

where  o is  the  applied  stress,  0 the  fiber  angle, and  y is  defined  In 
terms  of  the  effective  matrix  Poisson's  ratio, 

Y = — 5-  (139) 

1 - ve  + ve 

The  nonlinear  initial  compliance,  5^(0),  in  [9]  was  written  In 
terms  of  the  linear  viscoelastic  creep  compliance  and  a scalar  fac- 
tor, gQ,  as  defined  in  Equation  (68).  In  a similar  manner  the 
creep  coefficient  was  also  defined  using  the  linear  viscoelastic 
creep  coefficient  and  the  scalar  factor  (g1g2/aan).  The  latter 
factor  was  found  to  roughly  follow  a hyperbolic  sine  function  as 
in  Equation  (72b)  with  0 = xQ  and  am  = aQ  (Figure  96). 

An  attempt  was  made  to  model  the  early  loading  cycles  using 
this  approach;  however.  Figures  97  and  98  show  that  gQ  is  a strong 
function  of  fiber  angle  and  tq  does  not  collapse  the  curves  into  a 
single,  invariant  function  for  either  test  temperature.  The  creep 


octahedral  stress 


octahedral  stress 


•rai'npiain,-- 


233 


E 


coefficient  appears  to  fit  the  hyperbolic  sine  function  fairly 
well.  Unfortunately,  several  points,  especially  for  e = 90°, 
are  not  very  close  to  the  analytical  function.  The  creep  coeffic- 
ient appears  to  depend  primarily  on  the  shear  stress,  t,  existing 
parallel  to  the  fibers  (see  Figure  99).  As  might  be  expected,  this 
arises  from  the  matrix  contribution.  However,  this  does  not  explain 
the  e = 90°  behavior  where  the  shear  stress  is  zero.  The  initial 
creep  compliance  appears  to  be  controlled  primarily  by  the  stress 
normal  to  the  fibers  through  the  nonlinear  parameter,  gQ,  as  shown 
in  Figures  100  through  105.  The  behavior  of  the  two  high  stress 
level  points  for  the  20°  fiber  angle  shown  in  Figures  104  and  105 
is  not  understood  at  this  time.  We  shall  point  out  later  how  these 
curves  may  be  used  as  correction  factors  in  the  actual  motor  case 
analysis. 

The  effect  of  multiple  cycling  on  the  creep  compliance  of  the 
case  material  for  the  ±70°  fiber  angle  is  shown  in  Figures  106  and 
107  for  two  stress  levels.  Remember  that  these  materials  have  al- 
ready been  stressed  to  a high  level  during  the  actual  motor  hydro- 
testing. As  a result,  the  creep  compliances  for  both  stress  levels 
are  considerably  higher  than  would  be  predicted  by  linear  viscoelas- 
tic analysis.  However,  the  trends  from  cycle  to  cycle  as  well  as 
stress  level  follow  the  same  pattern  established  for  the  tensile 
coupons  cut  from  the  plates.  The  net  creep  compliances  are  shown 
in  Figures  108  and  109  and  are  found  to  obey  the  power  law 
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Figure  99.  Multiple  loading  effects  on  creep  coefficient  showing  dependence  on 
the  average  shear  stress,  t (T  = 75°F). 
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property,  gQ,  as  s function  of  the  stress 
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Figure  105.  Initial  nonlinear  property,  g0,  as  a function  of  the  stress 
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Figure  106 


Temperature  75°F 
Fiber  Angle  ±70° 
Stress  300  psi 


0 1 
LOG  TIME  (minutes) 

Effect  of  multiple  loading  on  the  creep 
compliance  of  the  case  material  (e  = ±7i 
T = 75°F,  a = 300  psi). 


Temperature  75  F 
Fiber  Angle  ±70° 
Stress  300  psl 


O Cycle  1 
• Cvcle  2 


LOG  TIME  (minutes) 

Figure  108.  Net  creep  compliance  for  case  material  sub- 
jected to  multiple  loading  (e  = ±70°, 

T * 75°F,  a = 300  psl). 
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Figure  109-  Net  creep  compliance  for  case  material 

subjected  to  multiple  loading  (e  = +70°, 
T = 75°F,  a = 1000  psi). 
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with  n s 0.19.  The  data  from  the  plate  material,  therefore,  yields 
comparable  results  to  that  seen  in  the  filament  wound  motor  case 
materials. 


Lamina  Angular  Dependence  and  Principal  Creep  Compliances 

The  angular  dependence  of  the  compliances  are  given  by 
Equation  (89a).  The  uniaxial  compliance,  Sn(e),  may  be  expressed 
In  terms  of  the  two  principal  compliances  Sn  and  S22  and  the  uni- 
axial compliance  S)i(ei)  where  9i  t 6.  By  eliminating  the  term 
containing  S66  and  solving  for  Sn(o)  we  find  that 


where 


■> 


and 


Sn(e)  = SnUj)  + 


i2  m2 


- 


2 ' 

-JSii 


iz mzm? 


’22 


i - cose  m = -sine 

= cosOj  mj  = -sinOj 


(140) 


(141) 


s»  ” V 

(142a) 

■ s90“ 

(142b) 
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Therefore,  In  order  to  predict  the  angular  dependence  of  the  creep 
compliance  from  experimental  data,  one  only  needs  experimental  com- 
pliances for  Sq0,  SgQo  and  one  more  angle,  say  s|j(e).  For  the 
present  analysis  the  measured  compliance  S^o  was  used  to  predict 
the  other  compliances.  For  this  special  case  Equation  (140)  with 
0 i - 45°  becomes 

S I ! ( o)  = 4n?mzS45o  * (x.l‘-«.2m2)S0O  + (m4-Jt.?m2)Sg0o  (143) 

The  initial  compliance  and  the  one-hour  creep  compliance  have 
been  calculated  from  Equation  (143)  for  different  fiber  angles  for 
the  three  test  temperatures  and  are  shown  in  Figures  110  through 
112.  The  values  for  Sg0 , Sgg0  arid  S^o  were  taken  from  the  linear 
viscoelastic  creep  compliances  shown  in  Table  12.  The  experiment- 
ally determined  compliances  for  the  last  cycle  of  the  low  stress 
level  tests  are  shown  on  the  figures  for  both  the  unidirectional 
and  angle-ply  glass/epoxy  composites.  The  predictions  agree  reason- 
ably well  with  the  experimental  points  at  20° F and  75°F.  The  one- 
hour  curves  for  the  20°F  and  140°F  predictions  were  also  derived 
by  applying  the  aT  factor  from  the  epoxy  data  to  the  75°F  creep 
coefficient  in  the  S^o  and  Sgg0  compliances.  The  20°F  prediction 
was  essentially  the  same;  however,  the  140°F  indicates  what  appears 
to  be  considerable  crack  growth  nonlinearity.  It  should  also  be 
realized  that  a value  of  rt  - 0.19  was  used  for  the  predictions  and 
it  was  noted  earlier  that  n > 0.19  for  most  of  the  140°F  data  as  a 


dependence  of  the  linear  viscoelastic  creep  compliance. 


dependence  of  the  linear  vfscoelas 
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result  of  the  apparent  crack  growth.  In  the  absence  of  the  crack 
growth,  it  appears  that  the  power  law  representation  adequately 
predicts  the  creep  compliance;,. 

The  principal  creep  compliances,  necessary  for  the  complete 
definition  of  the  fibrous  composite  under  conditions  of  plane 
stress,  are  given  by  Equation  (86a)  for  a transversely  isotropic 
material.  They  are  Sn,  S22,  S12  and  S6G.  The  values  of  Sn  and 
52?,  which  are  the  compliances  for  e = 0°  and  o = 90°,  respecti vely , 
have  already  been  given  in  Table  12  for  the  linear  viscoelastic 
region.  The  principal  creep  compliance  S12  is  given  by  Eouation 
(83a) 

■VO  ] 


-Vio 

S|2  = -----  = 

Eli  E23 


(83a) 


which  may  be  rewritten  us 


S 1 2 = = “V-.  i$ 


2 2 


(144) 


Goth  relationships  assume  that  the  compliance  matrix  is  symmetric, 
i.e.,  ji:.  - b;. , . Experimentally,  S t ■ may  be  found  from  the  trans- 
verse strain,  , measured  duriny  a creep  test  on  a specimen  with 
c = 0"  under  an  axial  stress,  e , 


12 


The  experimental  values  were  found  to  be  independent  of  time  and 


are  given  In  Table  14  along  with  the  values  of  v12  and  v21  measured 
during  the  creep  of  a 0°  and  90°  sample,  respectl vely. 


Table  14.  Measured  Properties  S12,  vj2  and  v2J 
for  S -901  Glass/Shell  58-60R  Epoxy 


Temperature 

(°F) 

$12 

(x  10"6)(ps1-1) 

\>i2 

v2l 

20 

-.0329 

.272 

.151 

75 

-.0330 

.273 

.096 

140 

-.0359 

.283 

.051 

The  values  of  S12  and  v12  are  reasonable;  however,  the  values 
of  \j21  , when  used  with  the  S22  compliance  in  Equation  (144)  give 
values  of  S2i  which  differ  significantly  from  S12.  It  is  felt  that 
most  of  the  error  is  due  to  the  small  transverse  strain  associated 
with  v21  measurements  and  possibly  microstructural  damage  at  the 
140°F  temperature.  The  average  value  of  S21  over  the  temperature 
range  is  within  15%  of  the  S12  values.  All  subsequent  analysis 
assumed  that  Si2  = S21  and  the  values  of  S12  given  in  Table  14 
were  used. 

The  last  principal  compliance,  S66,  can  be  calculated  trom  the 
transformation  relation.  Equation  (89a),  from  the  known  values  of 
$n,  $22 » $12  the  compliance  from  one  other  angle,  sji(e). 

Using  the  45°  creep  compliance,  $450,  yields  the  following 
relationship. 
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$66  ~ ^45°  * 2S j 2 - S 1 1 - S22  (146) 


The  experimental  creep  compliance,  S66,  was  calculated  using  the 
linear  viscoelastic  creep  compliance  for  S^go,  Sq„(Su)»  sgQo(S22) 
and  S12  (Table  14).  The  time  dependence  of  S66  is  shown  in 
Figure  113.  The  solid  curve  represents  the  predicted  behavior  if 
one  were  to  use  the  epoxy  ay  shift  factors  on  the  75°F  creep  co- 
efficient to  obtain  the  creep  coefficients  at  20°F  and  140°F.  The 
initial  compliance  remains  unchanged.  The  net  compliance, 

S66(t)  - S66(0),  shown  in  Figure  114  can  be  seen  to  obey  the  power 
law  in  time  with  n = 0.19  as  expected.  The  ay  sh'ft  factors  deter 
mined  by  a horizontal  shifting  of  ,iS66  along  the  time  scale  are 
given  in  the  figure. 

All  of  the  four  principal  compliances  have  been  determined 
experimentally,  with  S66  predicted  from  experimental  data  and  not 
measured  directly.  Our  attention  will  now  be  turned  toward  the 
prediction  of  the  principal  compliances  using  the  various  micro- 
mechanics theories. 


Comparison  with  Micromechanics  Theory 


ii 

SJ 


The  "rule  of  mixtures"  representation,  viz.,  Equations  (114) 
and  (117),  was  used  to  predict  the  major  Poisson's  ratio,  v12,  and 
longitudinal  modulus,  En-  For  the  glass/epoxy  composite-  studied 
here  the  last  terms  in  Equations  (113)  and  (116)  may  be  dropped 


Figure  113.  Temperature  dependence  of  principal  shear  compliance. 
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with  only  negligible  error  in  order  to  arrive  at  the  "rule  of  mix- 
tures" relationships.  This  error  is  on  the  order  of  0.1  and  2A% 
for  the  respective  equations.  The  fiber  properties,  E^  and 
given  by  Equation  (129)  were  considered  independent  of  time  and 
temperature.  A fiber  volume  fraction,  vf,  of  0.616  was  used  for 
the  analysis  (see  Table  3 for  experimental  values  of  various 
plates).  The  matrix  volume  fraction  was  0.384.  The  viscoelastic 
behavior  of  En  and  vJ2  was  readily  found  by  employing  the  quasi- 
elastic approach  given  by  Equations  (118a)  and  (118b).  The  matrix 
modulus  was  determined  using  the  quasi -elastic  approach 


Emlt>  ‘ roty 


(147) 


where  the  values  of  Dm(t)  are  given  in  Table  11. 

The  principal  creep  compliances,  Su  and  S12»  are  expressed  in 
terms  of  Elt  and  v12  through  Equation  (83a)  and  were  predicted  with 
the  use  of  a computer  for  several  decades  of  time.  Calculation  of 
the  other  two  principal  creep  compliances,  S22  and  S6e,  is  more 
involved.  Several  micromechanics  theories  were  considered;  however, 
the  Halpin-Tsai  Equations  given  by  Equations  (119)  and  (120)  for  the 
elastic  and  quasi-elastic  cases,  respectively,  are  relatively  pop- 
ular. The  Halpin-Tsai  relations  are  a semi -empirical  approach  which 
has  been  used  by  several  Investigators  [29,  33,  62],  The  adjust- 
ment factors,  and  cG»  are  selected  to  give  the  best  fit  with 
experimental  data.  The  factors  are  known  to  be  functions  of  fiber 
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geometry  and  spacing;  their  time-dependence.  If  any,  has  not  been 
Investigated  by  others. 

The  initial  principal  creep  compliances  which  were  measured 
are  summarized  in  Table  15  for  future  reference. 


Table  15.  Experimental,  Initial  Principal  Creep  Compliances 
for  S “901  Glass/Shell  58-68R  Epoxy 


Temperature 

(°F) 

Su  x 10-6 
(psi"1) 

S12  x 10'6 
( PS i “ 1 ) 

S22  x 10“6 
(psi"1) 

S66  x 10"6 
(psi-1) 

20 

.121 

-.0329 

.3115 

.8417 

75 

.121 

-.0330 

.3315 

.9175 

140 

.127 

-.0359 

.3490 

1 .0054 

In  order  to  determine  the  optimum  values  for  ^ and  cG,  Equations 
(147)  and  (120)  were  solved  for  several  values  of  the  adjustment 
factors  using  initial  properties.  Equation  (83a)  was  used  to  con- 
vert E22  and  G17  to  the  compliances  S22  and  S66.  The  analytical 
values  of  S22  and  S66  were  plotted  as  a function  of  the  adjustment 
factor.  It  was  found  that  = 3.14  gave  the  best  fit  over  the 

entire  temperature  range,  never  exceeding  an  error  of  ±4%. 

The  temperature  dependence  of  the  initial,  principal  creep 
compliances  is  shown  in  Figure  115  along  with  the  corresponding 
analytical  predictions.  It  Is  evident  that  the  ‘rule  of  mixtures" 
and  the  Halpin-Tsai  relations  give  excellent  results.  The  dotted 
1 ir.e,  shown  for  S22  and  S66,  represents  the  predicted  Halpin-Tsai 
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compliances  when  ^ = 2 and  = 1 are  used.  The  specific  values 
are  based  on  numerical  solutions  [2]  but  it  is  evident  that  they 
do  not  give  good  predictions  for  the  glass/epoxy  composite.  We 
shall  comment  later  concerning  the  use  of  these  specific  values. 

The  reduced  stiffnesses,  Q..,  are  related  to  the  S,  . by  Equation 

* w ' J 

(128)  and  are  shown  in  Figure  116.  Again,  the  agreement  with  the 
predicted  values  (solid  line)  is  good. 

The  time  dependence  of  the  principal  compliances  was  found  by 
using  a quasi-el astic  analysis  and  the  linear  viscoelastic  creep 
compliance  for  the  epoxy  matrix  (see  Table  11).  Data  over  several 
decades  of  time  was  obtained  by  using  the  linear  viscoelastic  power 
law  coefficients  and  extrapolation  to  times  beyond  the  actual  exper 
i mental  time.  The  "rule  of  mixtures",  Equations  (114)  and  (117), 
and  Halpin-Tsai  Equations  (120)  were  solved  for  several  decades  of 
time.  Equation  (83a)  was  used  to  convert  the  engineering  prop- 
erties, Eu,  etc.,  to  the  principal  creep  compliances,  S... 

' J 

Figures  117,  118  and  119  show  the  time  and  temperature  dependence 
of  the  predicted  (solid  and  dotted  lines)  and  experimental  (open 
circles)  compliances.  All  of  the  curves  are  based  on  the  last 
cycle  of  creep.  Except  for  140°F  where  a significant  amount  of 
crack  growth  occurs,  the  agreement  is  good  over  several  decades  of 
time.  At  140°F  the  strong  time  dependence  exhibited  by  the  S22 
compliance  is  evident.  The  applied  stress  is  normal  to  the  fibers 
and,  therefore,  does  not  provide  any  known  mechanism  for  crack 
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arrest  In  the  absence  of  any  shear  stress  component.  The  opening 
mode  crack  growth  Is  hypothesized  to  predominate  very  strongly. 
Further  evidence  of  crack  growth  as  a result  of  the  normal  stress 
can  be  seen  in  the  behavior  of  the  creep  compliance  for  the  e = 90° 
specimen.  Although  the  initial  compliance  appears  to  be  reasonable, 
the  experimentally  determined  creep  coefficient  appears  to  be  high. 
As  an  example,  from  Figure  75  we  can  determine  the  value  of  log  a^ 
at  140°F  to  be  about  -3.25.  Using  this  value  of  and  the  creep 
coefficient  at  75°F,  the  140°F  creep  coefficient  is  predicted  to 
be  0.0336  (10~6).  The  experimental  value  shown  in  Table  12,  for  the 
last  creep  cycle  at  a stress  level  of  475  psi  is  more  than  4 times 
greater  than  predicted.  At  140°F  crack  growth  is  significant  and 
the  stress  level  used  to  define  the  linear  viscoelastic  creep  com- 
pliance, namely  475  psi,  already  is  creating  a strong  nonlinearity. 

The  corresponding  reduced  stiffnesses,  Q.  .,  found  by  using 

' J 

Equation  (128)  are  shown  in  Figures  120,  121  and  122.  Again,  the 
agreement  is  good  except  at  140°F. 

It  was  noted  earlier  that  the  values  of  ^ = 2 and  = 1 are 
often  used  in  the  Halpln-Tsai  relations.  These  values  appear  to 
work  reasonably  well  when  the  ratio  cf  fiber  to  matrix  properties 
is  high,  e.g.,  graphite/epoxy  and  boron/epoxy.  A comparison  of  the 
relative  stiffness  for  these  two  materials  gives  a range  of  50  to 
120  typically.  Glass/epoxy  composites  generally  lie  in  the  range 
from  20-25  at  the  lower  end  of  the  spectrum.  Several  investigators 
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Figure  122.  Linear  viscoelastic  reduced  stiffness  of 
glass/epoxy  composite  (T  = HOnF). 
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[33,  62]  further  justify  the  use  of  these  specific  values  in  order 
to  determine  an  "effective"  resin  compliance.  This  procedure  in- 
volves the  solution  of  the  Halpin-Tsai  Equaiions  (119)  or  (120) 
in  terms  of  the  "effective"  resin  properties  Em  and  G^.  The  argu- 
ment for  this  approach  is  based  on  a consideration  that: 

1.  The  bulk  epoxy  resin  properties  are  generally 
not  available. 

2.  The  mechanical  properties  of  the  bulk  resin 
and  the  resin  in  the  composite  may  be  different 
because  of  curing  conditions  and  chemical  compo- 
si tion. 

3.  The  effect  of  stress  concentrations  around  the 
fibers  may  affect  the  resin  properties. 

Sims  and  Halpin  [62]  used  this  approach  on  an  E-glass/epoxy 
and  a graphi te/epoxy  composite  with  fiber  volume  fractions  on  the 
order  of  0.56.  A comparison  of  the  "effective"  resin  properties 
with  measured  bulk  resin  properties  showed  that  the  "effective" 
properties  were  within  ±10%.  A careful  examination  of  Figures  117, 
118  and  119  shows  that  the  time-dependence  of  the  experimental  data, 
especially  for  S66,  agrees  more  closely  with  the  Halpin-Tsai  re- 
lations when  one  uses  = 2 and  = 1 • 

The  primary  difference,  with  the  exception  of  S22  at  140°F, 
appears  to  be  in  the  value  of  the  Initial  compliances.  In  order  to 
check  this  behavior,  the  Halpin-Tsai  Equation  (120)  was  solved  for 
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the  "effective"  resin  compliance  using  c . - 2 and  r,^  - 1 . Only  the 
initial  compliance,  DQ  of  the  epoxy  resin  was  determined.  The 
solution  of  Equation  (120)  yields  two  real  roots,  one  positive  and 
the  other  negative.  Only  the  positive  value  of  Dq  was  considered 
on  the  basis  of  physical  reasons.  Figure  123  compares  the  experi- 
mental initial  compliances  with  the  "effective"  compliances  pre- 
dicted from  the  S22  and  SG6  principal  compliances.  .Note  that  there 
is  a considerable  difference  in  the  S22  and  SGG  predictions;  how- 
ever, the  S22  compliance  is  not  nearly  as  sensitive  to  changes  in 
the  matrix  Dq  as  Sr,c*  The  "effective"  resin  properties  indicate 
approximate ly  a twofold  increase  in  stiffness  from  the  bulk  material. 
This  appears  to  be  unreasonably  high  and  it  is  felt  that  the  fibers 
are  creating  most  of  the  stiffening  effect.  However,  a closer 
examination  of  the  theoretical  bounds  on  the  principal  compliances 
will  be  discussed  later  in  order  to  justify  this  hypothesis. 

The  time-dependence  of  522  and  S66  using  the  "effective"  SG6 
resin  compliances  and  the  Halpin-Tsai  relations  with  = 2 and 
.'G  = 1 are  shown  in  Figures  124,  125  and  126.  Remember  that  only 
the  value  of  Dq  was  changed,  retaining  the  original  D1  values  found 
experimentally.  The  time-dependence  of  S66  agrees  remarkably  well 
over  several  decades  of  time.  The  predicted  S22  response  is  in 
close  agreement  except  where  considerable  crack  growth  is  present 
at  the  1 40° F temperature. 

The  theoretical  upper  and  lower  bounds  for  S66  and  S22  have 


Figure  125.  Principal  compliances,  S22  and  S66,  predicted 
equations  using  the  SG6  "effective"  resin  com 
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been  established  by  Hashin  and  Rosen  [189].  The  "composite  cylin- 
der assemblage"  (CCA)  model  corresponds  to  the  upper  bounds  on 
compliance.  The  CCA  model  is  obtained  by  progressively  filling  out 
a cylindrical  specimen  with  composite  circular  cylinders.  Each 
cylinder  consists  of  a circular  fiber  and  concentric  matrix  shell 
having  the  same  fiber  volume  fraction  as  the  overall  composite  and 
is  filled  out  in  such  a manner  that  the  assemblage  is  statistically 
homogeneous  and  transversely  isotropic.  The  upper  and  lower  bounds 
for  the  shear  modulus  in  the  plane  parallel  to  the  fibers,  G12»  are 
given  by  (after  some  rearrangement) 


G12 (+) 


rGf(l-vm)  + G(l+vj1 
f v m m m' 

g;o+v r+"TTT-*vT 

L T nr  nr  nr  J 


(148a) 


Gu(~) 


’ Gf(l+Vf)  + Gm(l-vf) 

070 -vf)  + G-r+vf7 


(148b) 


where  the  subscripts  f and  m refer  to  the  fiber  and  matrix,  res- 
pectively, and  G and  v are  the  shear  modulus  and  volume  fraction, 
respectively. 

Determination  of  the  bounds  on  the  transverse  modulus,  E22 
(or  $22).  is  more  Involved  since  it  requires  the  evaluation  of  the 
effective  plane  strain  transverse  bulk  modulus,  k*  as  well  as 
bounds  on  the  transverse  shear  modulus,  G*.  The  effective  plane 
strain  transverse  bulk  modulus  is  given  by 
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t m 1 + vm 

kf  - k k + 6 
T m m m 


where  k,  the  bulk  modulus  of  the  fiber  (f)  and  matrix  (m)  is 
defined  by 


In  Equation  (150),  E is  the  tensile  modulus  and  v is  the  Poisson's 
ratio  of  the  individual  constituents. 

The  upper  and  lower  bounds  on  the  transverse  shear  modulus 


are  given  by 


Gt(+)  = Gf + p — \ 

LGm  ' Gf  + 2WGfT- 


G*(-) 


' 1 

-Gf  " 5n 


(k  +2G  jv 
v m nr  m 

20r+G"T 
m m nr  - 


The  upper  and  lower  bounds  on  the  transverse  modulus,  E22  , may  now 
be  expressed  in  terms  of  k*,  G*  and  the  major  Poisson's  ratio,  v,,, 
and  longitudinal  modulus,  Eu,  given  by  Equations  (117)  and  (114), 
respectively.  These  bounds  become 


En(+)  - 


4kt  Gt(+) 


+G*(+J 


‘t  V1 2/  1 


(152a) 


E22^  " k*  + 


4k*  SJ(-) 


'■t  vi£'tiiJ 


(152b) 


Using  Equation  (83a)  to  convert  E22  and  G12  to  S22  and  S6C- 
respectively,  the  temperature  dependence  of  the  initial  compliances, 
r.ow  using  the  experimental  Dq  and  not  the  "effective"  initial  com- 
pliance, is  compared  with  the  experimental  data  as  shown  in  Figures 
12/  and  128.  The  previous  predictions  for  the  Halpin-Tsai  rela- 
tions are  also  shown  for  reference.  It  should  be  noted  that  the 
upper  bounds  (+)  on  the  compliances  are  the  values  determined  from 
the  lower  bounds  (-)  on  the  moduli  due  to  the  inversion.  The  oppo- 
site is  true  for  the  lower  (-)  compliance  bounds.  The  upper  bound 
on  S66(+)  also  corresponds  to  using  = 1 in  the  Halpin-Tsai  rela- 
tionships. 

Time-dependence  may  be  represented  by  using  the  associated 
relationships  and  the  quasi-elastic  analysis  approach.  Figures 
129  through  131  show  that  the  response  of  S22  and  S6$  is  well  within 
the  upper  and  lower  bounds  with  the  exception  of  S22  at  140°F. 

Rather  than  force  the  Halpin-Tsai  relations  to  fit  the  experimental 
data  by  requiring  the  matrix  compliance  to  change  significantly, 
it  seems  more  reasonable  to  adjust  the  factor  ^ and  Cq  according 
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to  the  best  fit  over  the  desired  time  scale  of  engineering  Interest. 
The  experimental  data  falls  between  the  bounds,  suggesting  that 
upper  bounds  for  S66(cG  = 1)  does  not  represent  the  actual  model  for 
the  g'lass/epoxy  composite  within  the  theoretical  constraints  orig- 
inally imposed. 

Constant  Crosshead  Rate  Tests  on  Glass/Epoxy  Laminates 
Effect  of  Temperature 

The  effects  of  temperature  on  the  multiple  cycling  of  a 45  and 
±45°  glass/epoxy  composite  tensile  coupon  is  shown  in  Figures  132 
through  135.  For  these  figures  only,  the  stress-strain  curve  has 
been  drawn  for  each  specimen  with  the  strain  always  referenced  to 
the  sample  length  at  the  start  of  each  cycle,  i.e.,  the  strain  at 
the  beginning  of  each  cycle  is  rezeroed  (e  = 0).  The  net  result, 
as  we  noted  earlier,  causes  the  material  to  appear  to  become  stiffer 
with  each  succeeding  cycle.  This  effect  is  misleading  and  is  cited 
here  omy  to  show  the  trend  in  cycle-to-cycle  behavior.  In  Figures 
132  and  133  the  glass/epoxy  appears  to  become  nonlinear  as  a result 
of  irreversible  damage  at  a stress  level  near  4000  psi  (75°F). 
However,  at  140°F  it  appears  that  this  damage  occurs  at  extremely 
low  stress  levels  for  the  45°  unidirectional  composite  (Figure  134), 
approaching  zero  stress.  The  ±45°  laminated  composite,  possibly 
because  of  the  layering  effect  on  crack  arrest,  becomes  nonlinear 
near  2000  psi  (Figure  135). 
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STRAIN,  c {%) 

Figure  132.  Effect  of  repeated  tensile  loading  on  stress- 
strain  response  (o  = 45°,  T = 75°F) . 
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At  75°F  it  appears  that  the  degree  of  crack  growth  and  arrest 
are  on  the  same  order  of  magnitude  for  the  second  and  subsequent 
cycles  for  both  the  45  and  ±45°  composites.  This  is  supported  by 
the  relative  change  from  the  first  to  second  cycles  as  seen  in 
Figures  132  and  133,  along  with  the  similarity  in  stress  levels^ 
at  which  the  two  cycles  begin  to  differ,  viz.  4000  ps i . Earlier  it 
was  reported  in  Table  12  that  the  initial  creep  compliance  for  the 
±45°  laminated  composite  was  slightly  lower  (6,.5%)  than  the  45° 
unidirectional  composite.  The  differences  in  fiber  volume  fraction, 
although  in  the  correct  direction  of  the  experimental  results,  are 
only  on  the  order  of  1%  (TAMU  data)  to  3.5%  (SCI  data).  A compar- 
ison of  the  two  tensile  stress-strain  curves  also  shows  the  same 
trend  in  the  stiffness. 

The  layering  effect  seems  to  cause  the  material  to  become 
stiffer  at  both  75  and  140°F.  At  14C°F  crack  arrest  occurs  more 
rapidly  in  the  +45°  laminated  composite  which  displays  the  same 
behavior  in  the  second  and  third  cycle.  The  45°  unidirectional 
composite,  on  the  other  hand,  still  shows  significant  change  be- 
tween the  second  and  third  cycle.  The  effect  of  temperature  appears 
to  strongly  affect  the  crack  growth  behavior  in  the  unidirectional 
45°  composite  more  than  the  ±45°  laminated  composite. 

Effect  of  Multiple  Cycling  and  Stress  Level 


After  the  initial  set  of  experiments  were  completed,  it  was 
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decided  to  conduct  another  series  of  tests  on  45,  0/90  and  ±45° 
glass/epoxy  tensile  coupons  In  order  to  further  study  the  effects 
of  multiple  cycling  and  stress  level,  Including  the  unloading 
phase.  All  of  these  tests  were  conducted  at  75°F  using  a ramp 
loading  and  unloading  rate  of  0.02  inches  per  minute  followed  by  a 
recovery  period  (no  load)  equal  to  two  times  the  total  time  under 
load.  The  stress  histories  used  for  this  series  of  tests  are 
shown  in  Table  16. 

The  first  series  of  tests  were  carried  out  on  the  0°/90° 
laminated  composite.  The  results  of  these  tests  are  shown  in 
Figures  136,  137  and  138.  Only  the  first  cycle  is  shown  for  the 
stress  level  of  15,000  psi  since  there  was  no  visible  evidence  of 
irreversible  damage  present.  All  three  cycles  exhibited  nearly 
elastic  behavior  within  the  time  scale  of  the  experiment.  Consid- 
ering the  maximum  strain  level  induced  in  the  90°  layers,  one  would 
not  expect  any  significant  damage,  at  least  not  visible  damage 
in  this  type  of  test.  The  total  experimental  time  for  one  complete 
cycle  (neglecting  recovery)  was  less  than  three  minutes;  therefore, 
the  viscoelastic  effects  for  the  0°/90°  fiber  angle  are  negligible. 

The  next  two  cycles  at  the  42,000  psi  level  reflect  the  break- 
down of  the  90°  layers  at  a strain  level  close  to  0.3%.  Once  the 
layers  are  broken  the  succeeding  cycles  (Figure  13C)  show  a perm- 
anent strain  and  little  if  any  hysteresis.  This  condition  is  prob- 
ably represontati ve  of  the  behavior  in  the  barrel  region  of  the 
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Table  16.  Stress  History  of  Multiple  Cycling, 
Constant  Crosshead  Rate  Tests 


Fiber  Angle 

Sample 

Cycle 

Maximum 

(0) 

Number 

Number 

Stress  (psi) 

0°/90° 

090.11 

1,2,3 

15,000 

A, 5 

42,000 

45° 

45.13 

1,2,3 

8000 

4,5 

5000 

6 

8000 

7,8 

10,000 

9(broke) 

12, 0G0 

±45° 

AP45.18 

1,2, 3, 4 

11,500 

5(broke) 

10,500 

±45° 

AP45.11 

1,2 

5000 
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Figure  137.  Effect  of  cyclic  loading  on  C°/9C'>  tensile 

coupon;  - 42,000  psi  (Cycle  4). 
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Figure  138.  Effect  of  cyclic  loading  on  0°/90c  tensile 
coupon;  a = 42,000  psi  (Cycle  5). 
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motor  case  during  the  first  cycle  (hydrotest)  and  subsequent  cycles 
(motor  handling  and  firing). 

In  the  initial  test  series  on  45°  unidirectional  glass/epoxy 
laminates  the  maximum  stress  was  about  6000  psi.  This  series  was 
conducted  to  study  the  effects  of  successively  higher  stress  levels 
up  to  complete  fracture.  The  first  three  cycles  shown  in  Figures 
139  and  140  indicate  the  amount  of  irreversible  damage  as  a result 
of  crack  growth.  The  accumulation  of  permanent  strain  and  the  de- 
crease in  the  hysteresis  loop  is  clearly  shown.  Figures  141  and  142 
show  the  two  cycles  conducted  at  a stress  level  of  10,000  psi.  The 
relative  change  in  the  hysteresis  loop  appears  to  be  on  the  same 
order  of  maqnitude  as  the  cycles  conducted  at  the  lower  stress 
levels. 

The  last  test  series  was  run  on  two  +45°  laminated  composites 
using  two  different  stress  histories.  Sample  AP45.18  was  tested 
completely  at  a high  stress  level  without  cycling  at  lower  levels. 

A maximum  stress  of  11,500  psi  was  selected  for  this  sample.  The 
stress  history  for  sample  AP45.11  consisted  of  several  cycles  at 
5000  psi  followed  by  several  additional  cycles  at  increasingly  high- 
er stress  levels  until  fracture  occurred.  Figures  143,  144  and  145 
show  the  rather  dramatic  changes  which  occur  through  the  first  four 
cycles  of  loading.  It  is  evident  that  a disproportionate  amount 
of  damage  is  done  in  the  first  cycle.  This  fact  agrees  with  the 
results  shown  earlier  for  the  creep  and  recovery  tests. 
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Close  comparison  of  the  second  and  fourth  cycle  indicates  that  the 
shape  of  the  hysteresis  loop  is  almost  identical,  differing  only  by 
a slight  shift  along  the  strain  axis.  This  corresponds  primarily 
to  slight  changes  in  the  initial  compliance  with  the  net  creep  com- 
pliance remaining  essentially  constant.  Crack  growth,  therefore, 
appears  to  be  arrested  very  rapidly  in  the  ±45°  glass/epoxy  even  at 
extremely  high  stress  levels.  Figure  146  was  drawn  to  graphically 
show  the  effect  of  the  cycle-to-cycle  damage  by  superimposing  the 
first  and  fourth  cycle  stress-strai n curves. 

The  last  series  of  tests,  on  sample  AP45.11,  was  conducted  to 
show  the  effects  of  multiple  cycles  at  low  stress  levels  followed 
by  cycles  at  sequentially  higher  stress  levels  until  ultimate 
fracture.  Figure  147  shows  that  cyclic  response  up  to  5000  psi 
stress  levels  are  relatively  insensitive  to  cycle-to-cycle  damage. 
Figures  148  and  149  represent  the  first  and  third  cycle  at  a max- 
imum stress  of  12,000  psi.  However,  the  sample  has  already  been 
subjected  to  fourteen  previous  cycles  at  stress  levels  below  12,000 
psi  as  rioted  in  Table  16.  Let  us  now  compare  the  bahavior  of  the 
two  samples,  AP45.11  and  AP45.18.  Remember  that  one  sample  has 
seen  only  four  high  stress  levels  near  11,500  psi  while  the  other 
sample  has  seen  a total  of  seventeen  cycles,  all  at  or  above  5000 
psi.  The  last  cycles  for  the  two  samples,  conducted  at  or  near 
a maximum  stress  of  11,500  psi,  are  shown  superimposed  in  Figure  150. 
Except  for  slight  differences  in  maximum  stress  level  and  sample 
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Figure  146.  Superimposed  stress-strain  curves  (Figures  143 
and  145)  showing  cycle-to-cycle  damage  effects 
on  ±45°  glass/epcxy  tensile  coupon 
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Figure  149.  Effect  of  cyclic  loading  on  ±45e  tensile 


coupon;  - 12,000  psi  (Cycle  17). 
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Figure  150.  Damage  effects  as  a result  of  different  stress  histories. 

Open  circle  (o)  represents  multiple  cycles  at  sequent- 
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variability the  stress-strain  behavior  appears  to  be  a function  of 
the  maximum  stress  level  achieved  during  the  sample's  history.  This 
behavior  is  typical  of  that  reported  by  Farris  et  al.  [20,  128] 
for  propellant  materials  although  he  proposed  the  dependence  on 
strain  using  Lebesgue  norms,  Schapery  [12,  20]  has 
proposed  the  use  of  Lebesgue  norms  of  stress  rather  than  strain  to 
model  this  behavior.  The  functional  behavior  of  the  Lebesgue  norm 
is  such  that  it  appears  ideally  suited  for  the  representation  of  the 
response  seen  in  Figure  150;  namely,  that  it  is  more  influenced  by 
the  maximum  stress  level  than  the  entire  past  history. 

Four-Point  Beam  Bending  Tests 

The  previously  described  uniaxial  tensile  tests  were  conducted 
in  order  to  investigate  the  time  and  temperature  dependence  of  the 
glass/epoxy  composites.  The  studies  have  shown  that  the  glf>ss/epoxy 
materials  are  nonlinearly  viscoelastic,  in  part,  probably  due  to 
the  time-dependent  development  and  growth  of  microscopic  cracks. 

The  degree  of  nonlinearity  depends  on  many  factors  which  already 
have  been  cited,  namely,  stress  level,  temperature,  fiber  angle, 
etc.  However,  the  uniaxial  tensile  test  differs  significantly 
from  a bending  test  in  that  a strain  gradient  exists  across  the 
the  thickness  which  is  not  present  in  the  former  one  (neglecting 
a slight  amount  of  bending  due  to  grip  effects).  This  series  of 
tests  was  conducted  in  order  to  determine  the  effect  of  the  strain 
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gradient  on  material  nonlinearity. 

Creep  and  Recovery  Tests 

The  results  from  some  of  the  creep  and  recovery  tests  are 

s » 

shown  in  Figures  151  through  155.  The  data  is  plotted  as  usual  for 
the  previous  tests  of  this  type  except  that  strain,  e,  rather  than 
compliance  is  used.  Strictly  speaking,  we  cannot  assume  that  all  of 
the  layers  of  the  beam  are  within  the  linear  viscoelastic  range, 
particularly  in  the  presence  of  a significant  strain  gradient. 
Neglecting  this  situation  for  the  moment,  we  shall  study  some  of  the 
behavior  shown  by  the  beam  tests.  Remember  also  that  the  strain 
used  here  represents  the  output  from  the  strain  gage  on  the  tension 
side  of  the  beam. 

In  general,  the  beam  cycle-to-cycle  variation  behaves  in  the 
same  manner  as  the  tensile  coupons,  showing  a disproportionate 
change  between  the  first  and  second  loading  cycles.  It  has  been 
found  that  the  strain  can  be  represented  by  the  power  law  in  time 
with  n = 0.19  for  all  of  the  tests,  regardless  of  fiber  angle  or 
stress  level,  with  excellent  correlation  at  each  cycle. 

With  the  behavior  pattern  now  established,  let  us  assume  that 
the  stress,  a,  at  the  outer  fibers  of  the  beam  is  still  within  the 
linear  range.  Consequently,  we  can  use  the  flexure  formula 
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Figure  152.  Effect  of  multiple  cycling  on  the  creep  strain  for  the  45 
glass/epoxy  beam  for  a 55.5  in-lb  applied  moment  (T  = 75°F 
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154.  Effect  of  multiple  cycling  on  the  creep  strain  for  the  +70° 
case  material  beam  for  a 1.5  in-lb  applied  moment  (T  = 75°F). 
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where  M Is  the  applied  moment  (constant  for  the  creep  test),  I is 
the  moment  of  inertia  of  the  cross-section  about  the  neutral  axis 
and  C is  the  distance  from  the  neutral  axis  to  the  outer  fibers. 

We  are,  of  course,  assuming  linear  elastic  (viscoelastic)  behavior 
although  at  the  higher  stress  levels  this  will  not  be  true.  The 
creep  compliance  may  be  uetermined  in  the  same  manner  as  we  have 
in  the  previous  tests  for  the  constant  applied  stress.  Table  17 
shows  the  effective  creep  compliance  found  by  using  the  stress 
defined  by  Equation  (153). 

The  initial  creep  compliance  for  the  +45°  glass/epoxy  tested 
at  the  bending  moment  of  12.5  in-lbs  is  within  2 % of  the  value 
shown  in  Table  12  for  the  uniaxial  creep  and  recovery  tests.  The 
creep  coefficient  which  does  not  change  significantly  with  sub- 
sequent cycles  also  agrees  with  the  uniaxial  data.  At  the  higher 
stress  levels  (moments)  the  ±45°  glass/epoxy  exhibits  the  soft- 
ening effect  seen  in  the  tensile  coupons. 

Consider  the  ±20°  glass/epoxy  case  material  next.  Earlier  we 
predicted  the  angular  dependence  of  the  linear  viscoelastic  creep 
compliances  in  Figure  111.  From  the  predicted  curve  based  on  third 
cycle  data  we  find  the  initial  compliance  to  be  0.187  x 10"6(psi"J), 
which  agrees  reasonably  well  with  the  beam  values.  The  agreement 
is  within  4.5%  after  correcting  for  the  difference  in  the  fiber 
volume  fraction  using  the  experimental  data  in  Table  3.  The  third 
cycle  beam  data  was  corrected  by  simply  multiplying  by  the  ratio  of 
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the  fiber  volume  fractions,  l.e.,  0.635/0.616,  to  obtain  an  esti- 
mated initial  creep  compliance. 

The  ±70°  case  material  does  not  show  good  agreement  with  the 
initial  linear  creep  compliance  of  0.348  x 10"F'  ( psl “ 1 ) which  would 
be  predicted  from  Figure  111,  being  in  error  by  about  100%.  How- 
ever, the  strain  which  exists  at  the  outer  fibers  of  the  beam  is 
on  the  same  order  of  magnitude  as  the  earlier  tensile  creep  and 
recovery  tests  which  had  also  indicated  significant  softening. 

It  is  not  unreasonable  to  expect  severe  softening  after  one  con- 
siders the  prior  history  of  the  case  as  a result  of  hydrotesting. 

There  is  a considerable  number  of  surface  cracks  present  in 
the  case  materials  as  a result  of  the  earlier  hydrotesting  (see 
Figure  5b).  The  material  is,  therefore,  much  softer  in  the  dome 
| region  of  the  motor  case  where  the  fiber  angles  are  greater  ( ;20° 

to  ±70°)  than  In  the  barrel  section  (0°/90°). 

Another  interesting  feature  is  the  relative  growth  in  the 
initial  creep  compliance  for  the  ±70°  glass/epoxy  with  each  cycle 
while  the  creep  coefficient  remains  virtually  unchanged  from  cycle 
to  cycle.  The  normal  stress  is  considerably  higher  than  the  shear  *■ 

stress  for  the  ±70°  fiber  angle  and  would  tend  to  explain  this  [ 

behavior  in  terms  of  sudden  (rather  than  slow)  crack  growth.  |i 
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effects  of  multiple  cycling  under  controlled  displacement  condi- 
tions. The  test  has  been  observed  to  be  essentially  a constant 
strain  rate  test  for  the  sample  stiffnesses  encountered.  A plot  of 
the  outer  fiber  strain  (tension  side)  as  a function  of  time  indi- 
cates that  the  strain  rate  is  constant  (Figure  156).  The  effects 
of  multiple  cycles  and  applied  moment  are  shown  in  Figures  157 
through  168  for  fiber  angles  of  ±45.  0/90,  +20  (case)  and  ±70° 
(case). 

In  all  situations  the  strain  shown  is  the  out-DUt  from  the 
tension  side  of  the  beam.  Additional  tests  conducted  using  strain 
gages  on  both  sides  of  the  beam  (tension  and  compression)  indicated 
the  compression  strain  to  be  less  than  the  tension  strain  for  the 
0°/90°  specimen.  The  neutral  axis,  therefore,  shifts  slightly 
toward  the  compression  side  as  a result  of  the  softening  which 
occurs  under  tension  in  the  90°  layers.  This  effect,  however, 
appears  to  be  very  small  ar.d  is  not  considered  further  ir,  the  ana- 
lysis of  the  beam  data. 

Therefore,  by  assuming  that  the  strain  is  a linear  function  of 
the  distance  from  the  centroidal  axis  (plane  sections  remain  plane) 
and  that  the  effect  of  tension  and  compression  stresses  is  the 
same,  we  can  predict  the  linear  end  nonlinear  behavior  by  using 
the  data  from  constant  crosshead  rate  tensile  stress-strain  curves. 
The  moment-strain  response  for  the  ±45°  g^ss/epoxy  beam  was  pre- 
dicted by  using  the  constant  rate- to-fai lure  stress-strain  curve 
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Figure  161.  Effect  of  cyclic  loading  of  0n/90° 
glass/epoxy  beam  (Cycles  1 and  3). 
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Figure  166.  Effect  of  cyclic  loading  of  ±20°  beam  cut 
from  motor  case  material  (Cycles  5 and  6) 


Figure  167.  Effect  of  cyclic  loading  of  ±70°  beam  cut  from  motor  case  material 


material 


for  the  45°  unidirectional  glass/epoxy  composite  (Figure  30). 

The  0°/90°  beam  behavior  was  predicted  by  using  the  corres- 
ponding 0°  and  90°  unidirectional  stress- strain  curves  (Figures  28 
and  31,  respectively).  The  moment  was  determined  by  employing  the 
following  procedure: 

1.  Select  several  strain  levels,  e.g.,  1000,  2000, 

....  10,000  ye,  etc  , corresponding  to  the  max- 
imum outer  fiber  strain.  The  data  for  the  i45° 
and  0°/90°  specimens  indicates  the  output  from 
both  gages  (tension  and  compression  surfaces)  is 
relatively  close;  therefore,  only  the  tension 
strain  gage  output  is  used  here. 

2.  Using  the  assumption  of  a linear  strain  distri- 
bution, determine  the  strain  at  the  mid-point  of 
each  layer  and  the  corre;  ponding  moment  arm  meas- 
ured from  the  neutral  axis  of  the  beam  (assumed 
fixed  for  all  strain  levels). 

3.  Use  the  strains  establishec  for  the  individual 
layers  to  determine  the  average  stress  in  the 
layer  from  the  stress-strain  curve  representing 
that  particular  fiber  angle.  (Some  of  the  stress- 
strain  curves,  e.g.,  45°,  hud  to  be  extrapolated 
to  higher  strain  levels). 

The  stress  acting  on  a given  lviyer,  multiplied 
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times  the  area  of  the  layer  cross-section  and 
the  moment  arm  yields  the  incremental  moment 
for  that  layer. 

5.  Add  all  cf  the  incremental  moments  to  give  the 
total  moment  for  a specified  strain  level. 

The  linear  theory  is  shown  in  the  ±45°  and  0°/90°  beam  results. 
The  agreement  is  very  good  over  a Urge  portion  of  the  range.  In 
fact,  the  presence  of  the  strain  gradient  greatly  reduces  the  soft- 
ening effect.  In  Figure  160  it  can  be  seen  that  the  nonlinear 
theory  predicts  a significant  decrease  in  the  applied  moment.  The 
45°  data  used  to  predict  the  bending  behavior  results  in  conserva- 
tive predictions  of  moment.  Remember  that  we  have  already  noted 
significant  crack  growth  during  the  first  few  cycles  in  th.  45° 
unidirectional  compos Ue  specimen  with  crack  arrest  probably 
occurring  only  when  the  crack  reaches  another  fiber.  The  +45° 
tensile  specimen  exhibited  more  crack  arrest  during  the  first  few 
cycles  and  this  was  attributed  to  the  layering  effect. 

The  beam  appears  much  stiffer  than  nonlinear  theory  predicts 
(Figure  160).  In  the  beam  crack  growth  must  occur  at  the  outer 
fibers  where  the  stress  is  highest  and  then  work  inward  toward  the 
center.  The  cracks  are  immediately  retarded  by  the  interface  at 
the  second  layer  and  must  be  redi rented  before  proceeding.  However, 
the  cracks  have  now  moved  into  a region  of  lower  stresses  which 
reduces  the  growth  rate.  It  could  also  be  argued  from  a statistical 
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viewpoint,  that  there  Is  a much  larger  percentage  of  critical 
cracks  at  a given  stress  level  in  a tensile  specimen  than  a beam. 

The  cyclic  loading  of  the  motor  case  materials  appears  to  have 
reached  an  equilibrium  behavior.  The  hysteresis  loop  remains 
essentially  unchanged  from  cycle  to  cycle  except  at  the  highest 
moments  (Figure  168).  However,  by  using  Figure  111,  we  can  pre- 
dict the  linear  behavior  using  the  linear  viscoelastic  creep  compli- 
ance. The  material  exhibits  considerable  softening  in  the  initial 
compliance  (see  Table  17)  although  the  creep  coefficient  which  re- 
flects time-dependent  crack  growth  appears  to  remain  unchanged 
with  successive  cycles. 

Plate  Twist  Tests 

The  case  of  plate  bending  under  the  action  of  twisting  moments 
is  discussed  by  Timoshenko  and  Woi nowsky-Krieger  [183]  and 
Lekhnitski  [25].  The  specific  test  outlined  earlier  in  Section  III 
is  commonly  called  the  "plate  twi^t  test".  Several  authors  [184, 
191-193]  have  applied  the  test  to  anisotropic,  and  more  specific- 
ally, orthotropic  materials  such  as  plywood  [191]  and  fiber-rein- 
forced  composites  [192,  193],  Whitney  [184]  and  Tsai  [193]  have 
used  the  test  to  determine  the  elastic  material  properties  of 
glass/epoxy  and  boron/epoxy  composites  using  linear  theory.  The 
purpose  of  these  tests  in  this  particular  program  was  to  investi- 
gate the  effects  of  multi  axial i ty  on  the  behavior  of  the 
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glass/opoxy  composite. 

Creep  and  Recovery  Tests 

Since  only  one  plate  was  available  for  the  plate  twist  tests, 
low  load  level  creep  and  recovery  tests  were  conducted  first.  Two 
load  levels  were  used  as  a check  for  linearity.  The  average  strain 
along  the  two  diagonals  (using  the  absolute  value)  was  found  to  be 
very  small  for  the  load  levels  used.  Again,  the  strain  was  found 
to  fit  a power  law  with  n = 0.19.  Bryzgalin  [148,  194]  conducted 
similar  tests  on  a glass/epoxy  plate  and  found  the  power  law  be- 
havior. In  order  to  accurately  measure  the  effects  of  damage 
using  this  test  it  would  be  necessary  to  use  higher  loads,  consid- 
erably longer  time  periods  and  higher  temperatures.  At  the  load 
levels  which  we  tested  there  was  no  measured  amount  of  damage.  The 
surface  strains  were  below  600  ye,  which,  from  Figure  30,  can  be 
seen  to  be  well  within  the  linear  region.  Further  creep  testing 
was  discontinued  and  attention  was  shifted  to  conducting  the  con- 
stant rate  tests  to  higher  load  levels. 

Constant  Crosshead  Rate  Tests 

The  constant  crosshead  rate  tests  were  conducted  to  loads 
which  caused  large  deflections  at  the  plate  tips.  In  addition  to 
the  large  deflections  encountered,  the  plate  was  loaded  in  a manner 
different  from  that  which  is  normally  treated  by  the  analysis. 
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The  analysis  used  within  the  linear  region  assumes  that  the  plate 
Is  loaded  at  the  corner.  Experimentally  this  causes  localized 
deformations  at  the  tip  due  to  the  concentrated  load  at  the  edge. 
The  test  also  requires  excellent  balancing  between  the  four  load 
points  to  prevent  the  plate  from  slipping  off  the  edges.  It  was 
also  noted  that  the  tests  are  often  conducted  with  the  loads  con- 
centrated a short  distance  in  from  the  edge,  but  using  the  original 
analysis  [195],  Although  not  evident,  it  is  believed  that  this  is 
common  practice  because  of  the  experimental  difficulties. 

In  the  Appendix  the  plate  twist  test  is  studied  analytically 
by  using  an  energy  method  and  a technique  is  proposed  to  take  into 
account  the  corner  effect.  The  incorporation  of  large  deflection 
theory  in  the  energy  formulation  results  in  a nonlinear  load- 
deflection  relationship,  Equation  (A-13),  which  predicts  that  the 
load  depends  on  the  cube  of  the  deflection  in  the  nonlinear  range. 

Equation  (A-13)  represents  the  case  of  an  isotropic  material. 
The  load-deflection  data  for  the  aluminum  plate  is  plotted  on  a 
log-log  scale  in  Figure  169.  The  cubic  behavior  in  wQ  is  not 
clearly  evident  until  the  linear  contribution  is  subtracted  out  as 
in  Figure  170,  after  which  the  data  clearly  indicate  the  cubic 
behavior.  In  fact,  apart  from  the  differences  due  to  the  material 
properties,  the  glass/epoxy  plate  would  also  display  the  same  func- 
tional behavior  in  wQ.  This  cubic  variation  results  from  the  form 
of  assumed  displacement  function  as  long  as  the  material  properties 
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do  not  depend  on  strain. 

The  load-deflection  curve  for  the  aluminum  is  shown  in  Figure 
171.  The  linear  theory  is  given  by  the  first  term  in  Equation 
(A-13)  using  the  correction  for  the  corner  effect.  Linear  theory 
predicted  by  this  method  shows  good  agreement  up  to  deflections  on 
the  order  of  one-half  the  plate  thickness.  This  is  the  normal  limit 
for  linear  theory  In  plates  and  shells  [183]  work.  The  nonlinear 
theory  predicted  by  Equation  (A-13)  exhibits  considerably  more 
stiffness  than  noted  experimentally.  This  is  probably  the  result  of 
the  choice  of  the  displacement  functions,  Equation  (A-l).  It  would 
be  possible  to  obtain  better  agreement  by  selecting  more  suitable 
displacement  functions  compatible  with  loads  removed  from  the  cor- 
ners a finite  distance.  The  plate  also  exhibits  an  instability 
which  will  te  discussed  later.  This  behavior  creates  larger  de- 
flections for  a given  load  and  the  displacements  appear  to  reflect 
single  curvature.  Therefore,  the  assumed  displacement  functions 
are  valid  only  in  a very  small  range. 

The  load  deflection  curves  for  the  ±45°  glass/epoxy  are  shown 
in  Figures  172  through  174  for  increasing  load  levels.  The  linear 
theory  for  the  ±45°  fiber  angle  is  given  by  Tsai  [193]  as 


*o  = 4F3- 


(154) 


where  P is  the  applied  corner  load,  £ is  the  length  of  the  side, 
h is  the  plate  thickness  and  SQ  is  given  as 


-deflection  behavior  of  t45°  glass/epoxy 


342 


0 » +45°  SQ  = 2(S2?.-S12)  (155a) 

6 = -45°  SQ  = Z(Sn-S12).  055b) 

The  ±45°  glass/epoxy  plate  consisted  of  15  layers;  therefore, 
the  average  value  of  SG  was  used  in  Equation  (154).  The  elastic 
values  of  Sn,  S!2  and  S22  are  given  in  Table  15.  The  predicted 
linear  theory  is  shown  on  the  figures  after  correcting  for  the 
corner  effect  as  we  did  before.  Linear  theory  agrees  reasonably 
well  up  to  plate  deflections  on  the  order  of  one-half  the  plate 
thickness. 

At  the  higher  loads  the  load-deflection  curves  indicate  what 
appears  to  be  a softening  effect.  When  plotted  as  a load-strain 
curve,  as  in  Figures  175  through  177,  the  same  effect  is  seen  upon 
unloading.  However,  in  view  of  the  surface  strain  levels  which  are 
seen  to  be  less  than  1400  ue,  one  would  not  expect  significant  non- 
linear behavior  such  as  this.  By  plotting  the  load-deflection 
curves  for  cycles  1 and  6 on  log-log  paper  (Figure  178),  it  can  be 
seen  that  there  are  three  linear  regions.  This  behavior  has  been 
described  qualitatively  by  Foye  [196]  to  consist  of  a linear  range 
followed  by  the  nonlinear  region  which  possesses  an  unstable  branch. 
After  reaching  this  unstable  point,  the  load-deflection  relationship 
again  becomes  linear.  Closely  examine  the  edges  of  the  plates 
shown  in  Figures  179  and  180.  The  low  level  creep  and  recovery 
test  shows  that  the  plate  edge  is  essentially  straight. 
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Figure  175.  Load  versus  average  surface  strain  for  ±45°  glass/epoxy 


Figure  176.  Load  versus  average  surface 
(Cycle  3). 


Finure  176.  Load  versus  average  surface 
(Cycle  3). 


versus  average  surface  strain  for  ±45°  glass/epoxy 
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In  Figure  180  the  edge  shows  considerable  curvature  as  a result  of 
the  high  load  which  has  caused  the  plate  to  become  unstable  with 
respect  to  the  twisting  moments.  Also  noticeable  is  the  lack  of 
double  curvature  (saddle  effect)  after  the  plate  becomes  unstable. 

Once  instability  has  occurred  the  plate  should  be  treated 
differently  analytically.  The  plate  could  possibly  be  considered 
as  a triangular  plate  loaded  at  the  corner  as  a cantilever  beam. 

A single  parabolic  displacement  function  would  appear  to  be  a 
reasonable  choice. 
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SECTION  V 

APPLICATION  TO  MOTOR  CASE  DESIGN 

The  observations  and  results  found  In  the  previous  sections 
on  the  effects  of  stress,  temperature  and  multiple  cycling  on  the 
glass/epoxy  composite  strongly  suggest  a means  for  improved  design 
methodology.  In  this  section,  we  shall  first  review  the  present 
problem  associated  with  solid  rocket  motor  case  design  and  analysis 
and  existing  technology.  Lebesgue  norms  are  used  to  predict  the 
effects  of  stress  and  multiple  cycling  for  conditions  which  are 
typical  of  the  motor  case  as  a result  of  hydrotesting  and  subse- 
quent mission  loads. 

Review  of  Present  Technology 

The  present  design  and  analysis  of  fiber-reinforced  composite 
solid  rocket  motor  cases  is  firmly  entrenched  in  the  use  of  linear, 
orthotropic  elastic  analysis  coupled  with  the  application  of  "exper- 
ience" factors  to  account  for  effects  such  as  temperature,  case 
thickness,  hydrotesting,  etc.  [13-17].  It  should  be  noted,  however, 
that  in  the  absence  of  any  sound  theoretical  baseline,  experience 
carries  considerable  weight.  Under  certain  conditions,  such  as 
low  stress  levels,  short  time  periods  and  moderate  temperatures, 
a linear,  orthotropic  elastic  analysis  using  either  measured  or 
predicted  (micromechanics)  properties  may  be  adequate.  In  this 
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same  context  one  may  extend  the  approach  to  linear  viscoelasticity 
in  the  absence  of  significant  crack  growth. 

However,  Rawe  [197]  has  pointed  out  the  need  to  assess  the 
damage  created  during  hydrotesting  of  filament-wound  motor  cases 
such  as  that  shown  in  Figure  5b.  Rawe  observed  that  the  cracks 
become  visible  at  pressures  which  were  only  40%  of  the  required 
burst  pressure.  The  number  of  cracks  was  found  to  increase  with 
pressurization  and  became  a permanent  feature  of  the  chamber's 
appearance  upon  release  of  the  pressure.  He  also  noted  that  a dis- 
proportionate amount  of  cracking  occurred  during  the  first  (hydro- 
test) cycle.  In  addition  to  the  obvious  concern  for  the  softening 
effect  created  by  the  crack  growth,  the  presence  of  cracks  leads  to 
eventual  reduction  in  strength  as  a result  of  moisture  penetration 
during  prolonged  storage.  We  have  also  noted  the  effect  of  moisture 
on  the  volumetric  expansion  (see  Figure  26)  and  the  evidence 
suggests  [198]  that  moisture  is  present  in  the  composite  as  the 
result  of  small  microcracks. 

Crownover  [14]  suggests  that  a pressure  level  of  80%  of  the 
burst  pressure  causes  significant  damage  to  the  case  material.  His 
results  showed  that  the  strength  reduction  was  a function  of  both 
the  pressure  level  and  the  duration  of  the  holding  period.  Although 
his  observations  were  qualitative,  they  are  consistent  with  the  use 
of  Lebesgue  norms  to  describe  the  time-dependent  crack  growth 
effects . 
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In  the  Mlnuteman  III  motor  case  the  design  loads  are  nominally 
determined  for  a safety  factor  of  1 . ii5  [17].  The  hydrotest  pressure 
Is  5%  higher  than  the  mean  operating  pressure  and  80S&  of  the  min- 
imum burst  pressure.  For  design  considerations  the  factor  of 
safety  is  applied  to  the  mean  operating  pressure  and  not  the  hydro- 
test pressure.  The  properties  used  in  the  motor  analysis  are  based 
on  a fiber  volume  content  of  0.657  which  Is  slightly  higher  than 
that  reported  in  Table  3.  In  order  to  account  for  damage  due  to 
the  hydrotest  cycle,  secant  moduli  determined  from  strain  measure- 
ments on  the  surfaces  of  several  6~1nch  diameter  and  full  scale 
pressure  chambers  were  compared  by  ASPC  [17]  with  the  moduli  pre- 
dicted from  the  "rule  of  mixtures"  and  laminated  plate  theory. 

These  correction  factors  on  stiffness  were  found  to  be  on  the  order 
of  0.4-0. 7,  depending  on  the  load  (strain)  level  and  motor  location 
(fiber  angle  and  layer  sequence). 

Temperature  degradation  factors  were  determined  by  comparing 
the  flexural  stiffness  of  glass/epoxy  specimens  at  various  temper- 
atures. As  expected,  these  "temperature  degradation  factors" 
actually  reflect  the  viscoelastic  temperature  dependence  of  the 
composite  as  well  as  crack  growth  behavior  as  the  strain  level  in- 
creases. These  factors  on  stiffness  were  found  to  be  approximately 
0. 7-1.0  for  the  specified  operating  temperatures.  Certain  locations 
of  the  motor  case  actually  reach  temperatures  approaching  250-300°F 
for  short  periods  of  time.  It  is  entirely  possible  that  the 
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time-dependent  growth  of  cracks  during  this  period  is  more  signif- 
icant than  the  effect  of  hydrotesting.  However,  the  effects  of  the 
higher  temperatures  are  not  known. 
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The  present  approach,  therefore,  is  based  on  prior  motor  case 
experience  and  utilizes  an  empirical  method  for  correcting  temper- 
ature and  multiple  cycle  effects.  No  effort  has  been  made,  to  the 
author's  knowledge,  to  correct  for  more  than  one  loading  cycle 
(hydrotest)  or  the  time-dependent  effects  due  to  crack  growth. 

Prediction  of  Multiple  Cycling  Effects 

We  have  already  shown  that  the  linear  viscoelastic  properties 
can  be  predicted  reasonably  well  using  the  Halpin-Tsai  Equations 
(120)  and  the  "rule  of  mixtures"  Equation  (118)  in  the  absence  of 
crack  growth.  However,  we  also  noted  that  the  crack  growth  is 
significant  even  at  low  stress  levels  when  the  temperature  increases. 
Experimental  evidence  shows  that  there  is  considerable  softening  of 
the  off-angle  glass/epoxy  composite  which  must  be  considered  in  the 
dome  areas  of  the  motor  case  due  to  the  fiber  angles  employed. 
Multiple  cycling  effects  were  noted  in  Figures  100  through  105 
where  it  was  found  that  the  nonlinear  behavior  depended  primarily 
on  the  stress  normal  to  the  glass  fibers,  a . However,  the  nonlin- 
ear property,  g , reflects  stress  dependence  and,  as  shown  in 
the  figures,  does  not  explicitly  account  for  the  cycle-to-cycle 
damage. 
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In  Section  II  the  constitutive  relationships  were  given  for 
materials  exhibiting  damage  due  to  the  time-dependent  growth  of 
microcracks,  viz..  Equations  (102),  (106)  and  (107).  The  consti- 
tutive theory  is  based  on  fracture  mechanics  and  Involves  the 
Lebesgue  norm  functions  of  stress  given  by 


(156) 


i,  L. 

where  IMIp  is  normally  called  the  ptn  order.  Lebesgue  norm  or  Lp 
norm.  In  Equation  (105)  the  value  of  p was  defined  in  terms  of  the 
matrix  exponent,  n.  Fracture  mechanics  implies  the  use  of  the 
norm  when  the  crack  tip  velocity  is  proportional  to  the  square  of 
the  stress  intensity  factor  [199].  Both  norms  were  studied  in 
order  to  determine  which  one  fit  the  observed  data. 


Using  the  form  of  Equation  (106),  the  constitutive  relation- 
ship for  the  glass/epoxy  may  be  written  as 


e 


(157) 


where  S(t)  is  the  linear  viscoelastic  creep  compliance  for  the 
specific  fiber  angle,  0,  and  oe^  is  the  "effective  stress"  given  by 

uef  = o[l  + f(||o||)]  . 


(107b) 
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The  function  f is  dependent  upon  the  stress  level  through  Lebesuwe 
norms.  The  problem  involves  measuring  the  experimental  strains,  t, 
due  to  an  applied  constant  stress,  a,  and  then  determining  the 
function  f. 

The  form  of  Equation  (157)  may  be  inverted  in  the  usual  sense 
to  obtain  (after  using  the  quasi -elastic  approach). 
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(158) 


where  we  can  now  consider  the  measured  strains  as  input  strains  by 
fitting  the  experimental  data  to  a power  law 

t = + £jtni  (15S) 


where  > ^ and  i;  have  the  usual  definition  and  reflect  the  stress 

dependence  due  to  crack  growth.  Also,  m = n in  the  linear  visco-  j 

| 

elastic  region  but  rn  > 0.19  in  the  presence  of  crack  growth.  j 

The  reciprocal  of  $(t),  which  obeys  the  power  law  form,  is  j 

i 

found  to  be  approximated  by  ' - \ 
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using  the  form  of  S(t)  given  by  Equation  (135). 

Equation  (158)  was  solved  for  ogf.  using  the  experimental  data, 
Equation  (159),  as  the  strain  input  and  the  approximation  given  by 
Equation  (161)  to  obtain 
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by  neglecting  the  terms  of  order  t2n  and  higher  in  Equation  (161). 
(Note:  In  the  actual  reduction  of  the  data,  all  of  the  terms  shown 
inequation  (161)  were  used;  however,  the  difference  was  found  to  be 
insignificant).  The  function,  f,  or  more  exactly,  the  function, 
(1+f),  was  found  by  using  Equation  (107b). 

Data  from  the  20,  45  and  90°  fiber  angles  at  the  75°F  and  1 40° F 
temperatures  was  used  to  examine  the  multiple  cycle  effects.  The 
data  is  shown  in  Figures  181  through  185  as  a function  of  the  Lz 
norm.  The  l2  norm  was  found  to  describe  the  behavior  better  than 
the  norm  using  n - 0 19  defined  by  Equation  (105).  The  data  clearly 
indicated  the  softening  effect  in  f through  several  cycles.  How- 
ever, it  appears  that  the  1.2  norm  does  not  completely  separate  out 
the  effects  of  stress  level  which  appear  to  reflect  more  crack 
growth  effects  in  the  initial  loading.  Evidence  of  this  is  re- 
flected by  the  small  jumps  within  each  cycle. 

Examination  of  the  curves  suggests  that  they  may  be  approx- 
imated by  a straight  line  given  by 
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Figure  182.  Softening  function,  1 +■  f,  for  = 45’’  at  75°F . 


183.  Softening  function,  1 + f,  for  3 = 90°  at  75°F. 


Figure  184.  Softening  function,  1 + f,  for  • = 20°  at  140°F. 
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Fiqure  185.  Softening  function,  1 + f,  for  0 = 45°  at  140°F. 
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1 + f - A(o)  + B(o) | |o | | 2 (163) 

where  A (a)  and  B (a)  exhibit  stress  dependence  and  ||o|j2  Is  the  L2 
norm.  As  a result  of  the  earlier  studies  we  have  found  that  A(o) 
and  B(o)  depend  primarily  on  the  stress  normal  to  the  fiber,  o . 

The  functions  are  shown  in  Figures  186  through  189.  The  values 
for  A(c>n)  and  B(an)  are  given  in  Table  18. 


Table  18,  Softening  Parameters  for  S-901 
Glass/Shell  58-68R  Epoxy 


Temperature 

(°F) 

A(on) 

B(o„) 

75 

1 + 38.5(10-6)an 

6.41(10-i°)on 

140 

1 + 288.3(10"6)an 

222.0(l0'10)an 

At  low  stress  levels  the  value  of  f approaches  unity  as  pre- 
dicted by  linear  viscoelastic  theory  in  the  absence  of  crack  growth. 
The  curves  also  show  the  strong  dependence  on  temperature. 

At  75°F  the  softening  arises  mainly  from  the  A(on)  contribu- 
tion. The  slope  of  the  curves,  B(cn),  at  75°F  could  be  considered 
constant;  however,  there  is  not  enough  data  available  for  various 
fiber  angles  and  stress  levels  to  assess  this  assumption.  The  data 
at  140°F  suggests  that  the  B(an)  contribution  Is  considerably  more 
significant. 

The  softening  effect  which  occurs  in  the  dome  region  of  the 


Figure  186.  Normal  stress  dependence  of  the  softening  parameter. 


motor  case,  as  a result  of  the  stress  levels  and  fiber  angles, 
can  be  accounted  for  by  using  the  L2  norm  and  the  softening  para- 
meters given  by  Table  18  in  Equation  (163).  The  effects  of  multiple 
cycles  such  as  hydro  testing , motor  handling  and  subsequent  firing 
can  be  predicted  through  this  procedure.  It  is  interesting  to  note 
that  the  L2  norm  given  by 
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is  actually  the  root  mean  square  (RMS)  value  of  the  stress,  a, 
times  t This  implies  that  the  motor  case  analysis  can  incor- 
porate the  use  of  the  local  RMS  value  of  the  stress,  or  the  normal 
stress,  into  existing  analyses  to  determine  the  amount  of  softening 
with  each  loading  cycle.  Laboratory  data  from  tensile  coupons, 
such  as  used  here,  can  be  used  to  determine  the  functional  relation- 
ships between  the  desired  L norm  and  stress. 
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SECTION  VI 
CONCLUSIONS 

The  Shell  58-68R  epoxy  resin  was  found  to  exhibit  linear 
viscoelastic  behavior  for  the  stress  levels  studied,  which  cover 
the  range  from  5 to  50%  of  ultimate.  Isothermal  data  indicated 
the  epoxy  to  be  thermorheologically  simple  except  for  temperature 
dependence  of  the  initial  compliance  with  the  creep  compliance 
represented  by  a power  law  in  time  with  a constant  exponent 
n = 0.19. 

Evaluation  of  the  S-901  Glass/Shell  58-58R  epoxy  composite 
under  uniaxial  tensile  loading  showed  that  the  creep  compliance 
could  also  be  represented  by  a power  law  in  time  in  the  range  of 
linear  viscoelastic  behavior  with  the  same  exponent  n = njg. 
However,  the  glass/epoxy  composite  was  found  to  exhibit  a strong 
nonlinearity  with  both  stress  and  temperature.  Temperature  depen- 
dence analogous  to  that  of  the  epoxy  was  found.  At  the  higher 
temperatures  and  stress  levels,  the  nonlinearity  was  attributed 
primarily  to  crack  growth  in  addition  to  the  normal  reversible 
nonlinearity.  During  the  first  few  loading  cycles  the  exponent  in 
the  power  law  representation,  n,  was  found  to  be  considerably 
greater  than  0.19.  The  value  of  the  exponent  decreased  with  each 
subsequent  loading  cycle.  The  rate  of  decrease  depended  upon  the 
particular  stress  level  as  well  as  the  temperature.  The  nonlinear 
properties  were  found  to  depend  primarily  on  the  stress  normal  to 
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the  fiber  rather  than  the  shear  stress,  which  implies  agreement 
with  the  crack  opening  mode  of  growth  behavior;  for  reversible 
nonlinearity  dependence  on  the  octahedral  shear  stress  was 
found  [9]. 

The  time-dependence  of  the  principal  creep  compliances,  Sn, 
Si;.,  S.;,  and  S(lli,  was  determined  from  experimental  data  using 
fourth-order  tensor  transformations . Using  these  tensor  relations, 
the  angular  dependence  of  the  compliances  was  also  predicted  as  a 
function  of  fiber  angle  and  found  to  agree  with  the  experimental 
results  for  all  but  the  hignest  temperature  of  140°F;  at  tin's 
temperature  the  disagreement  appears  to  be  a result  of  significant 
crack  growth  and  possible  reversible  nonlinearity  even  at  low 
stress  levels. 

Crack  growth  during  the  first  cycle  seems  to  cause  a dispro- 
portionate amoun,.  of  damage  compared  to  subsequent  cycles.  This 
behavior  is  reflected  in  both  the  unidirectional  and  laminated  com- 
posites. However,  it  appears  that  there  is  a difference  in  the 
crack  growth  and  arrest  in  the  laminated  composites  in  tension. 

Both  the  creep  and  recovery  and  constant  crosshead  rate  tests  in- 
dicate that  the  angle-ply  ('!•)  composites  probably  further  reduce 
crack  growth  as  a result  of  the  barriers  created  by  the  layer  inter- 
faces. Consequently , the  angle-ply  composites  were  found  to  be 
stiffer  than  their  corresponding  unidirectional  counterparts  in 


tension. 
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The  motor  case  materials  were  also  found  to  obey  the  power  law 
creep  compliance  with  n " 0.19.  They  exhibited  the  same  character- 
istics as  the  composites  fabricated  in  plate  form.  Experimental 
data  indicate  that  the  off-angle  composites  exhibit  considerable 
softening  due  to  crack  growth.  This  behavior  will  affect  the 
stiffness  of  the  dome  regions  of  the  motor  case  appreciably  more 
than  the  barrel  section. 

The  "rule  of  mixtures"  and  Halpin-Tsai  relations,  Equations 
(118)  and  (120),  respectively,  were  found  to  describe  the  principal 
creep  compliances  in  the  linear  viscoelastic  range  reasonably  well 
in  the  absence  of  significant  crack  growth.  The  adjustment  factors 
were  found  to  be  approximately  equal  with  a value  of  = Cg  = 3.14. 
Effective  epoxy  resin  compliances  predicted  by  using  r = 2 and 
Cg  = 1 were  found  to  be  considerably  lower  than  the  bulk  epoxy 
compliances.  Determination  of  the  theoretical  upper  and  lower 
bounds  on  compliance  suggests  that  the  way  in  which  the  fibers  are 
distributed  in  the  cross-section  contributes  most  of  the  additional 
stiffness  to  tee  Halpin-Tsai  model  rather  than  the  in-situ  matrix 
as  previously  suggested  by  Sims  and  Halpin  [62].  The  experimental 
data  fall  within  the  established  bounds,  except  for  S22  at  NOT. 
Opening  mode  crack  growth  at  low  stresses  appears  to  exist  at  that 
temperature . 

Bending  tests  conducted  using  beams  and  plates  indicate  a 
strong  influence  of  the  strain  gradient  on  reducing  the  softening 


of  the  material.  The  beam  tests  show  that  linear  theory  agrees 
very  well  for  laminated  composites.  The  surface  cracks  are 
apparently  arrested  rapidly  by  the  layer  interfaces  and  concurrent 
movement  into  a region  of  lower  stress  (strain).  Nonlinear  beam 
theory,  based  on  the  uni di rectional  data  underpredicts  the  actual 
stiffness  cons iderably . The  bending  tests  imply  that  stretenmg, 
ratner  than  bending,  creates  significantly  more  softening  in  the 
compos i te. 

The  experimental  plate  data  show  good  agreement  with  linear 
theory  for  small  deflections.  A nonlinear  theory,  developed  using 
the  energy  method  for  large  deflections,  shows  that  the  load- 
deflection  relationship  exhibits  cubic  dependence  on  the  deflection. 
The  experimental  results  indicate  that  an  instability  occurs  in  the 
plate  twist  test  at  higher  loads  which  agrees  qualitatively  with 
Foye’s  results  [196].  Creep  and  recovery  bending  tests  on  the 
plate  in  the  linear  range  also  follow  the  power  law  behavior. 

Multiple  cycling  effects  on  the  glass/epoxy  composite  showed 
that  the  response  was  independent  of  earlier  lower  stresses. 
Viscoelastic  fracture  mechanics  theory,  together  with  the  observed 
behavior,  suggested  the  use  of  Lebesgue  norms  to  model  multiple 
cycling  effects.  The  second-order,  or  U norm,  was  used  to  repre- 
sent the  softening  effect  for  various  stress  levels  and  numerous 
loading  cycles.  Higher  order,  L|;..,  norms,  were  examined  but  the 
best  agreement,  was  obtaii.ed  using  the  L-,  norm.  The  softening 
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effect  was  found  to  depend  primarily  on  the  stress  normal  to  the 
fiber.  The  L2  norm,  which  is  equal  to  the  root  mean  square  (RMS) 
of  the  stress  multiplied  by  the  square  root  of  time,  may  be  used 
to  determine  the  softening  at  various  locations  within  the  motor 
case.  The  use  of  the  L2  norm  to  model  time-dependent  crack  growth 
softening  during  the  hydrotesting  and  subsequent  motor  case  loads 
is,  therefore,  proposed. 

Finally,  it  was  found  that  several  experimental  factors  should 
be  considered  in  terms  of  refining  the  determination  of  the  non- 
linear effects.  The  errors  resulting  from  gage  misalignment, 
assumed  fiber  orientation  (due  to  sample  fabrication),  and  general 
sample- to-sample  variability  may  easily  be  corrected  by  normaliza- 
tion. Thus  to  achieve  greater  correlation  with  the  data,  all 
specimens  should  first  be  tested  at  a low  stress  level  (without 
crack  growth  if  possible)  in  order  to  determine  the  linear  visco- 
elastic compliance  for  the  specific  specimen.  The  nonlinearity 
resulting  from  stress  dependence  in  the  absence  of  cracks  must  be 
determined  after  crack  growth  has  substantially  arrested.  Non- 
linearity due  to  crack  growth  can  be  determined  more  accurately 
by  subtracting  out  the  effects  due  to  stress  dependence  on  the 
viscosity.  Healing  effects  were  not  investigated  but  previous 
work  [9]  showed  that  they  may  be  significant  for  long  rest 
periods. 
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APPENDIX 


Analysis  of  Plate  Twist  Test 

In  this  Appendix  the  energy  method  is  used  to  obtain  a large 
deflection  solution  for  the  plate  twist  test  [191-193]  for  the 
case  of  an  isotropic  material.  The  analysis  is  then  corrected  for 
the  experimental  application  of  loads  at  points  removed  from  the 
corners  of  the  plate.  Finally,  extension  of  the  approach  to  include 
orthotropic  materials  is  indicated. 

Consider  a square  plate  with  the  coordinate  system  (x,  y)  as 
shown  in  Figure  A-l  where  a is  one-half  the  length  of  the  diagonal. 

A pure  twisting  moment  is  imposed  on  the  plate  by  loading  all  four 
corners  with  equal  forces.  The  forces  are  perpendicular  to  the 
plate  with  those  forces  at  the  first  and  third  (diagonal)  corners 
being  downward  and  the  other  two  forces  upward.  The  corner  loads 
cause  the  square  plate  to  assume  a hyperbolic  paraboloid  of  saddle- 
shaped  surface  [183,  184]. 

Displacement  functions  are  assumed  to  take  the  form 


w 


w = -f  (x2-y2) 

(A-l a) 

a 

u = Ay2x  + Cx 

(A-l b) 

v = -Ax2y  + Cy  . 

(A-l c) 
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The  deflection,  w,  follows  the  form  for  linear  theory  given  by 
Timoshenko  and  Woinowsky-Krieger  [183]  and  reduces  to  the  tip  de- 
flection, w , when  x = ±a  or  y - ±a.  The  displacements  u and  v are 
even  functions  of  y and  x,  respectively,  as  a result  of  symmetry. 
The  coefficients  A and  C will  be  found  by  minimizing  the  potential 
energy  [183], 

The  equilibrium  equation  for  small  deflections  is  given  by 

[183] 


S4W 


34W 


+ 2 -&L-  + 

L p 0 L 

ax  ax  ay  ay 


0 . 


' A-2 ) 


The  assumed  displacement  functions  automatically  satisfy  equilibrium 
for  small  deflections.  For  large  deflections,  in-place  forces  N , 

A 

N and  N are  developed  which,  if  exact,  will  satisfy  the  following 

y ™ 

equilibrium  equations, 


aNv  3N 

_2i+  0 

ax  ay 


(A- 3a) 


aN  aN 

0 . 

ax  ay 


(A-3b) 


The  strain-displacement  relationships  for  large  deflections 


must  include  the  strain  in  the  middle  surface  of  the  plate  during 
bending.  The  strain  components  are 
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= 3_U  1 / 3w\2 

ax  2 yaxy 

_ 3_V  1 /3w\2 

" ay  ?\ay/ 


- + av  + i w aw 

Yxy  'iy  3x  3x3y 


(A-4a) 

(A-4b) 

(A-4c) 


For  an  isotropic  material  the  constitutive  aquation  is  given  by 


ex  = H^V'V 


(A-5a) 


cy  = h^Ny“vNx] 


1 


Yxy  = FIT  Nxy 


( A- 5b ) 


(A-5c) 


where  h is  the  plate  thickness,  E is  the  elastic  modulus  and  v 
is  Poisson's  ratio. 

The  strain  energy  arises  from  the  pure  bending  contribution 
and  the  additional  energy  which  results  from  the  stretching  of  the 
middle  surface.  The  bending  strain  energy  is  given  by 


'B 


= § (\(~  + 3jj)2  ~ 2(l-v)  [ ^ ^ - I—)2]  i dA  (A-6) 

CX  (\3x2  ay  / L3x  ay  Vaxay/J) 


where  D is  the  flexural  rigidity  given  by 

Ds 

12(1 -v2) 


(A-7) 


The  strain  energy  due  to  stretching  is  given  by  the  expression 


V 


S 


2vExGy 


+ ^->*2xy] 


dA 


(A-8) 


The  total  strain  energy  is  the  sum  of  and  V<-.  The  area 
integration  may  be  performed  over  one  quadrant  to  obtain  the  strain 
energy  in  terms  of  the  three  unknown  parameters,  wQ,  A and  C,  by 
applying  Equations  (A-l),  (A-4),  (A-6)  and  (A-8).  The  total  poten- 
tial energy  of  the  system  is 

U * VB  + V$  - Pwo  (A-9) 


where  the  last  term  represents  the  potential  energy  due  to  the 
applied  load  acting  at  the  corner  with  the  tip  deflection  wQ.  The 
coefficients  may  now  be  determined  by  minimizing  the  total  energy 
with  respect  to  the  coefficients 


A » 


all 

3* 

= 0 

(A-lOa) 

all 

aC 

= 0 

A-lOb) 

aU 

3wo 

= 0 . 

(A-lOc) 

are 

found  to 

be 

V 

" 7 - llv* 

(A-l 1 ) 

a1*” 

1 5 - 11 v 

22w„2  r . 

C f- 

3a2  Ll 5 - llvj 


(A-12) 


Equation  (A-lOc)  yields  the  load-deflection  relationship  after  some 
rearrangement 


Eh3  w + 4Eh(43-t-33v) 

3a2(l+v) J 0 45a2(l+v) (15-11 


The  first  term  represents  the  linear  contribution  while  the  second 
term,  cubic  in  wQl1s  the  membrane  or  stretching  contribution. 

Recall  that  Equation  (A-13)  gives  the  corner  load  and  tip 
deflection  (wQ)  when  the  plate  is  loaded  at  the  corner.  Using 
Equation  (A-la)  we  can  write  the  deflection  w1  at  a distance  a' 


along  the  diagonal  as 


- ■ (f)\ 


Now  let  us  examine  the  last  term  in  Equation  (A-9)  more  closely. 
We  may  rewrite  this  term  using  Equation  (A-14)  to  obtain 


*0  * p(r-)2“' 


which  gives  the  same  energy  contribution  if  the  load  acting  at  a1 


were  given  by 


■ ■ i-)2  • 


(A-16) 


Consequently,  Equation  (A-14)  may  be  rewritten  in  terms  of  the 
experimental  load,  P ' , and  corresponding  deflection,  w',  acting  at 
point  a'  on  the  plate's  diagonal, 


P' 


* — — 1 w 1 x + [ — l(w'x)3  (A-17) 

_3aZ(l+u) J L45a2(l+v)(15-llv)  . 


where 


A * (fr)*  • (*-18) 

Following  this  technique  we  can  approximate  the  load  P*  and  corres- 
ponding deflection  wQ‘  without  having  to  assume  a more  complicated 
set  of  displacement  functions  to  account  for  the  contribution  of  the 
region  past  the  loads. 

Note  that  the  procedure  for  transferring  the  point  of  load 


application  is  completely  independent  of  the  type  of  material, 
differing  only  in  the  coefficients  in  the  load-deflection  terms. 
Con  segue n t ly,  we  can  apply  the  procedure  equ ally  as  well  to  aniso- 
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